A numerical method for solving inhomogeneous nonlocal time fractional convection-diffusion-reaction equations with variable coefficients has been developed. The fractional time operator is taken in the sense of the modified ...