Rev 7692 | 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 |
4161 | bpr | 4 | |
20 | reyssat | 5 | slib_parms=3\ |
6 | ,[data]\ |
||
7 | ,a\ |
||
8 | ,b |
||
9 | |||
4330 | bpr | 10 | slib_author=Sophie, Lemaire; Bernadette, Perrin-Riou |
20 | reyssat | 11 | slib_example=[0.2,0.4,0.7],6,8 |
12 | !exit |
||
13 | |||
14 | :proc |
||
15 | slib_out= |
||
16 | |||
4501 | bpr | 17 | !distribute item $wims_read_parm into slib_data,slib_a,slib_b |
20 | reyssat | 18 | |
4501 | bpr | 19 | slib_a=$[$slib_a] |
20 | slib_b=$[$slib_b] |
||
20 | reyssat | 21 | slib_data=!declosing $slib_data |
22 | |||
23 | |||
24 | slib_out= |
||
25 | !for slib_x in $slib_data |
||
13586 | bpr | 26 | !if $slib_x<=0 |
27 | slib_out=!append item 0 to $slib_out |
||
28 | !else |
||
29 | slib_out=!append item $[$slib_a*$slib_b*exp(-$slib_b*($slib_x)^($slib_a))*$slib_x^($slib_a-1)] to $slib_out |
||
30 | !endif |
||
20 | reyssat | 31 | !next slib_x |
32 | |||
33 | slib_out=!trim $slib_out |
||
34 | slib_data=$empty |