CPlusPlusVsMyData.py

import numpy as np

from matplotlib import pylab as pl

import matplotlib.pyplot as plt

def get_data_HC_vs_T(x):
    data = np.loadtxt(str(x) +'x' + str(x) + '.dat')
    data_C = np.loadtxt('H:/ImperialChem-Year3CMPExpt-master/C++_data/' + 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))
    
    T_C = data_C[:,0]

    C_C = data_C[:,5]
    
    return T, C, T_C, C_C
    

T, C, T_C, C_C = get_data_HC_vs_T(2)
    
pl.plot(T, C, label = str(2) + 'x' + str(2))
pl.plot(T_C, C_C, label = str(2) + 'x' + str(2) +'C++ data')

plt.legend()
pl.title('Heat Capacity vs Tempurature')
pl.xlabel('Tempurature')
pl.ylabel('Heat Capacity')

pl.show()