Rev 18195 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed
Rev 18195 | Rev 18466 | ||
---|---|---|---|
Line 26... | Line 26... | ||
26 | slib_out=!makelist $slib_n for x=1 to $slib_m |
26 | slib_out=!makelist $slib_n for x=1 to $slib_m |
27 | !exit |
27 | !exit |
28 | !endif |
28 | !endif |
29 | slib_u=!random 0,1 repeat $slib_m |
29 | slib_u=!random 0,1 repeat $slib_m |
30 | slib_C=!values x for x=0 to $[$slib_n] |
30 | slib_C=!values x for x=0 to $[$slib_n] |
31 | slib_out=!exec pari {slib_g(p,x,m)=local(L,u,q,i,k) ; |
31 | slib_out=!exec pari {slib_g(p,x,m)=local(L,u,q,i,k) ;L=List([]);\ |
32 | for(k=1,m,u=[$slib_u][k];q=p[1];i=1;\ |
32 | for(k=1,m,u=[$slib_u][k];q=p[1];i=1;\ |
33 | if(u<q,listput(L,x[1]),while(u>=q,i=i+1;q+=p[i]);listput(L,x[i])));L};\ |
33 | if(u<q,listput(L,x[1]),while(u>=q,i=i+1;q+=p[i]);listput(L,x[i])));L};\ |
34 | {slib_bino(n,p,m)=local(q,b,prob); |
34 | {slib_bino(n,p,m)=local(q,b,prob);b=List([]);\ |
35 | q=p/(1-p);prob=(1-p)^n;listput(b,prob);\ |
35 | q=p/(1-p);prob=(1-p)^n;listput(b,prob);\ |
36 | for(i=1,n,prob=prob*(n-i+1)/i*q;listput(b,prob));\ |
36 | for(i=1,n,prob=prob*(n-i+1)/i*q;listput(b,prob));\ |
37 | Vec(slib_g(b,[$slib_C],m))};\ |
37 | Vec(slib_g(b,[$slib_C],m))};\ |
38 | slib_bino($slib_n,$slib_p,$slib_m) |
38 | slib_bino($slib_n,$slib_p,$slib_m) |
39 | 39 |