Numerical Method for The Time Fractional Fokker-Planck Equation