pendulum

Objective: clear all close all clc %initial condition l=1 theta_0=[0 5]; %time points5 t_span=linspace(0,10,100); %solve ode [t,results] = ode45(@(t,theta) ode_func(t,theta),t_span,theta_0); ct=1; for i=1:length(t_span) j=results(i,1) a0=0; b0=0; a1=l*sind(j); b1=-l*cosd(j); figure(2) plot([a0,a1],[b0,b1],'marker','o','markersize',5,'markerfacecolor','r');…