E2_MPvsHB.py

import matplotlib.pyplot as plt
import seaborn
import numpy as np
data = np.genfromtxt("E2_HB.dat")
x=data[:,0]
y=data[:,1]
yerr=data[:,2]
plt.errorbar(x,y,yerr=yerr,fmt='x')
plt.scatter(x,y,marker='x',label="heatbath")

data = np.genfromtxt("E2_MP.dat")
x=data[:,0]
y=data[:,1]
yerr=data[:,2]
plt.errorbar(x,y,yerr=yerr,fmt='x')
plt.scatter(x,y,marker='x',label="Metropolis")
plt.xlim(0,)
plt.xlabel("T")
plt.ylabel(r"$E^2$")

plt.legend()
plt.show()