Rev 4330 | Rev 10075 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed
Rev 4330 | Rev 9638 | ||
---|---|---|---|
Line 5... | Line 5... | ||
5 | slib_parms=2\ |
5 | slib_parms=2\ |
6 | ,[data]\ |
6 | ,[data]\ |
7 | 1,p |
7 | 1,p |
8 | slib_author=Sophie, Lemaire; Bernadette, Perrin-Riou |
8 | slib_author=Sophie, Lemaire; Bernadette, Perrin-Riou |
9 | 9 | ||
10 | slib_example=[2,4,7],6. |
10 | slib_example=[2,4,7],6.1\ |
- | 11 | 0,1\ |
|
- | 12 | [20,21,22,23,24,25,50],16 |
|
11 | slib_require=pari |
13 | slib_require=pari |
12 | !exit |
14 | !exit |
13 | 15 | ||
14 | :proc |
16 | :proc |
15 | 17 | ||
Line 21... | Line 23... | ||
21 | !exit |
23 | !exit |
22 | !endif |
24 | !endif |
23 | 25 | ||
24 | slib_data=!declosing $slib_data |
26 | slib_data=!declosing $slib_data |
25 | slib_cnt=!itemcnt $slib_data |
27 | slib_cnt=!itemcnt $slib_data |
- | 28 | pari_precision=100 |
|
26 |
|
29 | !!slib_out=!exec pari slib_data=[$slib_data];\ |
27 | 30 | !!print(vector($slib_cnt, x,if(slib_data[x] >= 0, incgam(floor(slib_data[x])+1,$slib_p)/gamma(floor(slib_data[x])+1)))) |
|
28 | 31 | ||
- | 32 | slib_out=!exec octave poisscdf([$slib_data],$slib_p) |
|
- | 33 | ||
- | 34 | slib_out=!words2items $slib_out |
|
29 | slib_out=!trim $slib_out |
35 | slib_out=!trim $slib_out |
- | 36 | slib_out=!values x for x in $slib_out |
|
- | 37 |