Rep:DMS3053-7

CMP Programming Experiment - Daniel Spencer, 00736964
Python script for Determining the heat capacity

   from matplotlib import pylab as pl 

   import numpy as np

   def extract_data(file):
       data_array = np.loadtxt(file)
       temps = data_array[:,0]
       energies = data_array[:,1]
       squared_energies = data_array[:,2]
       variance = squared_energies - (energies*energies)
       heat_capacity = variance / (temps*temps)
       return temps, heat_capacity
   
   def plotdata(file1, file2, file3, file4, file5):
       file1_data = extract_data(file1)
       file2_data = extract_data(file2)
       file3_data = extract_data(file3)
       file4_data = extract_data(file4)
       file5_data = extract_data(file5)
       pl.ylabel("Heat Capacity")
       pl.xlabel("Temperature")
       pl.ylim([0.0, 1.25])
       pl.plot(file1_data[0], file1_data[1], file2_data[0], 
           file2_data[1], file3_data[0], file3_data[1],
           file4_data[0], file4_data[1], file5_data[0], file5_data[1])
       pl.legend(['2x2', '4x4', '8x8', '16x16', '32x32'], loc='upper center', bbox_to_anchor=(0.5, 1.10),
             ncol=5)
       pl.show()