The numerical results show that the approach is easy to implement and straightforward when applied to space-time fractional Fokker-Planck equations. The method provides a promising tool for solving space-time fractional partial differential equations.