a new title

Objective: clear allclose allclcb = 0.05g = 9.81l = 1m = 1theta_0 = [0;5]t_span = linspace(0,20,500)[t,results] = ode45( @(t,theta) ode_func1(t,theta,b,g,l,m),t_span, theta_0);subplot(5,4,[13,14,17,18]);plot(t,results(:,1),'color','b');hold onplot(t,results(:,2),'color','r');hold offsubplot(5,4,[15,16,19,20]);plot(results(:,1), results(:,2));xlabel('position-Radions');ylabel('Angular…