We have studied helix-coil transition in all-atom model of polyalanine in vacuum and distance-dependent dielectric function solvent model. Molecules of up to length 30 residues were investigated by Monte Carlo method. The single histogram method was performed to obtain thermodynamic quantities such as average energy and specific heat. We have also calculated probability distribution of energy and Binder cumulant as well as autocorrelation times for both case in order to identify the order of phase transition by using finite-size scaling arguments. Even though the probability distributions showed a single Gaussian peak, especially for this solvent model, the significant differences between the vacuum and solvent model were observed at their specific heat maximum values and Binder cumulant behaviors. While Binder cumulant has a deep valley for vacuum, it is not sensitive for the distance-dependent dielectric function solvent model. Our data are consistent with recently obtained results for polyalanine by multicanonical method. (C) 2006 Elsevier B.V. All rights reserved.