Subversion Repositories wimsdev

Rev

Rev 3265 | Rev 4352 | Go to most recent revision | Blame | Compare with Previous | Last modification | View Log | RSS feed

!if $wims_read_parm!=slib_header
!goto proc
!endif

slib_author= 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