Week 3 - Solving second order ODEs

L=1

m=1

b=0.05

g=9.81

%initial condition

theta_@=[0;5];

%Time points

t_span=linespace(0,10,500);

%solveode

[t,results]=ode30(@(t,theta)ode_function(t,theta,b,g,l,m)t,t_span,theta_@);

plot(t_result(:,1))

holdon

plot(t,result(:,2))