second order ODE

clear all

close all

clc

b=0.5;

L= 1;

m=1;

g=9.8;

%initial condition

theta_θ = [0;3];

%time plots

t_span= linspace(0,20,500);

%solve ODE

[t,results]= ode 45 (@(t,theta)ode_func(t,theta,b,l,m,g),t_span,theta_θ);

plot(t,results (:,1))

hold on

plot (t,result (:,2))

ylabel(\'theta\')

xlabel(\'time\')