Incomplete Fermi Dirac ?

[giorgio mugnaini]

Hello,
I'm with University of Pisa and I've found that algorithm ACM 745 
computes incomplete Fermi-Dirac integrals.

An implementation is available in fortran 77 at:

http://www.netlib.org/toms/745

I've thought that this implementation could be used in order to build a 
port in C that can be included in GSL.

The program works fine with Lahey Fortran and Portland Group Fortran, 
but gives wrong results with G77 (bugs?).
Curiously, if it is translated in C with f2c, it works fine also with GCC!