% 4-5-2008 % Higher order ODE example clear all; close all; % Parameter values for our ODE A = 2; % ODE solver parameters y0 = [0 1 0]; % Initial conditions tspan = [0 30]; % time interval we would like to solve for % Solve our ode using 4th Order Runge-Kutta [t,y] = ode45('high_order_rhs',tspan,y0,[],A); plot(t,y(:,1)); title('ODE Solution'); xlabel('t'); ylabel('y(t)');