Subversion Repositories wimsdev

Rev

Rev 10075 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 10075 Rev 13586
Line 1... Line 1...
1 
!if $wims_read_parm!=slib_header
 1 
!if $wims_read_parm!=slib_header
2 
!goto proc
 2 
  !goto proc
3 
!endif
 3 
!endif
4 
 
 4 
 
5 
slib_parms=2\
 5 
slib_parms=2\
6 
,[data]\
 6 
,[data]\
7 
1,p
 7 
1,p
Line 16... Line 16...
16 
:proc
 16 
:proc
17 
 
 17 
 
18 
!distribute  item $wims_read_parm into slib_data,slib_p
 18 
!distribute  item $wims_read_parm into slib_data,slib_p
19 
 
 19 
 
20 
!default slib_N=1000
 20 
!default slib_N=1000
21 
!if  $slib_p<=0
 21 
!if $slib_p<=0
22 
slib_out=0
 22 
  slib_out=0
23 
!exit
 23 
  !exit
24 
!endif
 24 
!endif
25 
 
 25 
 
26 
slib_data=!declosing $slib_data
 26 
slib_data=!declosing $slib_data
27 
slib_cnt=!itemcnt $slib_data
 27 
slib_cnt=!itemcnt $slib_data
28 
pari_precision=100
 28 
pari_precision=100
29 
!!slib_out=!exec pari slib_data=[$slib_data];\
 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))))
 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))))
31 
 
 31 
 
32 
slib_out=!exec octave poisscdf([$slib_data],$slib_p)
 32 
slib_out=!exec octave poisscdf([$slib_data],$slib_p)
33 
 
 33 
 
34 
slib_out=!words2items $slib_out
 34 
slib_out=!words2items $slib_out
35 
slib_out=!trim $slib_out
 35 
slib_out=!trim $slib_out
36 
slib_out=!values x for x in $slib_out
 36 
slib_out=!values x for x in $slib_out
37 
 
 -