Subversion Repositories wimsdev

Rev

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

Rev 13586 Rev 18193
Line 19... Line 19...
19
 
19
 
20
!exit
20
!exit
21
 
21
 
22
:proc
22
:proc
23
!reset slib_out
23
!reset slib_out
24
!distribute  item $wims_read_parm into slib_data,slib_n,slib_p
24
!distribute item $wims_read_parm into slib_data,slib_n,slib_p
25
!default slib_n= 5
25
!default slib_n= 5
26
!default slib_p=0.5
26
!default slib_p=0.5
27
!if $slib_p>=1 or $slib_p<=0
27
!if $slib_p>=1 or $slib_p<=0
28
	slib_out=0
28
	slib_out=0
29
	!exit
29
	!exit
Line 71... Line 71...
71
	!for slib_k=0 to $slib_n
71
	!for slib_k=0 to $slib_n
72
		!!slib_a=!item $[$slib_k+1] of $slib_bino
72
		!!slib_a=!item $[$slib_k+1] of $slib_bino
73
		slib_P=$[$slib_P+$slib_a]
73
		slib_P=$[$slib_P+$slib_a]
74
		slib_out=!append item $slib_P to $slib_out
74
		slib_out=!append item $slib_P to $slib_out
75
		!ifval $slib_P>= $slib_di
75
		!ifval $slib_P>= $slib_di
76
			slib_out=!append item  $slib_k to $slib_out
76
			slib_out=!append item $slib_k to $slib_out
77
			!break
77
			!break
78
		!endif
78
		!endif
79
		slib_a=$[$slib_a*$slib_pp*($slib_n-$slib_k)/($slib_k+1)]
79
		slib_a=$[$slib_a*$slib_pp*($slib_n-$slib_k)/($slib_k+1)]
80
	!next slib_k
80
	!next slib_k
81
!next slib_i
81
!next slib_i