Abstract: A new SISD iterative algorithm is put forward to effectively solve arbitrary consistent or inconsistent system of linear algebraic equations. It is proved that accurate theoretic or least square solutions of the system can be obtained throughmiterations. In addition, the computational complexity, the numerical stability and the inherent parallelism of the algorithm are discussed.