A recursive algorithm for the approximate solution of Volterra integral equations of the first kind of convolution type