In this paper, the parabolic approximation of wave equation will be solved by the method of least squares. At first, the radio wave propagation in homogeneous media will be considered. The electromagnetic field will be expanded by proper expansion functions, which satisfy the parabolic equation in homogeneous media. The expansion coefficients will be derived by the least square method through enforcing initial and boundary conditions. The least square functional satisfy the initial and boundary conditions.