以下代码旨在通过将每个重要的身体的影响结合到其他人身上来绘制太阳系.当然它应该导致预期的轨道.在最终的嵌入式函数中,GravityDE它无法从中读取值PlanetVec,因此,每次都无法生成正确的新结果.我们得到了错误
??? Undefined function 'GravityDE' for input arguments of type double.
任何有关如何解决这个问题的建议都是最受欢迎的!
function Gravity1()
clear;
format long eng;
load('solar_system_data.mat');
StartTime = 0;
TimeStep = 24 * 3600 * 10;
EndTime = 24 * 3600 * 100;
TVec = StartTime:TimeStep:EndTime;
TimeStepMin = StartTime:2:TimeStep;
%Column Vectors for initial conditions
SunVec = [xposition(1), yposition(1), vx(1), vy(1),mass(1),1];
MercuryVec = [xposition(2), yposition(2), vx(2), vy(2),mass(2),2];
VenusVec = [xposition(3), yposition(3), vx(3), vy(3),mass(3),3];
EarthVec = [xposition(4), yposition(4),vx(4), vy(4),mass(4),4];
MoonVec = [xposition(10), yposition(10), vx(10), vy(10),mass(10),10]; …