Add license header (tests)
[scilab.git] / scilab / modules / statistics / tests / unit_tests / cdfgam.tst
1 // =============================================================================
2 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
3 // Copyright (C) ????-2008 - INRIA
4 //
5 //  This file is distributed under the same license as the Scilab package.
6 // =============================================================================
7
8 // =============================================================================
9 // Tests for cdfgam() function
10 // =============================================================================
11
12 prec  = 1.e-5;
13
14 shape = 2;
15 scale = 3;
16 bn    = 1;
17
18 deff('[y]=Gamma(x)','y=bn*(x^(shape-1) * exp(-scale*x))');
19
20 bn     = intg(0,10^3,Gamma);
21 bn     = 1/bn;
22 x      = 1:10;
23 [P,Q]  = cdfgam("PQ",x,shape*ones(x),scale*ones(x));
24 P1     = [];
25
26 for xx=x
27         P1=[P1,intg(0,xx,Gamma)];
28 end
29
30 if norm(P1-P) > prec then pause,end
31
32 shape    = 2*x;
33 scale    = 3*x;
34 [P,Q]    = cdfgam("PQ",x,shape,scale);
35 [x1]     = cdfgam("X",shape,scale,P,Q);
36 [shape1] = cdfgam("Shape",scale,P,Q,x);
37 [scale1] = cdfgam("Scale",P,Q,x,shape);
38
39 if norm(shape1-shape) > 10*prec then pause,end
40 if norm(scale1-scale) > 10*prec then pause,end
41 if norm(x1-x)         > prec    then pause,end