Subversion Repositories wimsdev

Rev

Rev 10075 | Details | Compare with Previous | Last modification | View Log | RSS feed

Rev Author Line No. Line
20 reyssat 1 
!if $wims_read_parm!=slib_header
13586 bpr 2 
  !goto proc
20 reyssat 3 
!endif
4161 bpr 4 
 
20 reyssat 5 
slib_parms=2\
2073 bpr 6 
,[data]\
7 
1,p
4330 bpr 8 
slib_author=Sophie, Lemaire; Bernadette, Perrin-Riou
4161 bpr 9 
 
9638 bpr 10 
slib_example=[2,4,7],6.1\
11 
0,1\
12 
[20,21,22,23,24,25,50],16
10075 bpr 13 
slib_require=octave
20 reyssat 14 
!exit
15 
 
16 
:proc
17 
 
18 
!distribute  item $wims_read_parm into slib_data,slib_p
19 
 
20 
!default slib_N=1000
13586 bpr 21 
!if $slib_p<=0
22 
  slib_out=0
23 
  !exit
20 reyssat 24 
!endif
25 
 
26 
slib_data=!declosing $slib_data
27 
slib_cnt=!itemcnt $slib_data
9638 bpr 28 
pari_precision=100
29 
!!slib_out=!exec pari slib_data=[$slib_data];\
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))))
20 reyssat 31 
 
9638 bpr 32 
slib_out=!exec octave poisscdf([$slib_data],$slib_p)
33 
 
34 
slib_out=!words2items $slib_out
20 reyssat 35 
slib_out=!trim $slib_out
9638 bpr 36 
slib_out=!values x for x in $slib_out