International audienceMedical systems usually consider linear propagation of ultrasound, an approximation of reality. However, numerous studies have attempted to accurately simulate the nonlinear pressure wave distortion and to evaluate the contribution of harmonic frequencies. In such simulations, the computation time is very large, except for the method based on the angular spectrum scheme where the derivative order is reduced using the Fourier transform. However, the harmonic computation is usually limited to the second harmonic because of quasi-linear approximation. In this paper, a slowly varying envelope approximation (SVEA) is used in the Fourier domain to compute the entire nonlinear distortion induced, including high harmonics and ...