Abstract: | A recursive numerical integration scheme based on the method of Clenshaw and Curtis is proposed for the efficient implementation of the variation–iteration procedure for the computation of approximate energies and wave functions for atomic systems. Extensive numerical tests are carried out to assess the accuracy and efficiency of the method and inaccuracies in some earlier calculations are pointed out. |