Как имитировать орбиты

Question

Когда я запускаю свой код, земля, похоже, стремится к пределу, а не к орбите вокруг Солнца, мне не хватает уравнений или что-то не так с моим кодом?

Здесь ошибка я получаю:

Warning (from warnings module): 
    File "C:\Python32\lib\site-packages\visual\visual_all.py", line 52 
    return numpy(x) 
RuntimeWarning: invalid value encountered in sqrt 

Вот код, который я написал:

from visual import * 

def SUVAT(A,B): 
     global EarthFinalV 
     global Acceleration 
     EarthFinalV = sqrt((A) + 2*Acceleration*(B)) 

GravitationalConstant = 1 

Sun = sphere(pos=(0,0,0), radius=10, color=color.red, 
      make_trail=True) 

Earth = sphere(pos=(50,0,0), radius=5, color=color.yellow, 
       make_trail=True) 

Sun.mass = 50 
Earth.mass = 10 

EarthInitialV = vector(0,1000,0) 
EarthFinalV = vector(0,0,0) 

while True: 
    rate(1) 
    Distance = Earth.pos - Sun.pos 

    GravitationalEquation = (GravitationalConstant*Sun.mass*Earth.mass)/mag(Distance)**2 
    Acceleration = GravitationalEquation/Earth.mass 

    SUVAT(EarthInitialV,Distance) 

    Earth.pos = Earth.pos - EarthFinalV 

Answer 1

Я заметил ваше гравитация уравнение является неправильным. Вы должны разделить на квадрат радиуса, а не в два раза больше радиуса. Вместо mag(Distance)*2 попробуйте Distance * Distance, или Distance**2