Week 4.1 - Solving second order ODEs
clear all
close all
clc
b = 0.05;
l = 1;
m = 1;
g = 9.81;
ct = 1;
theta_0 = [0;3];
t_span = linspace(0,20,500);
[t,results] = ode45(@(t,theta) ode_func(t,theta,b,g,l,m),t_span,theta_0);
for i=1:500
plot(t(i),results(i,1),'-o')
hold on
pause(0.03)
M(ct) = getframe(gcf);
ct = ct+1;
end
movie(M)
videofile = VideoWriter('pendulum_motion.avi','Uncompressed AVI');
open(videofile)
writeVideo(videofile,M)
close(videofile)