An effective approach of calculating scattering angle and time integral in Monte Carlo simulation is described for the Ziegler et al. interatomic potential. The scattering angle is determined using the “magic” formula by Lindhard et al. and that by Biersack et al. in the corresponding regions of p/a and ϵ variation. To calculate the time integral, a simple and rather precise approximation formula was obtained.