Week 3 - Solving second order ODEs

clear all

close all

clc

b = 0.05;

g = 9.81;

l = 1;

m = 0.1;

% intial condition

theta_0 = [0;5];

% time points

t_span = linespace(0,10,500);

% solve ODE

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

plot(t,results(:,1))

hold on

plot(t,results(:,2))

ylabel('plot')

xlabel('time')