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

slib_parms=3\
,m\
,n\
,p

slib_author=Sophie, Lemaire; Bernadette, Perrin-Riou

slib_example=5,4,0.3\
[2,3],4,0.3
slib_require=octave pari
!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_test=!exec octave pkg list statistics
!if  Installation iswordof $slib_test
  slib_out=!exec octave nbinrnd($slib_a,$slib_b,$slib_data) + $slib_a
!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 statnegbinrandom($slib_a,$slib_b)+$slib_a
      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
