The results show that the method used in the paper is better than the traditional Krylov subspace methods including BICG, CGS, BICGSTAB, GMRES and QMR.