Calculation of particle dynamics in radiationless model of plasma

An implicit numerical scheme for the particle dynamics equations is considered taking into account the specifics of the Darwin model of plasma. The problems of it's consistency, accuracy and stability are considered, sowing the adequacy of transfer of basic physical effects related to the motion of charges in electromagnetic field. On the basis of the scheme proposed an optimal numerical algorithm, considering the problem of it's effectiveness.