Reciprocal function method for Cauchy problems with first-order poles

For the numerical solution of the Cauchy problem with multiple poles, we propose a reciprocal function method. In the case of first-order poles, it makes it possible to continue the solution through the poles and to determine the solution and the pole positions with good accuracy. The method allows one to employ conventional explicit and implicit schemes, for example, explicit Runge–Kutta schemes. A test problem with multiple poles is computed as an example. The proposed method is useful for construction of software for direct computation of special functions.