The Darcy's equation consists of the mass conservation equation and the Darcy's law that involves the hydraulic potential (or called pressure) and the fluid velocity, which governs the flow of an incompressible fluid through a porous medium. In this paper, we investigate the Legendre Galerkin spectral collocation least squares method for approximating the problem of Darcy flow in homogeneous medium and non-homogeneous medium, respectively. The proposed scheme can be solved the approximate solutions of the hydraulic potential and the average velocity of the fluid simultaneously. A symmetric positive definite coefficient matrix of the corresponding linear algebra equation is obtained by applying our scheme. Numerical examples are presented to validate the efficiency and accuracy of the proposed scheme.