function [] = radpendel_rumpf()
% Konstanten
a  = 1.00;
m1 = 0.10;
m2 = 0.15;

% Anfangs und Endzeitpunkt
t0 =   0.00;
dt =   0.01;
t1 =  10.00;

% Stationärer Punkt
x_stat = % HIER AUSFUELLEN

% Startwerte
x0 = x_stat;
v0 = 0.1;

[T   ,Y   ] = ode45(@myFun,   [t0:dt:t1],[x0 v0]);
[Tlin,Ylin] = ode45(@myFunLin,[t0:dt:t1],[x0 v0]);

figure; plot(T, [Y(:,1) Ylin(:,1)]); grid on;

end



function dY = myFun(t,Y)
% Konstanten
    a  = 1.00;
    r  = 0.20;
    m1 = 0.10;
    m2 = 0.15;
    J  = 0.5*m1*r*r;
    g  = 9.81;

% Nichtlineare DGL
    % HIER AUSFUELLEN
    % dY = 
    % ...
end



function dY = myFunLin(t,Y)
% Konstanten
    a  = 1.00;
    r  = 0.20;
    m1 = 0.10;
    m2 = 0.15;
    J  = 0.5*m1*r*r;
    g  = 9.81;

% Linearisierte DGL
    % HIER AUSFUELLEN
    % dY = 
    % ...
end


