Two-point boundary value problems associated to functional differential equations of even order solved by iterated splines

A new iterative numerical method to solve two-point boundary value problems associated to functional differential equations of even order is proposed. The method uses a cubic spline interpolation procedure activated at each iterative step. The convergence of the method is proved and it is tested on some numerical experiments. The notion of numerical stability with respect to the choice of the first iteration is introduced proving that the proposed method is numerically stable in this sense.