Self-consistent interaction between charged particles and waves in a plasma is described by a system of two coupled nonlinear (Fokker-Planck) partial differential equations. A finite difference scheme is presented for the numerical solution of this system, and logarithmic mesh points are adopted. The method proves to be efficient and accurate even with a relatively large time step. Test runs are carried out with both zero-flux boundary conditions and zero-value boundary conditions.