Numerical solution of the Fokker–Planck equation with variable coefficients