dydt.m

这是本人收集的一些程序源代码

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];