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

slib_parms=2\
  1,m\
  1,lambda
slib_author=Sophie, Lemaire; Bernadette, Perrin-Riou

slib_example= 5,2.5
slib_require=octave pari
!exit

:proc
slib_out=
!distribute item $wims_read_parm into slib_data, slib_lambda

!default slib_m=1
!default slib_lambda=1

!!slib_ex=$[exp(-$slib_lambda)]
!!slib_q=$slib_ex
!!
!!
!!slib_U=!random 0,1 repeat $slib_m
!!
!!!for slib_i=1 to $slib_m
!! slib_V=!item $slib_i of $slib_U
!! slib_q=$slib_ex
!! slib_p=$slib_ex
!! slib_j=0
!!
!!!while $slib_V >= $slib_q
!!  !advance slib_j
!!  slib_p=$[$slib_p*$slib_a/$slib_j]
!!  slib_q=$[$slib_q+$slib_p]
!!!endwhile
!!
!!slib_out=!append item $slib_j to $slib_out
!!!next slib_i


slib_data=!declosing $slib_data
slib_cnt=!itemcnt $slib_data
!if $slib_cnt=1
  slib_data=1,$slib_data
!endif

slib_test=!exec octave pkg list statistics
!if Installation iswordof $slib_test
  slib_out=!exec octave poissrnd($slib_lambda,$slib_data)
!else
  !readproc slib/stat/commun stat
  !for slib_x=1 to $(slib_data[1])
    !reset slib_line
    !for slib_xx=1 to $(slib_data[2])
      slib_res=!exec pari statpoissrandom($slib_lambda)
      slib_line=!append item $slib_res to $slib_line
    !next
    slib_out=!append line $slib_line to $slib_out
  !next
!endif
slib_out=!trim $slib_out
slib_out=!replace internal $\
$ by ; in $slib_out
slib_out=!words2items $slib_out
