Rep:HeatCapacityVsTemp.py

import numpy as np

from matplotlib import pylab as pl

import matplotlib.pyplot as plt

def plot_HC_vs_T(x):
    data = np.loadtxt(str(x) +'x' + str(x) + '.dat')
    
    T = data[:,0]
    
    E = data[:,1]
    E2 = data[:,2]

    VarE = E2 - E**2

    C = (VarE/(T**2*x*x))
    
    return T, C

for i in range (1,6):
    T, C = plot_HC_vs_T(2**i)

    pl.plot(T, C, label = str(2**i) + 'x' + str(2**i))
    plt.legend()
    pl.title('Heat Capacity vs Tempurature')
    pl.xlabel('Tempurature')
    pl.ylabel('Heat Capacity')
    
pl.show()