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