Numerical analysis of a high order method for state-dependent delay integral equations

In this paper we study the piecewise collocation method for a class of functional integral equations with state‐dependent delays that is, where the delays depend on the solution. It is well known that these equations typically have discontinuity in the solution or its derivatives at the initial point of integration domain. This discontinuity propagated along the integration interval giving rise to subsequent points, called ‘singular points’, which can not be determined priori and the solution derivatives in these points are smoothed out along the interval. Most of the known numerical methods for this type of equations are generally very sensitive to the singular points and therefore must have a process that detects these points and insert them into the mesh to guarantee the required accuracy. Here, we present a numerical algorithm based on the piecewise collocation method and an approach for tracking the singular points relying on the recent strategy for implicit delay differential equations which has been proposed by Guglielmi and Hairer in 2008. The convergence analysis of the method is investigated and some numerical experiments are presented for clarifying the robustness of the method.