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