Rev 7692 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed
Rev 7692 | Rev 13586 | ||
---|---|---|---|
Line 1... | Line 1... | ||
1 | !if $wims_read_parm!=slib_header |
1 | !if $wims_read_parm!=slib_header |
2 | !goto proc |
2 | !goto proc |
3 | !endif |
3 | !endif |
4 | 4 | ||
5 | slib_parms=3\ |
5 | slib_parms=3\ |
6 | ,[data]\ |
6 | ,[data]\ |
7 | ,a\ |
7 | ,a\ |
Line 23... | Line 23... | ||
23 | slib_u=!exec pari gamma($slib_a)^(-1) |
23 | slib_u=!exec pari gamma($slib_a)^(-1) |
24 | 24 | ||
25 | slib_out=!values $slib_b*exp(-$slib_b*x)*($slib_b*x)^($slib_a-1)*$slib_u for x in $slib_data |
25 | slib_out=!values $slib_b*exp(-$slib_b*x)*($slib_b*x)^($slib_a-1)*$slib_u for x in $slib_data |
26 | 26 | ||
27 | !for slib_r =1 to $slib_n |
27 | !for slib_r =1 to $slib_n |
28 | slib_s = !item $slib_r of $slib_data |
28 | slib_s = !item $slib_r of $slib_data |
29 | !if $slib_s <= 0 |
29 | !if $slib_s <= 0 |
30 | slib_out=!replace internal item number $slib_r by 0 in $slib_out |
30 | slib_out=!replace internal item number $slib_r by 0 in $slib_out |
31 | !endif |
31 | !endif |
32 | !next slib_r |
32 | !next slib_r |
33 | 33 | ||
34 | slib_out=!trim $slib_out |
34 | slib_out=!trim $slib_out |