Variational perturbation theory for Fokker-Planck equation with nonlinear drift

We develop a recursive method for perturbative solutions of the Fokker-Planck equation with nonlinear drift. The series expansion of the time-dependent probability density in terms of powers of the coupling constant is obtained by solving a set of first-order linear ordinary differential equations. Resumming the series in the spirit of variational perturbation theory we are able to determine the probability density for all values of the coupling constant. Comparison with numerical results shows exponential convergence with increasing order.