Rev 7692 | Blame | Compare with Previous | Last modification | View Log | RSS feed
!if $wims_read_parm!=slib_header
!goto proc
!endif
slib_parms=3\
,[data]\
,a\
,b
slib_author=Sophie, Lemaire; Bernadette, Perrin-Riou
slib_example=[2,4,7],0.5,4
slib_require=pari
!exit
:proc
slib_out=
!distribute item $wims_read_parm into slib_data,slib_a, slib_b
slib_data=!declosing $slib_data
slib_n=!itemcnt $slib_data
slib_u=!exec pari gamma($slib_a)^(-1)
slib_out=!values $slib_b*exp(-$slib_b*x)*($slib_b*x)^($slib_a-1)*$slib_u for x in $slib_data
!for slib_r =1 to $slib_n
slib_s = !item $slib_r of $slib_data
!if $slib_s <= 0
slib_out=!replace internal item number $slib_r by 0 in $slib_out
!endif
!next slib_r
slib_out=!trim $slib_out