In order to improve the convergence properties of the conjugate gradients-squared method, the author suggests making use of preconditioned matrices. Four examples of symmetric and nonsymmetric matrices are considered where modified incomplete Cholesky decomposition or incomplete LU-factorization is used. The calculation results show that the residual in this case decreases considerably faster and that this decrease is almost monotonous.