Subversion Repositories wimsdev

!if $wims_read_parm!=slib_header

!goto proc

!endif

slib_title=Probability density function of Weibull law

slib_parms=3\

,[data]\

,a\

,b

slib_author=Bernadette PERRIN-RIOU

slib_out=values of probability density function of Weibull law on data

slib_comment=data can be a list of reals in brackets; density is a*b*exp(-b*x^a)*x^(a-1) for x>=0

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_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