We propose an efficient numerical algorithm for solving integral equations of the theory of liquids in the reference interaction site model (RISM) approximation for infinitely dilute solution of macromolecules with a large number of atoms. The algorithm is based on applying the nonstationary iterative methods for solving systems of linear algebraic equations. We calculate the solvent-solute atom-atom correlation functions for a fragment of the B-DNA duplex d(GGGGG).d(CCCCC) in infinitely dilute aqueous solution. The obtained results are compared with available experimental data and results from computer simulations.