% Stabpendel mit großen Ausschlägen

clear all;

x0    = [3*pi/4 ; 0] ; % Anfangswerte  [phi0 ; omega0]
tspan = [0  10] ;      % Zeitintervall

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

%options = odeset ('MaxStep' , 0.1) ;
%[t x] = ode45 ('Stabpendel_Dgl', tspan , x0 , options) ;

% 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)

phiMin   = min (phii)
omegaMax = max (abs(omegai))
