Subversion Repositories wimsdev

Rev

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

Rev 7692 Rev 13586
Line 1... Line 1...
1 
!if $wims_read_parm!=slib_header
 1 
!if $wims_read_parm!=slib_header
2 
!goto proc
 2 
  !goto proc
3 
!endif
 3 
!endif
4 
 
 4 
 
5 
slib_parms=4\
 5 
slib_parms=4\
6 
,[data]\
 6 
,[data]\
7 
,N\
 7 
,N\
Line 32... Line 32...
32 
slib_PP=!values binomial($slib_K,x)*binomial($slib_N-$slib_K,$slib_n-x)*$slib_b for x = $slib_min to $slib_max
 32 
slib_PP=!values binomial($slib_K,x)*binomial($slib_N-$slib_K,$slib_n-x)*$slib_b for x = $slib_min to $slib_max
33 
!for slib_i in $slib_data
 33 
!for slib_i in $slib_data
34 
  !if $slib_i =1
 34 
  !if $slib_i =1
35 
   slib_out=!append item  $[$slib_max -$slib_min+1] to $slib_out
 35 
   slib_out=!append item  $[$slib_max -$slib_min+1] to $slib_out
36 
  !else
 36 
  !else
37 
	 slib_P=0
 37 
	  slib_P=0
38 
	 !for slib_k= 1 to $[$slib_max -$slib_min+1]
 38 
	  !for slib_k= 1 to $[$slib_max -$slib_min+1]
39 
		slib_a=!item $slib_k of $slib_PP
 39 
		  slib_a=!item $slib_k of $slib_PP
40 
		slib_P=$[$slib_P+$slib_a]
 40 
		  slib_P=$[$slib_P+$slib_a]
41 
		!if $slib_P>= $slib_i
 41 
		  !if $slib_P>= $slib_i
42 
			slib_out=!append item  $[$slib_k-1] to $slib_out
 42 
		  	slib_out=!append item  $[$slib_k-1] to $slib_out
43 
			!break
 43 
		  	!break
44 
		!endif
 44 
		  !endif
45 
	!next slib_k
 45 
	  !next slib_k
46 
 !endif
 46 
  !endif
47 
!next slib_i
 47 
!next slib_i
48 
 
 48 
 
49 
slib_out=!trim $slib_out
 49 
slib_out=!trim $slib_out