Solving 2nd order ODE using MATLAB and simulating the motion of pendulum

Objective: Link to youtube video of pendulum motion : https://youtu.be/NMAZ9e2sX0g function tdot = pendulum_ode(t,y) tdot = zeros(2,1); tdot(1) = y(2); L=1; m=1; b=0.5; g=9.81; tdot(2) = -((b*y(2)/m) + (g*sin(y(1))/L)); end close all clear all clc t = [0 20]; y_0 = [0 3]; [t, y] = ode45(\'pendulum_ode\',t,y_0); theta…