Porovnání fyzikálního a matematického kyvadla: Matematické kyvadlo:

d^2 theta/ dt^2 = -g . theta/l

Fyzikální kyvadlo:

d^2 theta/ dt^2 = -g . sin(theta)/l

g = 9.8; l = 2;
f = @(t,theta) [theta(2); -g*sin(theta(1))/l]; 
theta_ic = [pi/3; 0];
tspan = linspace(0,10,500);
[t, theta] = ode45(f, tspan, theta_ic);
plot(t, theta);
legend({'$\theta - fyz$', '$\dot{\theta} - fyz$'}, ...
 'Location', 'best', ...
 'Interpreter', 'latex', ...
 'FontSize', 16)
hold on

f = @(t,theta) [theta(2); -g*theta(1)/l];
[t, theta] = ode45(f, tspan, theta_ic);
plot(t, theta);