Subversion Repositories wimsdev

Rev

Rev 4330 | Rev 7692 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 4330 Rev 4500
Line 7... Line 7...
7 
5,n\
 7 
5,n\
8 
0.5,p
 8 
0.5,p
9 
 
 9 
 
10 
slib_author=Sophie, Lemaire; Bernadette, Perrin-Riou
 10 
slib_author=Sophie, Lemaire; Bernadette, Perrin-Riou
11 
 
 11 
 
12 
slib_example=[2,4,5],20,0.7
 12 
slib_example=[0,1,2,3],3,7/15\
- 
 
 13 
[2,4,5],20,0.7
13 
slib_require=pari
 14 
slib_require=pari
14 
!exit
 15 
!exit
15 
 
 16 
 
16 
:proc
 17 
:proc
17 
 
 18 
 
Line 26... Line 27...
26 
!endif
 27 
!endif
27 
slib_out=
 28 
slib_out=
28 
slib_data=!declosing $slib_data
 29 
slib_data=!declosing $slib_data
29 
slib_cnt= !itemcnt $slib_data
 30 
slib_cnt= !itemcnt $slib_data
30 
slib_data1=!values rint(x)  for x in $slib_data
 31 
slib_data1=!values rint(x)  for x in $slib_data
31 
 
 - 
 
32 
slib_n=$[rint($slib_n)]
 32 
slib_n=$[rint($slib_n)]
- 
 
 33 
slib_p=$[$slib_p]
33 
slib_p1=$[1-$slib_p ]
 34 
slib_p1=$[1-$slib_p ]
34 
slib_out=!exec pari {slib_data=[$slib_data1];\
 35 
slib_out=!exec pari {slib_data=[$slib_data1];\
35 
print(vector($slib_cnt,i,binomial($slib_n,slib_data[i])\
 36 
print(vector($slib_cnt,i,binomial($slib_n,slib_data[i])\
36 
*$slib_p^(slib_data[i])*$slib_p1^($slib_n-slib_data[i])));}
 37 
*$slib_p^(slib_data[i])*$slib_p1^($slib_n-slib_data[i])));}
37 
 
 38