tests reorganization :
[scilab.git] / scilab / modules / statistics / tests / unit_tests / cdfgam.tst
1 // =============================================================================
2 // Tests for cdfgam() function
3 //
4 // Scilab Team
5 // Copyright INRIA
6 //
7 // =============================================================================
8
9 prec  = 1.e-5;
10
11 shape = 2;
12 scale = 3;
13 bn    = 1;
14
15 deff('[y]=Gamma(x)','y=bn*(x^(shape-1) * exp(-scale*x))');
16
17 bn     = intg(0,10^3,Gamma);
18 bn     = 1/bn;
19 x      = 1:10;
20 [P,Q]  = cdfgam("PQ",x,shape*ones(x),scale*ones(x));
21 P1     = [];
22
23 for xx=x
24         P1=[P1,intg(0,xx,Gamma)];
25 end
26
27 if norm(P1-P) > prec then pause,end
28
29 shape    = 2*x;
30 scale    = 3*x;
31 [P,Q]    = cdfgam("PQ",x,shape,scale);
32 [x1]     = cdfgam("X",shape,scale,P,Q);
33 [shape1] = cdfgam("Shape",scale,P,Q,x);
34 [scale1] = cdfgam("Scale",P,Q,x,shape);
35
36 if norm(shape1-shape) > 10*prec then pause,end
37 if norm(scale1-scale) > 10*prec then pause,end
38 if norm(x1-x)         > prec    then pause,end