Abstract
We study the two-dimensional Hubbard model on a finite lattice. By numerically iterating the Hartree-Fock equations one obtains baglike polaron solutions for nearly half-filled-band cases. The profile of the polaron and the interaction potential between the polarons are calculated. Comparison with analytical results of Schrieffer, Wen, and Zhang is made. We also discuss a scenario in which the polarons are turned into the superconducting state by their mutual attraction.