Single three-dimensional polymers confined to a slab, i.e., to the region between two parallel plane walls, are studied by Monte Carlo simulations. They are described by N-step walks on a simple cubic lattice confined to the region 1less than or equal tozless than or equal toD. The simulations cover both regions DR-F (where R(F)similar toN(nu) is the Flory radius, with nuapproximate to0.587), as well as the cross-over region in between. Chain lengths are up to N=80 000, slab widths up to D=120. In order to test the analysis program and to check for finite size corrections, we actually studied three different models: (a) ordinary random walks (mimicking Theta polymers); (b) self-avoiding walks; and (c) Domb-Joyce walks with the self-repulsion tuned to the point where finite size corrections for free (unrestricted) chains are minimal. For the simulations we employ the pruned-enriched-Rosenbluth method with Markovian anticipation. In addition to the partition sum (which gives us a direct estimate of the forces exerted onto the walls), we measure the density profiles of monomers and of end points transverse to the slab, and the radial extent of the chain parallel to the walls. All scaling laws and some of the universal amplitude ratios are compared to theoretical predictions. (C) 2004 American Institute of Physics.