The demand for high performance computer grows to solve large linear systems of equations in such engineering fields - circuit simulation for VLSI design, image processing, structural engineering, aerodynamics, etc. Many various parallel processing systems have been proposed and manufactured to satisfy the demand. The properties of linear system determine what algorithm is proper to solve the problem. Direct methods or iterative methods can be used for solving the problem. In this paper, an intelligent parallel iterative method for solving linear systems of equations with large sparse matrices is proposed and its efficiency is proved through simulation.