Subversion Repositories wimsdev

!if $wims_read_parm!=slib_header

!goto proc

!endif

slib_title=Cumulative distribution function of chi2 

slib_parms=2\

,[data]\

1,degree

slib_author=Bernadette PERRIN-RIOU

slib_out=values of cumulative distribution function of chi2  on data

slib_comment=data can be a list of reals in brackets

slib_example=  [2,4,7] ,5

slib_require=pari

!exit

:proc

slib_out=

!distribute  item $wims_read_parm into slib_data,slib_degree

!default slib_degree=1

slib_data=!declosing $slib_data

slib_n=!itemcnt $slib_data

slib_out=!exec pari {slib_u=gamma($slib_degree/2);\

print(vector($slib_n,x,if([$slib_data][x]>0,incgamc($slib_degree/2,[$slib_data][x]/2)/slib_u,0)))}

slib_out=!trim $slib_out