A new density matrix renormalization group (DMRG) algorithm is presented wh
ich conserves the total spin. This is the first time that standard DMRG is
extended to exploit fully the symmetries of the Hamiltonian. The obtained g
eneral procedure works for any operator, even operators that do not commute
with the Hamiltonian. The new method gives substantial improvements over t
he standard DMRG.