Уважаемый AAXEE в своих рассуждениях вы глубоко ошибаетесь. Потому что постановка задачи ошибочная по своей физической сути. Вы предлагаете использовать закон всемирного тяготения, но дело в том что использовать этот закон в первоначальном виде для расчета параметров орбиты нельзя, так как закон всемирного тяготения только констатирует что любые два тела во вселеной притягиваються с силой равной F=G*m1*m2/r^2 где r не радиус орбиты как вы утверждаете а расстояния между двумя телами. Кроме того вы в своей задачи предлагаете подбирать G , на самом деле G является константой которая равна 6.6742(10)·10-11 м3кг-1c-2, для расчетов параметров орбиты можно использовать законы Кеплера. И исходя из законов Кеплера можно утверждать что параметры орбиты не меняются.
если масса одного объекта во много раз больше массы другого, то радиус орбиты и будет расстоянием между телами.А вобще не суть важно: пусть это будут два одинаковых объекта и пусть они крутятся вокруг друг друга. Расстояние между ними будет увеличиваться.
Насчет G. Естественно, можно выбрать G, совпадающее с реальным, но это не обязательно, это всего лишь коэффициент, характеризующий расстояние взаимодейстия. Если я хочу, что бы два шарика массами по 1 кг крутились на расстоянии метра достаточно быстро, я вполне могу подобрать G таким, чтобы это осуществить.
Если речь идёт о моделировании движения, то при решении дифференциального уравнения второго порядка итерационно численными методами естественно будет погрешность накапливаться.
Это вовсе не означает, что вычисленное значение будет отклоняться от "фактического" бесконечно в одну сторону, просто доверие к результату будет уменьшаться.
Расстояние между телами стабильно увеличивается всегда. сначала медленно , потом все быстрее и быстрее.
Кроме того проще вычислить параметры орбиты и подставлять время в готовую формулу, а не итерировать.
А еще проще просто заставить объект летать функционально от времени.