gsl-chisquare cumulative density function

Gilles Frising frising@mail.desy.de
Wed Dec 19 13:20:00 GMT 2001


Hello there,

I wanted to use the cumulative chisquare pdf of gsl:
`gsl_ran_chisq_pdf'.
I checked this function in maple,cernlib and gsl, and i found different
result for gsl.



MAPLE:

statevalf[cdf,chisquare[6]](3.);
      .1911531695 
       ==========  



CERNLIB:

function: http://wwwinfo.cern.ch/asdoc/shortwrupsdir/g100/top.html

  extern "C" float prob_(float&,int&);
  int n=6;
  float x=3.;
  cout<<"prob: "<<1.-prob_(x,n)<<endl;

prob: 0.191153 
      ========  




GSL:

gsl_ran_chisq_pdf(3.,6.):  0.125511
                           =========


Is this a bug or have i done sth. wrong?

ciao,

Gilles



More information about the Gsl-discuss mailing list