Numerical detection and study of singularities in solutions of differential equations

Abstract

New simple and robust methods are proposed for detecting singularities, such as poles, logarithmic poles, and mixed singularities, in systems of ordinary differential equations. The methods produce characteristics of these singularities with an a posteriori asymptotically precise error estimate. They are applicable in the case of an arbitrary parametrization of integral curves, including one in terms of the arc length, which is optimal for stiff and ill-conditioned problems. Following this approach, blowup solutions can be detected for a broad class of important nonlinear partial differential equations, since they are reducible by the method of lines to systems of ordinary differential equations of huge orders. The simplicity and reliability of the approach are superior to those of previously known methods.