A new code for Volterra integral equations based on natural Runge-Kutta methods

Abstract We describe some issues related to the development of a new code for numerical solution of systems of Volterra integral equations of the second kind. This code is based on A - and V 0 -stable natural Volterra Runge-Kutta method of order p = 4 implemented in a variable stepsize environment. This method was derived by Conte et al. (2014) [13] . The numerical experiments on many equations and systems of equations illustrate that the code achieves the expected order of accuracy, and confirm its efficiency and robustness.