% s: Anfangswert, tm: gewünschte Zeitdauer
% löst DGL D1
function  DglSolver( s, tm )
  
   tspan = [0,tm];  % Zeitintervall
   x0 = s;          % Anfangswert

   % Solver-Aufruf, Übergabe der DGL-Funktion
   [t,x] = ode45( @D1, tspan, x0 );

   % Ergebnis plotten
   plot( t, x );
   grid on;
end