说明：  用Monte-Carlo方法研究二维Ising模型的相变问题，求出临界温度T，比热C及磁化率χ
(With Monte-Carlo method of two-dimensional Ising model of phase-change problem, find critical temperature T, specific heat C and susceptibility χ)