An iterative approach to approximate global solutions of high order initial value problems (IVPs):
with is proposed. The approach extends previous results of the same authors for first order non linear IVPs [Nonlinear Anal., Theory Methods Appl. 63, No. 1 (A), 97–105 (2005; Zbl 1097.34005)].
Under continuity and uniform Lipschitz conditions on the nonlinear function , the IVP (1) is equivalent to an integral equation with some integral operator and the problem reduces to the computation of a fixed point of in a Banach space. For the approximation of such a fixed point a suitable Faber-Schauder basis of piecewise linear functions associated to a sequence of nodes with dense in is proposed such that with appropriately chosen tends to the exact solution. Several convergence results for an iterative to the solution as well as error estimates are proposed and finally the technique is applied to an IVP second order test problem to show the convergence of the method depending on the number of the nodes and the number of iterations.