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