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