Abstract
A detailed scattering model, based on the continuous slowing down and the Born approximations, is shown to give results in good accordance with experiments when used in conjunction with the Monte Carlo method. Relativistic elastic and inelastic cross-sections adopted are derived from an analytical fit to the Dirac-Hartree-Fock-Slater atomic density with Wigner-Seitz boundary conditions. Random tracks are simulated in a completely analytical way; this fact largely reduces the complexity of the Monte Carlo procedure.