dydt4.m

精通matlab

function varargout=DYDt4(t,Y,flag,G,ME,tspan,Y0)
switch flag
case ''
   varargout{1} = f(t,Y,G,ME);		
case 'init'  
  [varargout{1:3}] = fi(tspan,Y0);
otherwise
  error(['Unknown flag ''' flag '''.']);
end
% ------------------------------------------------------------------
function Yd = f(t,Y,G,ME)
X=Y(1:2);V=Y(3:4);r=sqrt(sum(X.^2));Yd=[V; -G*ME*X/r^3];
% ------------------------------------------------------------------
function [ts,y0,options] = fi(tspan,Y0)
ts=tspan;y0 = Y0;
% 采用向域直接赋值法，设置options属性。以供与odeset使用方法对照。
options.RelTol=1e-5;options.AbsTol=1e-4;
options.OutputFcn='odephas2';	%在积分进程中，绘制相平面图。
options.OutputSel=[1 3];%解向量的第1、3分量分别为相平面图的横、纵坐标量。