  We have simulated the three-dimensional Heisenberg model on simple cubic
lattices, using the single-cluster Monte Carlo update algorithm. The expected
pronounced reduction of critical slowing down at the phase transition is
verified. This allows simulations on significantly larger lattices than in
previous studies and consequently a better control over systematic errors. In
one set of simulations we employ the usual finite-size scaling methods to
compute the critical exponents $\nu,\alpha,\beta,\gamma, \eta$ from a few
measurements in the vicinity of the critical point, making extensive use of
histogram reweighting and optimization techniques. In another set of
simulations we report measurements of improved estimators for the spatial
correlation length and the susceptibility in the high-temperature phase,
obtained on lattices with up to $100^3$ spins. This enables us to compute
independent estimates of $\nu$ and $\gamma$ from power-law fits of their
critical divergencies.
