Specific heat
There appears to be a problem with the specific heat in our code. When we run a range of values for the parameter 1<\tau<4
, we obtain the following plots for the susceptibility and magnetisation as expected.
The specific heat however, does not look as it should. We calculated it with: Cv[i] = 1/(kB*beta**2*L**2)*np.var(E)
. This results in:
Is the reason for this that we should not use the variance in E in this way or might the error be elsewhere?