Fast Runge-Kutta methods for nonlinear convolution systems of volterra integral equations