1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
|
function [z]=RndIntT(n)
//------------------------------- OK
[lhs,rhs]=argn(0)
if rhs <= 0 ; n=10000;end
imin=-10;
imax=10;
y=RndInt(1,n,-10,10);
i=imin-2:imax+2;
z=[];
for i1=i, z=[z,prod(size(find(y==i1)))],end
plot2d3("onn",i',z'/n,[1,2],"151","Simulation ",[-12,0,12,0.1]);
i1=(imin:imax)';
plot2d1("onn",i1,ones(i1)/prod(size(imin:imax)),[-2,5],"100","Theory");
function [z]=RndDiscT(n)
//------------------------------- OK
[lhs,rhs]=argn(0)
if rhs <= 0 ; n=10000;end
x=[1,3,4,6,10,12];
pr=[0.1,0.2,0.3,0.2,0.1,0.1];
y=RndDisc(1,n,x,pr);
i=0:13
z=[];
for i1=i, z=[z,prod(size(find(y==i1)))],end
plot2d3("onn",i',z'/n,[1,3],"151","Simulation ",[0,0,14,0.5]);
plot2d1("onn",x',pr',[-2,6],"100","Theory");
function [zt]=BinomialT(n)
//------------------------------- OK
[lhs,rhs]=argn(0)
if rhs <= 0 ; n=10000;end
prb=0.5;
N=10;
y=Binomial(1,n,prb,N);
i=0:10;
z=[];
for i1=i, z=[z,prod(size(find(y==i1)))],end
plot2d3("onn",i',z'/n,[1,3],"161","Simulation");
deff('[y]=fact(n)','y=prod(1:n)');
deff('[z]=C(N,n)','z= fact(N)/(fact(n)*fact(N-n))');
i=0:N;
zt=[];
for j=i, zt=[zt, C(N,j)*prb^j*(1-prb)^(N-j)];end
plot2d1("onn",i',zt',[-2,6],"100","Theory");
function []=GeomT(n)
//------------------------------- OK
[lhs,rhs]=argn(0)
if rhs <= 0 ; n=10000;end
pr=0.2
y=Geom(1,n,pr)
N=20
i=0:N;
z=[];
for i1=i, z=[z,prod(size(find(y==i1)))],end
plot2d3("onn",i',z'/n,[1,3],"161","Simulation");
zt=[0];for i1=1:N; zt=[zt,pr*(1-pr)^(i1-1)];end
plot2d1("onn",i',zt',[-2,6],"100","Theory");
function [z]=PoissonT(n)
//------------------------------- OK
[lhs,rhs]=argn(0)
if rhs <= 0 ; n=1000;end
pmean=3;
y=Poisson(1,n,pmean);
N=20;
i=0:N;
z=[];
for i1=i, z=[z,prod(size(find(y==i1)))],end
plot2d3("onn",i',z'/n,1,"161");
deff('[y]=fact(n)','if n==0 then y=1;else y=n*fact(n-1);end');
zt=[];for i1=0:N; zt=[zt, exp(-pmean) *pmean^i1/fact(i1)];end
plot2d1("onn",i',zt',[-2,6],"100","Theory");
function []=ExpT(n)
//------------------------------- OK
// lambda exp(-lambda x) x>=0
// ---------------------------
[lhs,rhs]=argn(0)
if rhs <= 0 ; n=1000;end
lambda=3;
y=Exp(1,n,lambda);
histplot([0:0.1:10],y,[1,1],'051',' ',[0,0,5,3]);
deff('[y]=f(x)','y=lambda*exp(-lambda*x);');
x=[0:0.1:10]';plot2d(x,f(x),1,"000");
titre= 'macro histplot : Histogram plot';
xtitle(titre,'Classes','N(C)/Nmax');
function []=WeibullT(n)
//-------------------------------
[lhs,rhs]=argn(0)
if rhs <= 0 ; n=10000;end
y=Weibull(1,n,1,2)
histplot(20,y,[1,1],'061');
function []=HyperGeomT(n)
//-------------------------------
[lhs,rhs]=argn(0)
if rhs <= 0 ; n=10000;end
y=HyperGeom(1,n,1,10)
histplot([0:0.25:10],y,[1,1],'061',' ',[0,0,10,0.4]);
function [y]=ErlangT(n)
//-------------------------------
[lhs,rhs]=argn(0)
if rhs <= 0 ; n=10000;end
y=Erlang(1,n,10,1)
histplot(20,y,[1,1],'061');
|
