function T11_Beispiel_Punktpendel
close all
% Parameter
global g l
%m = 10;
%c = 100000;
%d = 70;
g = 9.81;
l = 1.;

t_end = 6; phi0 = pi/14; phi0_p = 0;

tic; [t,z]=ode45(@dgl,[0 t_end],[phi0 phi0_p],[]); zeit_runge_kutta=toc;
plot(t,z(:,1),'-*',t,z(:,2),'-o'); grid on
legend('Lage','Geschwindigkeit')
title(['ode45 mit Schrittweitenkontrolle, Rechenzeit: ' num2str(zeit_runge_kutta)])

tic; [t_l,z_l]=ode45(@dgl_lin,[0 t_end],[phi0 phi0_p],[]); zeit_runge_kutta=toc;
figure
plot(t_l,z_l(:,1),'-*',t,z(:,1),'-o'); grid on
legend('linearisiert','original dgl')
title(['ode45 mit Schrittweitenkontrolle, Rechenzeit: ' num2str(zeit_runge_kutta)])

end

function z_p=dgl(t,z)
   global g l
   z_p=[z(2);...
        -g/l*sin(z(1))];
end

function z_p=dgl_lin(t,z)
   global g l
   z_p=[z(2);...
        -g/l*z(1)];
end
