Numerical studies of the two-dimensional Hubbard model have shown that it exhibits the basic phenomena seen in the cuprate materials. At half-filling one finds an antiferromagnetic Mott-Hubbard groundstate. When it is doped, a pseudogap appears and at low temperature d-wave pairing and striped states are seen. In addition, there is a delicate balance between these various phases. Here we review evidence for this and then discuss what numerical studies tell us about the structure of the interaction which is responsible for pairing in this model. stripes, pseudogap behavior, and d x 2 −y 2 pairing. In addition, the numerical studies have shown how delicately balanced these models are between nearly degenerate phases. Doping away from half-filling can tip the balance from antiferromagnetism to a striped state in which half-filled domain walls separate π-phase-shifted antiferromagnetic regions. Altering the next-near-neighbor hopping t ′ or the strength of U can favor d x 2 −y 2 pairing correlations over stripes. This delicate balance is also seen in the different results obtained using different numerical techniques for the same model. For example, density matrix renormalization group (DMRG) calculations for doped 8-leg t-J ladders find evidence for a striped ground state.[12] However, variational and Green's function Monte Carlo calculations for the doped t-J lattice, pioneered by Sorella and co-workers, [23,24] find groundstates characterized by d x 2 −y 2 superconducting order with only weak signs of stripes. Similarly, DMRG calculations for doped 6-leg Hubbard ladders [14] find stripes when the ratio of U to the near-neighbor hopping t is greater than 3, while various cluster calculations [27,[30][31][32][33] find evidence that antiferromagnetism and d x 2 −y 2 superconductivity compete in this same parameter regime.These techniques represent present day state-of-the-art numerical approaches. The fact that they can give different results may reflect the influence of different boundary conditions or different aspect ratios of the lattices that were studied. The n-leg open boundary conditions in the DMRG calculations can favor stripe formation. Alternately, the cluster lattice sizes and boundary conditions can frustrate stripe formation. It is also possible that these differences reflect subtle numerical biases in the different numerical methods. Nevertheless, these results taken together show that both the striped and the d x 2 −y 2 superconducting phases are nearly degenerate low energy states of the doped system. Determinantal quantum Monte Carlo calculations [21] as well as various cluster calculations show that the underdoped Hubbard model also exhibits pseudogap phenomena. [27-32] The remarkable similarity of this behavior to the range of phenomena observed in the cuprates provides strong evidence that the Hubbard and t-J models indeed contain a significant amount of the essential physics of the problem. [34]In this chapter, we will focus on the one-band Hubbard model. Section 2 provides an