Solving ODE for Simple Pendulum

Objective: close all clear all clc %Defining Input parameters b=0.5; m=1; g=9.81; l=1; %input conditions theta_0=[0;5]; %time array till 20sec t_span=linspace(0,20,500); [t, results]= ode45(@(t,theta) ode_func(t, theta,b,g,l,m),t_span,theta_0); plot(t,results(:,1)) hold on plot(t,results(:,2)) xlabel(\'time\') ylabel(\'plot\') …