Coulomb collisional relaxation process of ion beams in magnetized plasmas

An orbit following code is developed to calculate ion beam trajectories in magnetized plasmas. The equation of motion (the Newton's equation) is solved including the Lorentz force term and Coulomb collisional relaxation term. Furthermore, a new algorithm is introduced by applying perturbation method regarding the collision term as a small term. The reduction of computation time is suggested.