function [] = pendel()
% Anfangs und Endzeitpunkt
t0 =   0.000;
dt =   0.001;
t1 =   6.000;

% Startwerte
phi0   = 0;
Omega0 = pi/4;

[T,Phi] = ode113(@myFun,[t0:dt:t1],[phi0 Omega0]);
% Winkel in Grad umberechnen
% 0° entspricht der Lage Senkrecht nach oben
Phi(:,1) = mod(180+180*Phi(:,1)/pi, 360)-180;
Phi(:,2) = 180*Phi(:,2)/pi;

figure; plot(T, Phi(:,1));
ylim([-180 180]);
set(gca,'YTick',[-180:30:180]);
grid on;

end


function dY = myFun(t,Y)
% Konstanten
    g  = 9.81;
    m  = 1.00;
    l  = 0.50;
    r  = 0.25; % entspricht 1/2 * l
    b  = 0.05;
    JS = m*(l*l+b*b)/12;
    JA = JS+m*r*r;

% Nichtlineare DGL
    dY = zeros(2,1);
    dY(1) = Y(2);
    dY(2) = ((m*g*r)/JA) * sin(Y(1));
end

