Subversion Repositories wimsdev

!if $wims_read_parm!=slib_header
  !goto proc
!endif
slib_parms= 2\
,[data]\
1,degree

slib_author=Sophie, Lemaire; Bernadette, Perrin-Riou

slib_example=[2,3,5],5
slib_require=pari
!exit

:proc
slib_out=

!distribute  item $wims_read_parm into slib_data,slib_degree
!default slib_degree= 1
slib_data= !declosing $slib_data
slib_n= !itemcnt $slib_data
slib_u= !exec pari gamma($slib_degree*0.5)^(-1)
!for slib_i in $slib_data
        !if $slib_i > 0
                slib_val=$[0.5*exp(-0.5*$slib_i)*(0.5*$slib_i)^( $slib_degree/2-1)*$slib_u]
        !else
                slib_val=0
        !endif
        slib_out=!append item $slib_val to $slib_out
!next i
slib_out=!nonempty items $slib_out
slib_out= !trim $slib_out