Subversion Repositories wimsdev

Rev

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) ;L=listcreate(m);\
 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);b=listcreate(n+1);\
 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