Abstract
In this paper a least square method is proposed for the solution of the 3-dimensional (3-D) vector parabolic wave equation. This method can correctly incorporate the impedance boundary conditions of the propagation environment. A proper basis function is introduced to expand electromagnetic field components in the successive computation steps. The proposed method needs much less memory size against its finite-difference counterpart.