from pylab import *
def es(t):
    e = 0.6112*exp(17.62*t/(243.12+t))
    return e
t = arange(0.0, 35.0, 0.5)
plot(t, es(t))
ytext = ylabel('saturated vapour pressure in kPa')
xtext = xlabel('temperature')
show()