Subversion Repositories wimsdev

!if $wims_read_parm!=slib_header
!goto proc
!endif
slib_title=Generation of lognormal random data
slib_parms=3\
,m\
,mean\
,sigma

slib_author=Bernadette PERRIN-RIOU
slib_out=random data of m reals following a lognormal law of parameters a and b . If m is a list [r,s], it returns an r by  s matrice. 
slib_comment=
slib_example= 20,0,1\
[6,5],0,2
!exit

:proc
slib_out=

!distribute  item $wims_read_parm into slib_data,slib_a,slib_b
slib_data=!declosing $slib_data
slib_cnt=!itemcnt slib_data
!if $slib_cnt=1
        slib_data=1,$slib_data
!endif  
slib_out=!exec octave lognormal_rnd($[exp($slib_a)],$[$slib_b^2], $slib_data)
slib_out=!words2items $slib_out

slib_out=!trim $slib_out