The Half-order Fermi-Dirac integral is needed extensively in solid state physics applications. Here, a trapezoidal integration with the correction for the poles of the integrand is used to evaluate the integrand. It is a beautiful application of the classical Cauchy theorem.
This is extremely well written code. I've been using the Purdue nanohub numerical integration code and your version runs circles around it. The accuracy is as good as the numerical integration. Well done, and thank you!
15 Jan 2008
Thank you so much! It works execellent, and I really need it