% Anhang B: Nichtlineares Anfangswertproblem mit Daempfung

clear all;

options = odeset('MaxStep' , 0.1); % ... zwingend, um ausreichende Genaigkeit zu erreichen

x0    = [2.5 ; 0] ;   % Anfangswerte  [x ; v]
tspan = [0 100]   ;   % Zeitintervall

% Integration des Anfangswertproblems:
[t x] = ode45 ('NlAWP_Dgl', tspan , x0 , options) ;

% Grafische Ausgabe von x und v, zunaechst "Sortieren der Ergebnisse":

xi = [1 0] * x' ;  % x' ist die transponierte Matrix der Ergebnisse,
vi = [0 1] * x' ;  % xi und vi sind Vektoren 

subplot(3,1,1) ; plot (t  , xi) , grid on , title  ('Weg-Zeit-Diagramm:')
subplot(3,1,2) ; plot (t  , vi) , grid on , title  ('Geschwindigkeits-Zeit-Diagramm:')
subplot(3,1,3) ; plot (xi , vi) , grid on , title  ('Geschwindigkeits-Weg-Diagramm:')
