Rev 4330 | Blame | Compare with Previous | Last modification | View Log | RSS feed
!if $wims_read_parm!=slib_header
!goto proc
!endif
slib_author=Adeline, Grelot; Bernadette, PERRIN-RIOU
slib_example=[2,-7,6]\
[2,1;-7,1;6,4]\
[2,-7,6],[1,1,4]\
[2,-7,6;1,1,4]
slib_require=pari
!exit
:proc
!readproc slib/stat/dataproc $wims_read_parm
!if $slib_cnt=0
slib_out=0
!else
!if $slib_weight =$empty
slib_out=!sum x^2 for x in $slib_data
slib_out=$[sqrt($slib_out/$slib_cnt)]
!else
slib_cnt=!itemcnt $slib_data
slib_out=!exec pari print(sqrt(sum(x=1,$slib_cnt,[$slib_weight][x]*([$slib_data][x])^2)/$slib_tw))
!endif weight
!endif
slib_out=!trim $slib_out