[Bug-gsl] possible bug in F distribution computation function gsl_ran_fdist_pdf

[Gunjan Verma]

Dear GSL Community,
    Ran into a mysterious problem.

When I run gsl_ran_fdist_pdf(.012, 8, 249) I get infinity.
But when I run

gsl_ran_fdist_pdf(.012, 8, 248) I get .000074, which is correct (checked
against MATLAB and R)

one d.f. change makes this difference - but I don't see why it should from
the formula for F pdf.

i am running a Dell Latitude E6500 machine with ubuntu 9 linux , 64 bit and
gcc 4.4.1   . I built the gsl library from source for my system.

Please help!

Gunjan