Rev 4330 | 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 | 0.5,p\ |
||
8 | 1000,M |
||
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],0.5 |
13 | !exit |
||
14 | |||
15 | :proc |
||
16 | |||
17 | !distribute item $wims_read_parm into slib_data,slib_p,slib_N |
||
18 | !default slib_p=0.5 |
||
19 | !default slib_N=1000 |
||
20 | !if $slib_p<=0 or $slib_p>=1 |
||
21 | slib_out=0 |
||
22 | !exit |
||
23 | !endif |
||
24 | slib_out= |
||
25 | slib_data=!declosing $slib_data |
||
26 | |||
27 | slib_cnt= !itemcnt $slib_data |
||
4501 | bpr | 28 | slib_p=$[$slib_p] |
20 | reyssat | 29 | slib_pp=$[1-$slib_p] |
30 | !for slib_i =1 to $slib_cnt |
||
31 | slib_di=!item $slib_i of $slib_data |
||
32 | slib_P=0 |
||
33 | slib_a=$slib_p |
||
34 | !for slib_k=1 to $slib_N |
||
35 | slib_P=$[$slib_P+$slib_a] |
||
36 | !if $slib_P>= $slib_di |
||
37 | slib_out=!append item $slib_k to $slib_out |
||
38 | !break |
||
39 | !endif |
||
40 | slib_a=$[$slib_a*$slib_pp] |
||
41 | !next slib_k |
||
42 | !next slib_i |
||
43 | slib_out=!trim $slib_out |