from pylab import *
def HortonInt(t):
    fu = 3.0
    f0 = 12.0
    k = 0.05
    Hi = fu*t+(f0-fu)/k*(1-exp(-k*t))
    return Hi
t = arange(0.0, 60.0, 0.5)
plot(t, HortonInt(t))
ytext = ylabel('Infiltration (mm)')
xtext = xlabel('Zeit (Min.)')
show()