Multiderivative methods for nonlinear beam problems

Second, fourth and sixth order methods are developed and analysed for the numerical solution of nonlinear fourth order boundary value problems with boundary conditions given as two function values together with two second order derivatives or two first order derivatives. The methods arise from a five‐point recurrence relation involving exponential terms, these being replaced by Padé approximants. The methods are tested on two problems from the literature.