Rev 4328 | 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 |
||
4 | |||
4330 | bpr | 5 | slib_author=Sophie, Lemaire; Bernadette, Perrin-Riou |
4161 | bpr | 6 | |
4328 | bpr | 7 | slib_parms=3\ |
8 | ,[data]\ |
||
4330 | bpr | 9 | 1,mu\ |
4328 | bpr | 10 | 1,sigma |
11 | |||
12 | slib_example=[2,4,7] ,3,4\ |
||
13 | 2,2,sqrt(2)\ |
||
14 | 2,2,2\ |
||
15 | 2,1,1 |
||
20 | reyssat | 16 | !exit |
17 | |||
18 | :proc |
||
19 | slib_out= |
||
20 | |||
4328 | bpr | 21 | !distribute item $wims_read_parm into slib_data,slib_mean,slib_var |
20 | reyssat | 22 | !default slib_mean=1 |
23 | !default slib_var=1 |
||
4328 | bpr | 24 | !! en fait c'est un écart type |
20 | reyssat | 25 | slib_data=!declosing $slib_data |
26 | |||
27 | slib_out= |
||
28 | !for slib_i in $slib_data |
||
29 | !if $slib_i > 0 |
||
30 | slib_out = !append item $[exp(-(log($slib_i)-$slib_mean)^2/(2*$slib_var^2))/($slib_i*$slib_var*(2*pi)^0.5)] to $slib_out |
||
31 | !else |
||
32 | slib_out = !append item 0 to $slib_out |
||
33 | !endif |
||
34 | !next slib_i |
||
35 | |||
36 | |||
37 | slib_out=!trim $slib_out |