from pylab import *

def interception(No, LAI):
    icm = 0.935+0.498*LAI-0.00575*LAI*LAI
    ic = -0.42+(0.245*No)+(0.2*LAI)+(0.0271*No*LAI)-(0.0111*No*No)-(0.0109*LAI*LAI)
    return ic
    
LAI = 3.0    
No = arange(0.0, 25.0, 1.0)
plot(No, interception(No,LAI))
xtext = xlabel('precipitation (mm/d)')
ytext = ylabel('Interception (mm)')
setp(xtext, size='medium', name='courier', weight='bold', color='g')
setp(ytext, size='medium', name='helvetica', weight='light', color='b')
show()