Subversion Repositories wimsdev

Rev

Rev 4349 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 4349 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=1\
 5 
slib_parms=1\
6 
,[data] or [data],[weight]
 6 
,[data] or [data],[weight]
7 
slib_author=Adeline, Grelot
 7 
slib_author=Adeline, Grelot
Line 27... Line 27...
27 
		slib_out=$[sqrt($slib_dev/$slib_cnt)],$[sqrt($slib_dev/($slib_cnt-1))]
 27 
		slib_out=$[sqrt($slib_dev/$slib_cnt)],$[sqrt($slib_dev/($slib_cnt-1))]
28 
		slib_out =!trim $slib_out
 28 
		slib_out =!trim $slib_out
29 
		!exit
 29 
		!exit
30 
	!else
 30 
	!else
31 
		slib_out=!exec pari slib_mean=([$slib_data]*[$slib_weight]~)/$slib_tw;\
 31 
		slib_out=!exec pari slib_mean=([$slib_data]*[$slib_weight]~)/$slib_tw;\
32 
slib_dev=sqrt(sum(x=1,$slib_cnt,[$slib_weight][x]*([$slib_data][x]-slib_mean)^2));\
 32 
    slib_dev=sqrt(sum(x=1,$slib_cnt,[$slib_weight][x]*([$slib_data][x]-slib_mean)^2));\
33 
print(slib_dev/sqrt($slib_tw)","slib_dev/sqrt($slib_tw-1))
 33 
      print(slib_dev/sqrt($slib_tw)","slib_dev/sqrt($slib_tw-1))
34 
slib_out =!trim $slib_out
 34 
    slib_out =!trim $slib_out
35 
	!endif weight
 35 
	!endif weight
36 
!endif
 36 
!endif
37 
 
 -