Angular and energy distributions due to multiple small angle scattering were calculated with different models, namely from the analytical Szilágyi theory, the Monte Carlo code MCERD in binary collision approximation, and the molecular dynamics code MDRANGE, for 2 MeV 4 He in Au at backscattering geometry and for 20 MeV 127 I recoil analysis of carbon. The widths and detailed shapes of the distributions are compared, and reasons for deviations between the different models are discussed.