Abstract
An importance-sampling Monte Carlo method is used to study N×N Ising square lattices with nearest-neighbor interactions and either free edges or periodic boundary conditions. The internal energy, specific heat, order parameter, susceptibility, and near-neighbor spin-spin correlation functions of the finite lattices are determined as a function of N and extrapolated to the corresponding infinite-system values. The effect of finite size is greater for free edges in all cases. The results agree well with predictions of finite size scaling theory and the shape functions as well as amplitudes of surface contribution terms are determined.