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