Subversion Repositories wimsdev

Rev

Rev 4352 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 4352 Rev 7692
Line 10... Line 10...
10 
slib_example=[2,3,5],5
 10 
slib_example=[2,3,5],5
11 
slib_require=pari
 11 
slib_require=pari
12 
!exit
 12 
!exit
13 
 
 13 
 
14 
:proc
 14 
:proc
15 
slib_out=
15 
slib_out=
16 
 
 16 
 
17 
!distribute  item $wims_read_parm into slib_data,slib_degree
 17 
!distribute  item $wims_read_parm into slib_data,slib_degree
18 
!default slib_degree= 1
 18 
!default slib_degree= 1
19 
slib_data= !declosing $slib_data
 19 
slib_data= !declosing $slib_data
20 
slib_n= !itemcnt $slib_data
 20 
slib_n= !itemcnt $slib_data
21 
slib_u= !exec pari gamma($slib_degree*0.5)^(-1)
 21 
slib_u= !exec pari gamma($slib_degree*0.5)^(-1)
22 
!for slib_i in $slib_data
 22 
!for slib_i in $slib_data
23 
	!if $slib_i > 0
 23 
	!if $slib_i > 0
24 
	 	slib_val=$[0.5*exp(-0.5*$slib_i)*(0.5*$slib_i)^( $slib_degree/2-1)*$slib_u]
 24 
	 	slib_val=$[0.5*exp(-0.5*$slib_i)*(0.5*$slib_i)^( $slib_degree/2-1)*$slib_u]
25 
	!else
25 
	!else
26 
		slib_val=0
 26 
		slib_val=0
27 
	!endif
 27 
	!endif
28 
	slib_out=!append item $slib_val to $slib_out
 28 
	slib_out=!append item $slib_val to $slib_out
29 
!next i
 29 
!next i
30 
slib_out=!nonempty items $slib_out
 30 
slib_out=!nonempty items $slib_out