% Stabpendel mit Runge-Kutta, Aufruf der Function rk4, der fast identisch 
% mit dem Aufruf der MATLAB-ode-Functions

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

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

% Runge-Kutta-Algorithmus -------------------------------------------------- 
[t , x] = rk4 ('BewDglPendel' , [t0 ; tEnd] , [phi0 ; omega0] , 200) ;
% Runge-Kutta-Algorithmus --------------------------------------------------

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)