We analyze the finite volume scheme for the Riemann-Liouville fractional differential equations. We produce error estimates in an energy norm and
$L_2$ norm for the approximated solution. Numerical results are supplied to justify the theoretical work.