\r probas.gp

ploth(x=-4,4,statnormalpdf(0,1,x));
m = 0; s = 1;
sample = vector(10^5,i,statnormalrandom(m,s));
data = vecsort(sample);
ploth(x=data[1]-1,data[#data]+1,[statempiricalcdf(data,x), statnormalcdf(m,s,x)]);
ploth(p=1.e-4,1-1.e-4,statnormalinvcdf(m,s,p));

lam = 3.2;
sample = vector(10^5,i,statpoissrandom(lam));
data = vecsort(sample);
ploth(x=data[1]-1,data[#data]+1,[statempiricalcdf(data,x), statpoisscdf(lam,x)]);
ploth(p=0,1-1.e-4,statpoissinvcdf(lam,p));

df = 4;
ploth(x=0,20,statchi2pdf(df,x));
sample = vector(10^5,i,statchi2random(df));
data = vecsort(sample);
ploth(x=data[1]-1,data[#data]+1,[statempiricalcdf(data,x), statchi2cdf(df,x)]);

s = 1.; c = 0.8;
sample = vector(10^4,i,statzgaussrandom(s,c));
data = vecsort(sample);
ploth(x=data[1]-1,data[#data]+1,[statempiricalcdf(data,x), statzgausscdf(s,c,x)]);

s = 5.; c = 0.8;
sample = vector(10^5,i,statzgaussrandom(s,c));
data = vecsort(sample);
ploth(x=data[1]-1,data[#data]+1,[statempiricalcdf(data,x), statzgausscdf(s,c,x)]);

s = 30.; c = 0.8;
sample = vector(10^5,i,statzgaussrandom(s,c));
data = vecsort(sample);
ploth(x=data[1]-1,data[#data]+1,[statempiricalcdf(data,x), statzgausscdf(s,c,x)]);

d = 3;
for(i=1,10, v = randominball(d); print(v,"  \t", norml2(v)));

randomhnf(5,144)

nf = nfinit(x^2+3);
pr = nfrandomprime(nf,20); [pr.p,pr.f,pr.e]

randomperm(10)
permorder(randomperm(100))