Numerical method for solving volterra integral equations with a convolution kernel

This paper presents a numerical method for solving the Volterra integral equation with a convolution kernel. The integral equation was first converted to an algebraic equation using the Laplace transform, after which its numerical inversion was determined by power series. The Padé approximants were effectively used to improve the convergence rate and accuracy of the computed series. The method is described and illustrated with numerical examples. The results revealed that the method is accurate and easy to implement.