Numerical solution of fractional order diffusion problems with Neumann boundary conditions