Using a newly developed density matrix renormalization group algorithm, we have computed exact solutions of the Schrodinger equation for water at two geometries in a basis of 41 orbitals. Calculations of this size cannot be carried out using any other method. We compare our results with high-order coupled cluster and configuration interaction calculations. (C) 2003 American Institute of Physics.