Subversion Repositories wimsdev

!if $wims_read_parm!=slib_header

 !goto proc

!endif

slib_parms=3\

 1,m\

  ,n\

 ,[p1,p2, .. pj]

slib_author=Sophie, Lemaire; Bernadette, PERRIN-RIOU

slib_example=3,6,[1/3,1/3]

!exit

:proc

!distribute item $wims_read_parm into slib_M, slib_nn

slib_q=!item 3 to -1 of $wims_read_parm

!default slib_M=1

slib_q=!declosing $slib_q

slib_t=!itemcnt $slib_q

slib_s=!sum x for x in $slib_q

!if $slib_s>1 

!exit

!endif

!if $slib_s>0 and $slib_s<1

slib_q=$slib_q, $[1-$slib_s]

!advance slib_t

!endif

slib_mult=

!for slib_r=1 to $slib_M

slib_s=1

slib_N=$slib_nn

slib_V=

!for slib_i=1 to  $[$slib_t-1]

slib_j=!item $slib_i of $slib_q

!readproc slib/stat/binomial  1,$slib_N, $[$slib_j/$slib_s]

slib_V=!append item $slib_out to $slib_V

!distribute item $[$slib_s-$slib_j], $[$slib_N-$slib_out] into slib_s,slib_N

!next slib_i

slib_V=!append item $slib_N to $slib_V

slib_mult=!append line $slib_V to $slib_mult

!next slib_r

slib_out=!nonempty line $slib_mult

slib_out=!trim $slib_out