ode

Objective: clear all close all clc b=0.05; g=9.81; l=1; m=0.1; %initial condition theta_0=[0;5]; %time points t_span=linspace(0,10,500); %solve ode [results]=lsode(\'ode_fun\',theta_0,t_span); plot(t,results(:,1)) hold on plot(t,results(:,2)) ylabel(\'plot\') xlabel(\'time\')