
where \(t\equiv \frac{\vert T-T_c\vert }{T_c}\)  and \(\alpha \)  is the specific heat critical exponent. In order to find the critical exponent, we focus on the region near the critical point and plot \(C_v(T)\)  in this region in FIG. REF . By fitting a curve with the numerical result, the value of the critical exponent is found to be \(\alpha =0.67\) . This value is in complete agreement with the one reported in [1]}, [2]}. To get this number, we also plot the linear log-log diagram for which the critical exponent is the slope of a line i.e. \(\log (C_v)\propto \alpha \log (t)\) . To report how well our result is, we calculate relative error (RE) and root mean square (RMS) which are defined as
\(&{\rm RE}=\frac{\alpha -2/3}{2/3}, \cr &{\rm RMS}=\sqrt{\frac{1}{N}\sum \limits _{i=1}^N\left(y_{\rm fitted}(i)-y_{\rm data}(i)\right)^2},\) 
