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 |
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 |
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 |