Numerical Integration of the ThomasFermi Equation from Zero to Infinity
Abstract
The ThomasFermi equation is expanded in a series at both x = 0 and x = ∞. Forward numerical integration from x = 0 determines an initial slope of 1.5880710226. Backward integration from x = ∞ depends on one parameter. By fitting the backward and forward integrations near x = 30, the parameter is determined and the numerical solution is obtained from zero to infinity with a high degree of accuracy.
 Publication:

Journal of Computational Physics
 Pub Date:
 August 1982
 DOI:
 10.1016/00219991(82)900833
 Bibcode:
 1982JCoPh..47..308K