Rev 4161 | Rev 7692 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed
Rev 4161 | Rev 4330 | ||
---|---|---|---|
Line 6... | Line 6... | ||
6 | ,[data]\ |
6 | ,[data]\ |
7 | ,N\ |
7 | ,N\ |
8 | ,K\ |
8 | ,K\ |
9 | ,n |
9 | ,n |
10 | 10 | ||
11 |
|
11 | slib_author=Sophie, Lemaire; Bernadette, Perrin-Riou |
12 | 12 | ||
13 | slib_example=[2,4,7],10,5,7 |
13 | slib_example=[2,4,7],10,5,7 |
14 | slib_require=pari |
14 | slib_require=pari |
15 | !exit |
15 | !exit |
16 | 16 | ||
17 | :proc |
17 | :proc |
18 | 18 | ||
19 | slib_out=$empty |
19 | slib_out=$empty |
20 | !distribute |
20 | !distribute item $wims_read_parm into slib_data,slib_N,slib_K,slib_n |
21 | 21 | ||
22 | slib_b= $[binomial($slib_N,$slib_n)^(-1)] |
22 | slib_b= $[binomial($slib_N,$slib_n)^(-1)] |
23 | slib_min= $[max(0,$slib_n+$slib_K-$slib_N)] |
23 | slib_min= $[max(0,$slib_n+$slib_K-$slib_N)] |
24 | slib_max= $[min($slib_n,$slib_N)] |
24 | slib_max= $[min($slib_n,$slib_N)] |
25 | !if $slib_max<$slib_min |
25 | !if $slib_max<$slib_min |