Methods for numerically approximating solutions of Thomas-Fermi equati
ons are considered, Iterative algorithms are described which bound the
solutions. New discretisations for the algorithms provide O(h2) accur
acy, which is then improved by extrapolation.