Subversion Repositories wimsdev

Rev

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