A Taylor series method for the numerical solution of two-point boundary value problems

Summary For the numerical solution of two-point boundary value problems a shooting algorithm based on a Taylor series method is developed. Series coefficients are generated automatically by recurrence formulas. The performance of the algorithm is demonstrated by solving six problems arising in nonlinear shell theory, chemistry and superconductivity.