We solve the dynamical mean-field theory equations for the Hubbard model away from the particle-hole-symmetric case using the density matrix renormalization group method. We focus our study on the region of strong interactions and finite doping where two solutions coexist. We obtain precise predictions for the boundaries of the coexistence region. In addition, we demonstrate the capabilities of this precise method by obtaining the frequency-dependent optical conductivity spectra.7512