% súbor ‘kyv23sol.m’ na riešenie DR kyvadla
Tspan = [1 5]; %t0 = 0; tfin = 5; èasovy interval na výpoèet
x0 = [1;0]; % poèiatoèné podmienky
opt=odeset('AbsTol',1,'RelTol',1,'InitialStep',0.1,'MaxStep',0.1);
[t,x] = ode23('kyvadlo',Tspan, x0, opt);
subplot(211), plot(t,x(:,1), ':',t,x(:,2),'-'); title('Kyvadlo'),grid
xlabel('t[s]'), ylabel('uhol [rad], rych [rad/s]');
subplot(212), plot(x(:,1), x(:,2)); grid