dydt.m

来自「这是本人收集的一些程序源代码」· 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];