Recently developed numerical methods have enabled the explicit construction of the superconducting state of the Hubbard model of strongly correlated electrons in parameter regimes where the model also exhibits a pseudogap and a Mott insulating phase. d x 2 −y 2 symmetry superconductivity is found to occur in proximity to the Mott insulator, but separated from it by a pseudogapped nonsuperconducting phase. The superconducting transition temperature and order parameter amplitude are found to be maximal at the onset of the normal-state pseudogap. The emergence of superconductivity from the normal state pseudogap leads to a decrease in the excitation gap. All of these features are consistent with the observed behavior of the copper-oxide superconductors. where ε p = −2t(cos p x + cos p y ) + 4t ′ cos p x cos p y an electron dispersion and U > 0 a local interaction which disfavors double occupancy of a site.In the years since Anderson's paper, the interplay of the pseudogap and superconductivity and the relation of both to the Hubbard model have been of central interest to condensed matter physicists. The existence of dwave superconductivity in the Hubbard model has been demonstrated by perturbative analytic calculations [4] (later improved by renormalization group methods [5, 6]) and by numerics [7, 8]. The issue of the pseudogap has been more controversial. It has been variously argued that the pseudogap is a signature of unusual superconducting fluctuations [9][10][11], of a competing nonsuperconducting phase or regime [3, 12], or of physics not contained in the Hubbard model [13]. Theoretical determination of the interplay of the pseudogap and superconductivity in the Hubbard model is important in helping resolve this controversy, and will provide insight into the pseudogap phenomenon and into strongly correlated superconductivity more generally, but this requires access to intermediate or strong couplings for which perturbation theory is inadequate. The development of cluster dynamical mean field theory [15] has provided important nonperturbative information about the Hubbard model. Dynamical mean field theory approximates the electron self-energy in terms of a finite number of auxiliary functions determined from the solution of an N -site quantum impurity model and becomes exact as N tends to infinity. In this Letter we use dynamical mean field methods to determine the interplay of superconductivity and the pseudogap in the Hubbard model. This is challenging because the theory