Week 3 - Solving second order ODEs

clear all
close all
clc

%input variables
b=0.05;
m=0.1;
g=9.8;
l=1;

%initial condition
theta_y=[0;5];

%time span
t_span=linspace(0,10,500);

% solve ODE
[t,results]=ode45(@(t,theta)ode_func(t,theta,l,b,g,m),t_span,theta_y);

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

xlabel('time')
ylabel('plot')