Subversion Repositories wimsdev

Rev

Rev 4161 | 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
,[data]\
7
1,p\
8
1000,N
9
 
4330 bpr 10
slib_author=Sophie, Lemaire; Bernadette, Perrin-Riou
4161 bpr 11
 
20 reyssat 12
slib_example=[0.2,0.4,0.7],6.1
13
!exit
14
:proc
15
 
16
!distribute  item $wims_read_parm into slib_data,slib_p,slib_N
17
 
18
!default slib_N=1000
19
!if  $slib_p<=0
20
slib_out=0
21
!exit
22
!endif
23
slib_out=
24
slib_data=!declosing $slib_data
25
 
26
slib_cnt= !itemcnt $slib_data
27
!for slib_i =1 to $slib_cnt
28
slib_di=!item $slib_i of $slib_data
29
slib_P=0
30
slib_a=$[exp(-$slib_p)]
31
!for slib_k=0 to $slib_N
32
slib_P=$[$slib_P+$slib_a]
33
!if $slib_P>= $slib_di
34
slib_out=!append item  $slib_k to $slib_out
35
!break
36
!endif
37
slib_a=$[$slib_p*$slib_a/($slib_k+1)]
38
!next slib_k
39
!next slib_i
40
 
41
slib_out=!trim $slib_out