Using an improved version of the projection quantum Monte Carlo technique, we study the square-lattice Hubbard model with nearest-neighbor hopping t and next-nearest-neighbor hopping t 0 by simulation of lattices with up to 20 3 20 sites. For a given R 2t 0 ͞t, we consider that filling which leads to a singular density of states of the noninteracting problem. For repulsive interactions, we find an itinerant ferromagnet (antiferromagnet) for R 0.94 (R 0.2). This is consistent with the prediction of the T -matrix approximation, which sums the most singular set of diagrams. [S0031-9007 (97)02456-3] PACS numbers: 75.10.Jm, 71.10.Fd, 71.27.+ a The understanding of itinerant ferromagnetism (FM) is a long-standing problem of solid-state physics [1]. In the search for a generic model of FM, the Hubbard model, describing electrons from a single band subject to a local electron-electron repulsion U, has been investigated extensively. Motivated by an exact result of Nagaoka [2], most papers studied the stability of a fully polarized state (believed to occur in the phase diagram of the Hubbard model at large interaction strength and close to half filling) against single spin flips. The results turned out to be strongly dependent on the quasiparticle spectrum [3], and, generically, unrealistically large U were required to stabilize a fully polarized state. Motivated by the recent proofs of FM in certain models with flat bands [4], in this Letter we take a complementary route, investigating the stability of the paramagnetic phase of a model with a high density of states against FM at weak coupling.We consider electrons on a square lattice with L l 3 l (even l) sites, described by the Hamiltonian H 2t X c y i,s c j,s 1 t 0