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=[0.2,0.4,0.7],6,8
!exit

:proc
slib_out=

!distribute item $wims_read_parm into slib_data,slib_a,slib_b

slib_a=$[$slib_a]
slib_b=$[$slib_b]
slib_data=!declosing $slib_data


slib_out=
!for slib_x in $slib_data
  !if $slib_x<=0
    slib_out=!append item 0 to $slib_out
  !else
    slib_out=!append item $[$slib_a*$slib_b*exp(-$slib_b*($slib_x)^($slib_a))*$slib_x^($slib_a-1)] to $slib_out
  !endif
!next slib_x

slib_out=!trim $slib_out
slib_data=$empty