damped motion of pendulum using 2nd ODE

Objective: I have attached her the youtube link of the pendulum. Please checkcc. % function for solving a 2nd ODE function[dtheta] = ode(theta,t) b = 0.1; % damping cofficient g = 9.81; % gavity l = 1; % length of string m = 0.5; % mass theta1 = theta(1); theta2 = theta(2); dtheta1 = theta2; dtheta2 = -(b/m)*theta2-(g/l)*sin(theta1);…