% Seite 500: Nichtlineares Anfangswertproblem

clear all;

x0    = [0 ; 0] ;   % Anfangswerte  [phi ; omega]
tspan = [0  10] ;   % Zeitintervall

% Integration des Anfangswertproblems (definiert im "Function File" S500_Dgl.m):
[t x] = ode45 ('S500_Dgl', tspan , x0) ;

% Grafische Ausgabe von phi und omega, zunaechst "Sortieren der Ergebnisse":

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

subplot(2,1,1) ; plot (t  , phii)   , grid on , title  ('phi-t-Diagramm:')
subplot(2,1,2) ; plot (t  , omegai) , grid on , title  ('omega-t-Diagramm:')

Vektorlaenge    = length(t)
ErsteZeitpunkte = t(1:12)

phiEnd   = phii  (Vektorlaenge)
omegaEnd = omegai(Vektorlaenge)
