Subversion Repositories wimsdev

Rev

Rev 4501 | Rev 7692 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 4501 Rev 7277
Line 7... Line 7...
7 
5,n\
 7 
5,n\
8 
0.5,p
 8 
0.5,p
9 
 
 9 
 
10 
slib_author=Sophie, Lemaire; Bernadette, Perrin-Riou
 10 
slib_author=Sophie, Lemaire; Bernadette, Perrin-Riou
11 
 
 11 
 
12 
slib_example=[0.2,0.4,0.6],5,0.5
 12 
slib_example=0.5,11,0.5\
13 
!exit
 13 
0.5,4,0.5\
- 
 
 14 
[0.5,0.3,0.6],101,0.5\
14 
 
 15 
0.3,800,0.61\
15 
:proc
 16 
0.3,900,0.34\
- 
 
 17 
[0.2,0.4,0.6],5,0.5\
16 
 
 18 
 
- 
 
 19 
!exit
- 
 
 20 
 
- 
 
 21 
:proc
- 
 
 22 
!reset slib_out
17 
!distribute  item $wims_read_parm into slib_data,slib_n,slib_p
 23 
!distribute  item $wims_read_parm into slib_data,slib_n,slib_p
18 
!default slib_n= 5
 24 
!default slib_n= 5
19 
!default slib_p=0.5
 25 
!default slib_p=0.5
20 
!if $slib_p>=1 or $slib_p<=0
26 
!if $slib_p>=1 or $slib_p<=0
21 
	slib_out=0
 27 
	slib_out=0
22 
	!exit
 28 
	!exit
23 
!endif
 29 
!endif
24 
slib_out=
 30 
slib_out=
25 
slib_data=!declosing $slib_data
 31 
slib_data=!declosing $slib_data
26 
slib_p=$[$slib_p]
 32 
slib_p=$[$slib_p]
- 
 
 33 
 
- 
 
 34 
!if $slib_p=0.5 and 0.5 isitemof $slib_data
- 
 
 35 
  !for slib_q in $slib_data
- 
 
 36 
   !if $slib_q=0.5
- 
 
 37 
     !if $[$slib_n%2]=0
- 
 
 38 
       slib_out=!append item $[$slib_n/2] to $slib_out
- 
 
 39 
     !else
- 
 
 40 
       slib_out=!append item $[($slib_n-1)/2] to $slib_out
- 
 
 41 
     !endif
- 
 
 42 
   !else
- 
 
 43 
     slib_res=!exec octave binoinv($slib_q, $slib_n, $slib_p)
- 
 
 44 
     slib_out=!append item $slib_res to $slib_out
- 
 
 45 
   !endif
- 
 
 46 
  !next
- 
 
 47 
  !goto end
- 
 
 48 
!endif
- 
 
 49 
 
- 
 
 50 
slib_out=!exec octave binoinv([$slib_data], $slib_n, $slib_p)
- 
 
 51 
 
- 
 
 52 
:end
- 
 
 53 
 
- 
 
 54 
slib_out=!words2items $slib_out
- 
 
 55 
slib_out=!trim $slib_out
- 
 
 56 
 
- 
 
 57 
!exit
- 
 
 58 
old
- 
 
 59 
 
27 
slib_p1=$[1-$slib_p]
 60 
slib_p1=$[1-$slib_p]
28 
slib_r=$[$slib_p/$slib_p1]
 61 
slib_r=$[$slib_p/$slib_p1]
29 
slib_a=$[$slib_p1^($slib_n)]
 62 
slib_a=$slib_p1^($slib_n)
30 
!!slib_bino=!values binomial($slib_n,x)*$slib_p^x*$slib_p1^($slib_n-x) for x =0 to $slib_n
 - 
 
31 
slib_cnt= !itemcnt $slib_data
 63 
slib_cnt= !itemcnt $slib_data
32 
slib_pp=$[$slib_p/$slib_p1]
 64 
slib_pp=$[$slib_p/$slib_p1]
- 
 
 65 
 
33 
!for slib_i =1 to $slib_cnt
 66 
!for slib_i =1 to $slib_cnt
34 
	slib_di=!item $slib_i of $slib_data
 67 
	slib_di=!item $slib_i of $slib_data
35 
	slib_P=0
 68 
	slib_P=0
36 
	slib_a=$[$slib_p1^($slib_n)]
 69 
	slib_a=$[$slib_p1^($slib_n)]
37 
	!for slib_k=0 to $slib_n
 70 
	!for slib_k=0 to $slib_n
38 
		!!slib_a=!item $[$slib_k+1] of $slib_bino
 71 
		!!slib_a=!item $[$slib_k+1] of $slib_bino
39 
		slib_P=$[$slib_P+$slib_a]
 72 
		slib_P=$[$slib_P+$slib_a]
- 
 
 73 
		slib_out=!append item $slib_P to $slib_out
40 
		!if $slib_P>= $slib_di
 74 
		!ifval $slib_P>= $slib_di
41 
			slib_out=!append item  $slib_k to $slib_out
 75 
			slib_out=!append item  $slib_k to $slib_out
42 
			!break
 76 
			!break
43 
		!endif
 77 
		!endif
44 
		slib_a=$[$slib_a*$slib_pp*($slib_n-$slib_k)/($slib_k+1)]
 78 
		slib_a=$[$slib_a*$slib_pp*($slib_n-$slib_k)/($slib_k+1)]
45 
	!next slib_k
 79 
	!next slib_k
46 
!next slib_i
 80 
!next slib_i
- 
 
 81 
!default slib_out=$slib_k
47 
slib_out=!trim $slib_out
 82 
slib_out=!trim $slib_out