state2.m

弹道仿真程序包

function dy = state2(t,y);

global a b c m_0 d L p g
dy=zeros(2,1);
% y(1)=y  y(2)=dy
dy(1)=y(2);
dy(2)=(b*t-((1+d)*m_0-c*t)*g-0.1*y(2)^2)/((1+d)*m_0-c*t);