monte carlo simulation of square ice model for residual entropy calculation(thermodynamic properties)