The practical method is proposed in this papers which is based on the iterative method for solving linear system of equations. The normal equations are solved directly by the algorithm and a minimum norm solution can be obtained when coefficient matrix of normal equation is singular, so that the sequence generated by the algorithm may be convergent. The numerical examples illustrate the algorithm is very efficient for singular problem or ill-condition problem.