Eqns - velocity = x(2 mass = x(3 gravity = g_zero(r_earth(altitude r_earth(r_earth(altitude r_earth density = rho_zero*exp-altitude/scale_hgt drag

Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
function [dxdt] = eqns(n, t, x) global g_sealevel %reference value sea level global r_earth %radius of earth global rho_zero %density at sea level global scale_hgt %scale height global s_ref %reference area global cd %drag coefficient global thrust %thrust global m_dot %rate of mass change two = 2.0; altitude = x(1);
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: velocity = x(2); mass = x(3); gravity = g_zero*(r_earth/(altitude + r_earth))*(r_earth/(altitude + r_earth)); density = rho_zero*exp(-altitude/scale_hgt); drag = density*velocity*velocity*refa*cd/two; dxdt(1) = velocity; dxdt(2) = (thrust - drag)/mass - gravity; dxdt(3) = m_dot;...
View Full Document

This note was uploaded on 07/12/2011 for the course ECON 011 taught by Professor Bob during the Spring '11 term at Tufts.

Ask a homework question - tutors are online