Abstract: A linear algorithm of simulation of high power laser beam propagation in nonlinear medium is given, based on the BespalovTalanov linearization theory of nonlinear Schrdinger equation. It is faster in 2～3 orders than splitstep Fourier method and more precise while the initial beam distribution is approximately uniform in near field propagation.