% Stabpendel mit dem ode45-Solver von MATLAB 

global pl g
pl     = 1 ;                               % Pendellaenge
g      = 9.81 ;                            % Erdbeschleunigung

t0   =  0 ;  phi0 = pi/2 ; omega0 = 0 ;    % Anfangswerte
tEnd = 10 ;                                % Ende des Integrationsintervalls

% MATLAB-Solver ode45 ------------------------------------------------------ 
[t , x] = ode45 ('BewDglPendel' , [t0 ; tEnd] , [phi0 ; omega0]) ;
% MATLAB-Solver ode45 ------------------------------------------------------

phi   = x(:,1) ;                           % Sortieren der Ergebnisse:
omega = x(:,2) ;                           % Spalte 1: phi, Spalte 2: omega

phiMin   = min(phi) 
omegaMax = max(abs(omega))                 % Extremwerte ---> Command Window
plot (t,phi)                               % Graphik-Ausgabe phi(t)