In this paper we examine $ \theta-$ method to solve Fractional Poisson Differential Equations $ (0\leq\theta\leq1) $. We use Riemann-Liouville derivative based on Shifted Grunwald estimates for fractional derivative. Consistency , stability and (therefore) convergence of the method is discussed. A numerical example is presented and compared with the exact analytical solution.
