import numpy as np
import pylab as pl

pl.rc('lines', linewidth=2)

def s2a(s):
    return np.array(map(float, s.replace(',', '.').split()))

pl.plot(s2a("""-4,74
-3,69
-3,47
-3,49
-3,41
-3,4
"""), color='b', label='O fcc Pt')

pl.plot(s2a("""-2,31
-1,68
-1,14
-1,17
-1,14
-1,14
"""), color='g', label='CO 111t Pt')

pl.plot(s2a("""-2,39
-1,88
-1,71
-1,78
-1,69
-1,67
"""), color='r', label='CO edge Pt')

pl.plot(s2a("""-3.41
-2.02
-2.69
-2.52
-2.34
-2.3
-2.34"""), color='b', ls='--', label='O fcc Au')

pl.plot(s2a("""-0.68
-0.21
0
0
0.09
0.13
0.07"""), color='g', ls='--', label='CO 111t Au')

pl.plot(s2a("""-0.74
-0.35
-0.36
-0.37
-0.26
-0.33
-0.31"""), color='r', ls='--', label='CO edge Au')


Nvalues = [13, 55, 147, 309, 561, 923, 1415]
#    pl.plot(E, lw=2, label=descriptions[i])
pl.xticks(range(len(Nvalues)), map(str, Nvalues))
pl.xlabel('Number of atoms in cluster')
pl.ylabel('Adsorption energy [eV]')
pl.legend(loc='lower right')

pl.show()