Summary: We investigate strategies for the numerical solution of the initial value problem with initial conditions
where . Here denotes the derivative of order (not necessarily ) in the sense of Caputo. The methods are based on numerical integration techniques applied to an equivalent nonlinear and weakly singular Volterra integral equation. The classical approach leads to an algorithm with very high arithmetic complexity. Therefore we derive an alternative that leads to lower complexity without sacrificing too much precision.