Week 3 - Solving second order ODEs

clear 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, results] = ode45(@(t,theta)ode_func(t,theta,b,g,l,m),t_span,theta_0);

plot(t,results(:,2))
hold on
plot(t,results(:,1))
ylabel('plot')