Gravity:计算多个物体在重力作用下的位置。-matlab开发

% 调用序列命令 % 计算几个物体在重力作用下的位置和速度。 %基本策略是时间积分的Newmark方法和非线性问题的Newton-Raphson方法 重力2 % 给出起始体、发射时间和发射方向， % 计算撞击目标的速度。 %基本的数值策略是求根的割线法% 发现% 解决方案可能不存在[iniv,iiter,error] = findv(3,7,0.1,0.5,600,[0;-1],1e-5,1e-7,30); % 给定起始体和发射时间，计算初始方向% 发射速度最小（发射速度 = 逃逸速度） % 基本的数值策略是求根的割线法% 解决方案可能不存在[alpha1, iniv1, error, kiter] = secant1(3,7,0.1,300,0.5,1e-5,1e-7,30,-pi/2);