Subversion Repositories wimsdev

Rev

Rev 4328 | Rev 7692 | Go to most recent revision | Details | Compare with Previous | Last modification | View Log | RSS feed

Rev Author Line No. Line
20 reyssat 1
!if $wims_read_parm!=slib_header
2
 !goto proc
3
!endif
4161 bpr 4
 
20 reyssat 5
slib_parms=3\
6
 1,m\
7
  ,n\
8
 ,[p1,p2, .. pj]
4330 bpr 9
slib_author=Sophie, Lemaire; Bernadette, PERRIN-RIOU
20 reyssat 10
 
11
slib_example=3,6,[1/3,1/3]
3307 bpr 12
!exit
20 reyssat 13
 
14
:proc
15
 
16
!distribute item $wims_read_parm into slib_M, slib_nn
17
slib_q=!item 3 to -1 of $wims_read_parm
18
!default slib_M=1
19
 
20
slib_q=!declosing $slib_q
21
slib_t=!itemcnt $slib_q
22
 
23
slib_s=!sum x for x in $slib_q
24
!if $slib_s>1 
25
!exit
26
!endif
27
!if $slib_s>0 and $slib_s<1
28
slib_q=$slib_q, $[1-$slib_s]
29
!advance slib_t
30
!endif
31
 
32
slib_mult=
33
!for slib_r=1 to $slib_M
34
slib_s=1
35
slib_N=$slib_nn
36
 
37
slib_V=
38
!for slib_i=1 to  $[$slib_t-1]
39
slib_j=!item $slib_i of $slib_q
40
!readproc slib/stat/binomial  1,$slib_N, $[$slib_j/$slib_s]
41
 
42
slib_V=!append item $slib_out to $slib_V
4330 bpr 43
!distribute item $[$slib_s-$slib_j], $[$slib_N-$slib_out] into slib_s,slib_N
20 reyssat 44
 
45
!next slib_i
46
slib_V=!append item $slib_N to $slib_V
47
 
48
slib_mult=!append line $slib_V to $slib_mult
49
!next slib_r
50
 
4328 bpr 51
slib_out=!nonempty line $slib_mult
20 reyssat 52
slib_out=!trim $slib_out