from pylab import *
def PhilippInt(t):
    S = 85.0
    K = 5.0
    Ip = S*t**(0.5)+K*t
    return Ip
t = arange(0.0, 60.0, 0.5)
plot(t, PhilippInt(t))
ytext = ylabel('Infiltrationsmenge (mm)')
xtext = xlabel('Zeit (Min.)')
show()