Week 3 - Solving second order ODEs

clear all
close all
clc

b=0.05
g=9.81
L=1
m=1

%time
t_span=[0,20];

%initial conditions
theta_0=[0;3];

%solve ODE
[t,results]=ode45(p(t,theta)ode_func(t,theta,b,g,L,m),t_span,theta_0);

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