We describe an algorithm, based on Euler's method, for solving Volterra
integro-differential equations. The algorithm approximates the relevant
integral by means of the composite Trapezium Rule, using the discrete nodes of
the independent variable as the required nodes for the integration variable. We
have developed an error control device, using Richardson extrapolation, and we
have achieved accuracy better than 1e-12 for all numerical examples considered.Comment: 11 page