Fermi-Dirac integral

From: Jacek M. Holeczek <holeczek_at_us.edu.pl>
Date: Thu, 20 Jan 2005 14:18:46 +0100 (MET)

Sorry for this unrelated to ROOT question, but I hope someone here could know the answer.

I need to calculate some Fermi-Dirac integrals (0 to infinity) for functions in form :
> x^n
> f(x,n,u) = ---------------
> (e^(x - u) + 1)

where n = 2,3,4 (natural) and u = something between 0.0 and 5.0 (real)

Analytical formulas, that I know, use either the Lerch transcendent, or the Jonquiere's (known also as the polylogarithm) function.

I would like to use functions coming from CERNLIB, but none of the above mentioned is implemented. I think there must be another magic function in CERNLIB which can be used to calculate this interval, I just miss the magic formula to do this ...

Thanks in advance,
Best regards,
Jacek. Received on Thu Jan 20 2005 - 14:18:54 MET

This archive was generated by hypermail 2.2.0 : Tue Jan 02 2007 - 14:45:04 MET