from pylab import * def f(t): d = 0.85 k= 1.9 et = (0.142*t+1.095)*(t+17.8)*k*d return et t = arange(0.0, 35.0, 0.5) plot(t, f(t)) ytext = ylabel('Evaporation') xtext = xlabel('temperature') setp(xtext, size='medium', name='courier', weight='bold', color='g') setp(ytext, size='medium', name='helvetica', weight='light', color='b') show()