Solving second order ODEs using MATLAB

Objective: Objective- To simulate the transient behaviour of a simple pendulum create an animation of it's motion MATLAB Code- function [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 b = 0.05;…