We determine the scaling functions describing the crossover from Ising-like
critical behavior to classical critical behavior in two-dimensional systems
with a variable interaction range. Since this crossover spans several decades
in the reduced temperature as well as in the finite-size crossover variable, it
has up to now largely evaded a satisfactory numerical determination. Using a
new Monte Carlo method, we could obtain accurate results for sufficiently large
interactions ranges. Our data cover the full crossover region both above and
below the critical temperature and support the hypothesis that the crossover
functions are universal. Also the so-called effective exponents are discussed
and we show that these can vary nonmonotonically in the crossover region.Comment: 24 pages RevTeX 3.0/3.1, including 22 PostScript figures. Uses
epsf.st