MATLAB行星运动仿真程序.docxVIP

clear all;close all;G=6.67*10^-11;% The universal gravitational constantm = [1.989e30,3.5844e23,4.89868e24,5.974e24,6.5714e23,1.89854e27,5.68725e26,8.72204e25,1.02753e26];% An array of massesn = length(m); Number_of_planets = n% The number of massesx_p = [0,58340100000,1.07705e11,1.4959e11,2.27377e11,7.77868e11,1.42709e12,2.87512e12,4.49668e12];% An array of x positionsy_p= [0,0,0,0,0,0,0,0,0];% An array of y positionsx_v = [0,0,0,0,0,0,0,0,0];% An array of x velocitiesy_v = [0,47856.46,34961.72,29780,23883.56,12924.52,9618.94,6789.84,5419.96];% An array of y velocitiesT