Rev 4330 | Details | Compare with Previous | Last modification | View Log | RSS feed
Rev | Author | Line No. | Line |
---|---|---|---|
20 | reyssat | 1 | !if $wims_read_parm!=slib_header |
13586 | bpr | 2 | !goto proc |
20 | reyssat | 3 | !endif |
4 | |||
4330 | bpr | 5 | slib_author=Sophie, Lemaire; Bernadette, Perrin-Riou |
4328 | bpr | 6 | slib_parms=3\ |
7 | ,[data]\ |
||
4330 | bpr | 8 | 1,mu\ |
4328 | bpr | 9 | 1,sigma |
4161 | bpr | 10 | |
4328 | bpr | 11 | slib_example=[2,4,7] ,3,4 |
12 | !exit |
||
20 | reyssat | 13 | |
14 | :proc |
||
15 | slib_out= |
||
16 | |||
17 | !distribute item $wims_read_parm into slib_data,slib_mean,slib_var |
||
18 | !default slib_mean=0 |
||
19 | !default slib_var=1 |
||
20 | slib_data=!declosing $slib_data |
||
21 | |||
22 | slib_out= |
||
23 | !for slib_x in $slib_data |
||
13586 | bpr | 24 | !if $slib_x <= 0 |
20 | reyssat | 25 | slib_out =!append item 0 to $slib_out |
26 | !else |
||
27 | slib_out = !append item $[(1+erf((log($slib_x)-$slib_mean)/(sqrt(2)*$slib_var)))*0.5] to $slib_out |
||
28 | !endif |
||
29 | !next slib_x |
||
30 | |||
31 | slib_out=!trim $slib_out |