Research Articles
Published 2021-08-28
Keywords
- Fermi-Dirac Integrals,
- Partial Derivatives,
- QUADPACK Adaptive Quadrature,
- Optimised Break Points
Copyright (c) 2019

This work is licensed under a Creative Commons Attribution-ShareAlike 4.0 International License.
Abstract
In this paper, using the Python SciPy module “quad”, a fast auto-adaptive quadrature solver based on the pre-compiled QUADPACK Fortran package, computational research is undertaken to accurately integrate the generalised Fermi-Dirac function and all its partial derivatives up to the third order. The numerical results obtained with quad method when combined with optimised break points achieve an excellent accuracy comparable to that obtained by other publications using fixed-order quadratures.
References
[1] J. McDougall and E.C. Stoner, Phil. Trans. Roy. Soc. London, A, Vol 231, 67, 1939.
[2] R.B. Larson, and P.R. Demarque, An application of Henyey's approach to the integration of the equations of stellar structure, The Astrophysical Journal, 140, 524, 1964.
[3] D.D. Clayton, Principles of stellar evolution and nucleosynthesis, McGraw-Hill, 1968
[4] J.P. Cox and R.T. Giuli, 1968, Principles of stellar structure, Gordon & Breach, 1968.
[5] L.D. Cloutman, Numerical evaluation of FDIs, Ap. J. Supp., 7, 677-699, 1989.
[6] F.X.Timmes, http:// cococubed.asu.edu/ code_pages/ fermi_dirac.shtml
[7] F.X.Timmes, http://cococubed.asu.edu/code_pages/eos.shtml[
[8] Z. Gong etal., “Generalized FD Functions and Derivatives: Properties and Evaluation”, arXiv:astro-ph/0102329v1.20Feb.2001.
[9] http://www.cpc.cs.qub.ac.uk/summaries/ADNX_v1_0.html
[10] http://nines.cs.kuleuven.be/software/QUADPACK/
[11] https:// github.com /scipy/scipy/ tree/master/scipy/ integrate/ quadpack
[12] J.D. Cook, https://www.johndcook.com/blog/2012/ 03/20/scipy-integration/
[2] R.B. Larson, and P.R. Demarque, An application of Henyey's approach to the integration of the equations of stellar structure, The Astrophysical Journal, 140, 524, 1964.
[3] D.D. Clayton, Principles of stellar evolution and nucleosynthesis, McGraw-Hill, 1968
[4] J.P. Cox and R.T. Giuli, 1968, Principles of stellar structure, Gordon & Breach, 1968.
[5] L.D. Cloutman, Numerical evaluation of FDIs, Ap. J. Supp., 7, 677-699, 1989.
[6] F.X.Timmes, http:// cococubed.asu.edu/ code_pages/ fermi_dirac.shtml
[7] F.X.Timmes, http://cococubed.asu.edu/code_pages/eos.shtml[
[8] Z. Gong etal., “Generalized FD Functions and Derivatives: Properties and Evaluation”, arXiv:astro-ph/0102329v1.20Feb.2001.
[9] http://www.cpc.cs.qub.ac.uk/summaries/ADNX_v1_0.html
[10] http://nines.cs.kuleuven.be/software/QUADPACK/
[11] https:// github.com /scipy/scipy/ tree/master/scipy/ integrate/ quadpack
[12] J.D. Cook, https://www.johndcook.com/blog/2012/ 03/20/scipy-integration/