We propose an integral method to deal with the spherically symmetric n
on-linear Fokker-Planck equation appearing in plasma physics. A probab
ility transition expression is obtained, which takes into account the
proper domain for the radial velocity component. The analytical and co
mputational results are new, and the time evolution is completely sati
sfactory. The main achievement of the method is conservation of both t
he initial norm and energy for unlimited times, which has not been att
ained in the differential approach to the problem.