Re: [Bug-gsl] Error in gsl_cdf_fdist_Pinv

[Brian Gough]

At Fri, 10 Nov 2006 12:08:31 +1100,
Daniel Falster wrote:
> 
> The function double gsl_cdf_fdist_Pinv (double P, double nu1, double
> nu2) is in GSL ver1.8 but hasn't been included in earlier versions. It
> returns the inverse of the cumulative distribution function for the
> F-distribution. The arguments nu1 and nu2 are the degrees of freedom
> (1,..., Inf). When nu1 =1, and nu2 greater than ~ 250 the functions
> fails. I have verified this for both windows and unix versions.
> For example, 
>  for (int y = 1; y<500; y+=20)
>                      cout<<y<<"
> \t"<<gsl_cdf_fdist_Pinv(0.95,1.0,y)<<endl;
> fails for y > 240 on windows XP.

Hello,

Thanks for your email. Could you send the whole example program (with
#include's etc) and its complete output showing how it fails, so I can
reproduce the problem correctly.

-- 
best regards,

Brian Gough
(GSL Maintainer)

Network Theory Ltd,
Publishing the GSL Manual - http://www.network-theory.co.uk/gsl/manual/