Для планеты массой м , на расстоянии r от корабля, корабль испытает ускорение:
a = k m / r ^ 2,
где k - некоторая константа, которая зависит от используемых вами единиц измерения.Ускорение будет направлено на планету.Может быть удобно разбить ускорение на составляющие по осям x и y (при условии, что вы работаете в двух измерениях).Если планета находится под углом тета в плоскости xy, относительно корабля,
a x = a cos ( theta )
a y = a sin ( theta )
Для нескольких планет вы можете просто добавить ускорения по компонентам.
Если корабль имеет начальную скорость v x во время t , тогда скорость на следующем шаге времени t + delta_t будет:
v x + a x * delta_t
Если этот корабль находится в начальной позиции p x в момент времени t, тогда позиция в t + delta_t будет:
p x + v x delta_t + a x delta_t ^ 2/2
См .: Уравнения движения