Fokker-Planck equations can be efficiently solved using the Chang-Cooper method combined with unconditionally positive and conservative Patankar-type time integration schemes, which preserve positivity and steady states.
