Este documento presenta un método numérico práctico, un esquema implícito de diferencias finitas para resolver una ecuación de Fokker-Planck fraccional en tiempo-espacio bidimensional con coeficientes variables y término fuente, que representa un modelo de una partícula browniana en un potencial periódico. Se consideran la derivada de Caputo y la derivada de Riemann-Liouville en las direcciones temporal y espacial, respectivamente. La derivada de Riemann-Liouville se aproxima mediante la aproximación estándar de Grünwald y la aproximación de Grünwald desplazada. Se discuten la estabilidad y convergencia del esquema numérico. Finalmente, se proporciona un ejemplo numérico para probar el análisis teórico.