矢量化 N 体方程：二阶和一阶形式的完全矢量化 n 体运动方程。-matlab开发

NBODYVECT - 完全矢量化的 n 体方程nbodyVect(x0,dx0,mus,t) 使用辛二阶积分器或任何标准内置或用户提供的一阶积分器对 n 体运动方程进行积分。 x0 - 3n x 1 堆叠的初始位置向量： [r1(1);r1(2);r1(3);r2(1);r2(2)r2(3);...;rn(1);rn(2);rn(3)] dx0 - 3n x 1 堆叠的初始速度向量mus - 引力参数 (G*m_i)，其中 G 是引力常数，m_i 是第 i 个物体的质量。 t - 时间跨度：初始和最终时间的 2 元素数组，或要输出的所有时间步长的数组。 所有单位必须是互补的，即，如果位置以 AU 为单位，时间以天为单位，则 dx0 必须以 AU/天为单位，mus 必须以 AU^3/天 ^2 为单位（这些是 solarSystemData.mat 中的单位）。 该数据是从喷气推进实验