%Ausführen der Differentialgleichungs Berechnung und anschliessendes
%Zeichnen des Füllstandverlaufes h nach der Zeit t

global syphon
syphon = 0;
t0=0; tf=10000;
h0=0;
[t,h]=ode45(@plinius, [t0, tf],h0);
plot(t,h)
clear global syphon