Simple Pendulum using ODE in MATLAB

Objective: 1)%ode function function [dtheta_dt]=odemain(t,theta,b,g,l,m) theta1=theta(1) theta2=theta(2) dtheta1_dt=theta2 dtheta2_dt=-(b/m)*theta2-(g/l)*sind(theta1) dtheta_dt=[dtheta1_dt;dtheta2_dt] 2)%Plotting and Animation clear all close all clc b=0.06 g=9.81 l=1 m=1.5 %Initial conditions theta_0=[0;3];…