About a numerical method of successive interpolations for two-point boundary value problems with deviating argument (Q544081)

The authors provide a numerical method for the boundary value problem \[ x''(t) = f(t, x(t), x(\phi(t))),\;x(0) = b,\;x(a) = c, \] where \(\phi\) is a Lipschitz function such that \(0 \leq \phi(t) \leq a\) for \(t \in [0,a]\). The function \(f\) is also assumed to satisfy a Lipschitz condition in each argument. The paper begins by stating sufficient conditions for the problem to have a unique solution. Then the numerical method is developed. This method is based on a uniform partition of the basic interval \([0,a]\) into \(n\) subintervals of length \(h = a/n\). The differential equation is rewritten in the form of an equivalent integral equation. Then a trapezoidal quadrature formula is used to approximate the integrals over the subintervals mentioned above. This approximation is used to compute the coefficients of a piecewise defined interpolation scheme of lacunary Birkhoff type. A brief convergence and stability analysis is provided together with a large number of numerical examples.