% DGL-Funktion für ode45: 
%   für SenkrFall
% function dy_dt = SenkrFall_Fun( t, y, delta )
%    t : Zeitvariable, y : Vektor aus z und v
function  dy_dt = SenkrFall_Fun( t, y, delta )
   g = 9.81;      % Erdbeschleunigung
   
   % dy_dt(1) = Ableitung von z = y(1) ist v = y(2) :
   dy_dt(1,1) = y(2);
   % dy_dt(2) = Ableitung von v = y(2): DGL
   dy_dt(2,1) = - g - delta * y(2);
end  