% 単振り子の運動
% 単振り子の運動方程式を数値的に解く

global len;
global mass;
global grav;

len = 2.0;
mass = 0.01;
grav = 9.8;

interval = 0:0.1:10; % 固定ステップ
qinit = [pi/3; 0]; % 初期値
[time, q] = ode45(@simple_pendulum, interval, qinit);

plot(time,q(:,1),time,q(:,2),'--');
title(['simple pendulum : l = ', num2str(len), ', m = ', num2str(mass)]);
xlabel('time');
ylabel('angle \theta, angular velocity \omega');
legend('\theta', '\omega');
