A multiquadric interpolation method for solving initial value problems (Q1387813)

The authors propose a new interpolation method for the numerical solution of the initial value problem for an \(n\)th order linear ordinary differential equation. The method is based on the application of the multiquadric scheme for global interpolation of the solution and then on the collocation to obtain the equations for the unknown coefficients. The resulting system of equations is solved by using Gaussian elimination with partial pivoting. Two numerical examples which, in the authors' opinion, show that the method under consideration offers a higher degree of accuracy than some of the known methods (in particular, such as the fourth- order Runge-Kutta method) are given.