import numpy as np
import matplotlib.pyplot as plt

x = np.loadtxt('convergence.txt')
ind = np.arange(0,len(x[:,0]))

plt.plot(ind,1e6*x[:,2])
plt.xlabel('Iteration')
plt.ylabel('RMS emittance (mm mrad)')
plt.savefig('convergence_rough.eps')
plt.show()