Lobatto-Obrechkoff Formulae for 2nd Order Two-Point Boundary Value Problems

A substantial increase in efficiency may be obtained by numerical integration methods which take advantage of the special second order forms y" = f(x, y) or y" = f(x, y, y') in systems of second order two-point boundary value problems, while retaining the MIRK structure. In particular, for these special second order equations, we derive high order methods which require considerably fewer function evaluations than are required by methods intended for general first order systems. Methods based on a Lobatto-MIRK formula for finding y'n and an Obrechkoff type formula for finding y n using the same values of y" are derived and these methods are of sixth and of eighth order accuracy.