Subversion Repositories wimsdev

!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