Subversion Repositories wimsdev

Rev

Rev 13586 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 13586 Rev 18193
Line 26... Line 26...
26
!endif
26
!endif
27
!if $slib_s>0 and $slib_s<1
27
!if $slib_s>0 and $slib_s<1
28
  slib_q=$slib_q, $[1-$slib_s]
28
  slib_q=$slib_q, $[1-$slib_s]
29
  !advance slib_t
29
  !advance slib_t
30
!endif
30
!endif
31
 
-
 
32
slib_mult=
31
slib_mult=
33
!for slib_r=1 to $slib_M
32
!for slib_r=1 to $slib_M
34
  slib_s=1
33
  slib_s=1
35
  slib_N=$slib_nn
34
  slib_N=$slib_nn
36
 
-
 
37
  slib_V=
35
  slib_V=
38
  !for slib_i=1 to  $[$slib_t-1]
36
  !for slib_i=1 to $[$slib_t-1]
39
    slib_j=!item $slib_i of $slib_q
37
    slib_j=!item $slib_i of $slib_q
40
    !readproc slib/stat/binomial  1,$slib_N, $[$slib_j/$slib_s]
38
    !readproc slib/stat/binomial 1,$slib_N, $[$slib_j/$slib_s]
41
 
-
 
42
    slib_V=!append item $slib_out to $slib_V
39
    slib_V=!append item $slib_out to $slib_V
43
    !distribute item $[$slib_s-$slib_j], $[$slib_N-$slib_out] into slib_s,slib_N
40
    !distribute item $[$slib_s-$slib_j], $[$slib_N-$slib_out] into slib_s,slib_N
44
 
-
 
45
  !next slib_i
41
  !next slib_i
46
  slib_V=!append item $slib_N to $slib_V
42
  slib_V=!append item $slib_N to $slib_V
47
 
-
 
48
  slib_mult=!append line $slib_V to $slib_mult
43
  slib_mult=!append line $slib_V to $slib_mult
49
!next slib_r
44
!next slib_r
50
 
45
 
51
slib_out=!nonempty line $slib_mult
46
slib_out=!nonempty line $slib_mult
52
slib_out=!trim $slib_out
47
slib_out=!trim $slib_out