I need a Crank-Nicolson based scheme for a 2D nonlinear Fokker-Planck equation. How can I start? Someone can give me some references?

Thanks
Evo