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