a new title

Objective: % ode function(b,g,l,m,theta1,theta2) % ode_func is a function of t,theta,b,g,lmfunction [dtheta_dt] = ode_func(t,theta,b,g,l,m)theta1 = theta(1)theta2 = theta(2)dtheta1_dt = theta2; dtheta2_dt = -(b/m)*theta2 - (g/l)*sin(theta1);dtheta_dt = [dtheta1_dt; dtheta2_dt];end% writing script for motion of simple pendulum.clear…