dydt.m

来自「MATLAB 6.5 从基础讲解其功能」· M 代码 · 共 8 行

function Yd=DYdt(t,Y) 
% t		一定是标量形式的自变量
% Y		必须是列向量
global G ME					%在函数中定义全局变量传递参数
xy=Y(1:2);Vxy=Y(3:4);		%前两个元素是"位移量"，后两个是"速度量"。
r=sqrt(sum(xy.^2));
Yd=[Vxy;-G*ME*xy/r^3];		%Yd必须按式(4.14.2.1-5)编写，是与Y同维的列向量。