Week 4.1 - Solving second order ODEs

l=1;

m=1

b=0.05;

g=9.81;

t_span=linspace(0,20,100);

theta0=[0;3];

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

plot(t,results)