Если у вас небольшое количество тел / звезд, вы можете просто применять силы тяжести к своим телам каждый раз перед вызовом функции Step ().
Но я думаю, что это станет очень медленным, если система станет большой. Если ваши звезды и планеты не движутся, вы можете предварительно рассчитать общую силу тяжести в каждой точке 2D-пространства. Сохраните эти предварительно рассчитанные значения силы тяжести в некотором двумерном массиве (с некоторым шагом), а затем примените силы из этой справочной таблицы