function hd = Plinius(t,h);

R      = 0.05;          % Radius des Bodens
r      = 0.005;          % Radius des Ausguss
hlo    = 0.025;         % Höhe der Ausflussöffnung
hhi    = 2;           % Höhe des Siphon
C      = 0.6;           % Abfluss Koeffizient
zulauf = 50e-6;         % Deciliter/s
g      = 9.81;          % Gravitationskonstante
global syphon           % Deklaration der global Funktion syphon

if h>=hhi,              %voller Brunnen
  syphon=1;
end;
                        
if h<=hlo,              %leerer Brunnen
  syphon=0;
end;

if syphon,
   v =C*sqrt(2*g*h);      %Funktion für den Ausfluss
   ausfluss = (-pi)*r^2*v;
else
   ausfluss = 0;        %kein Ausfluss
end;
hd =(zulauf+ausfluss)/(pi*R^2); %Stand des Wasserpegels zur Zeit t