second order ode

Objective: Program clear all close allclc%inputsb=.05;g=9.81;l=1;m=.1;x=0;y=0;ct=1;%cond. 1theta0=[0;5];%timetime=linspace(0,10,500);%calc[t,results]=ode45(@(t,theta)ode_func(t,theta,b,g,l,m),time,theta0);%simulationfor i = 1 :length(results) c = results(i,1); %defining co-ordinates y1 = -l*cos(c); x1 = -l*sin(c); %plot hold on plot([x…