Week 3 - Solving second order ODEs
first make ode fuction
function [dtheta_dt] = ode_func(t,theta,b,g,l,m)
theta1 = theta(1)
theta2 = theta(2)
dtheta1_dt= theta2;
dtheta2_dt = -(b/m)*theta2-(g/l)*sin(theta1);
dtheta_dt = [dtheta1_dt;dtheta2_dt];
end
SOLVED MAIN ODE
clear all
close all
clc
b = 0.5;
g = 9.81;
m = 0.1;
l = 1;
% intial condition
theta_0 = [0;5];
% time points
t_span = linspace(0,10,500);
% solve ode
[t_results] = ode45(@(t,theta) ode_func(t,theta,b,g,l,m),t_span,theta_0);