Week 3 - Solving second order ODEs

Objective: %solution of odeclose allclear allclc%inputb=0.5; %damping coefficient.g=9.81; %accelaration due to gravity in m/s^2l=1; % length of the pendulum in mm=1; %mass of the ball in kg,%initial conditiontheta_0=[0;3];%time limitt_span=linspace(0,20,500);%solving ode[t,results]=ode45(@(t,theta) ode1_fun(t,theta,b,g,l,m), t_span,theta_0);ang_disp=results(:,1);ang_vel=results(:,2);figure(1)hold…