Week 3 - Solving second order ODEs

cleaer all

close all

clc

b = 0.05;

g = 9.81;

l= 1;

m= 0.1;

%initial condition

theta_0 = [0;5];

%time points

t_span = linspace(0,10,500);

%solve ODE

[t,result] = ode45(@(t,theta) ode_func(t,theta, b,g,l,m),t_span, theta_0);

plot (t ,result (:,1))

hold on

plot (t,result (:,2))

ylabel('plot')

xlabel('time')