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=[ |
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 |