Subversion Repositories wimsdev

Rev

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.1
 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 
slib_out=!exec pari slib_data=[$slib_data];\
 29 
!!slib_out=!exec pari slib_data=[$slib_data];\
27 
print(vector($slib_cnt, x,if(slib_data[x] >= 0, incgam(floor(slib_data[x])+1,$slib_p)/gamma(floor(slib_data[x])+1))))
 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