On the inversion of sparse matrices
A. L. Dulmage, N. S. Mendelsohn · Mathematics of Computation · 1962
choice of E. The usual procedure is to divide the range of integration into two parts, integrate outwards for a solution satisfying one boundary condition, integrate inwards for a solution satisfying the other boundary condition, match the solutions at an intermediate point and adjust E so that the derivatives also agree [1], [2]. The inward integration may be avoided with the procedure described earlier. A convenient way of dividing the range is according to the sign of f(r). For some r, f(r) < 0 so that condition (ii) is not satisfied: the procedure described here is not always numerically stable when f(r) < 0 [3]; in fact, for some values of i, I d, I < 1. Of a series of standard methods, the Numerov method,