In this paper, cubic B-spline wavelets are used to obtain the numerical solution of Fokker-Plank equations of fractional order, in the Caputo sense. The operational matrices of fractional derivative and integral for cubic B-spline wavelets are constructed and utilized via Galerkin and collocation spectral methods to transform the main problem to some algebraic system. The suggested approach’s accuracy and applicability are demonstrated using an illustrative example.