Abstract
We use a histogram Monte Carlo simulation method to calculate the scaling functions of the existence probability Ep and the percolation probability P of the site percolation model on square lattices with free and periodic boundary conditions. We find that different boundary conditions give quite different scaling functions near the critical region. However, they give the consistent critical point, critical exponents, and the thermodynamic order parameter from renormalization-group calculations. Similar results are found for other percolation models. The implications of our calculated results for some theoretical problems of current interest are discussed.