轨道计算错误，位置不变

如何解决轨道计算错误，位置不变

我制作了一个程序来计算行星的轨道，但是当我看火星的x位置时，它并不依赖于它的-270192884240750.16单位长度。

我认为问题不在pr部分，问题在PE.py中。我以为位置r可能正在变化，但变化在0.000 ..部分。如果您查看火星的初始rx，然后在rx上方，您会发现它已经发生了非常大的变化。有关详细信息，您可以在这里了解我在做什么，并且可以帮助我:) Feynman Lectures On Physics： 您可以查看此链接的最后一页

事实上，我成功地模拟了地球和太阳，然后尝试推广该系统，您可以从here来访问它

这是我的代码： 我的星球对象在这里：星球代码

"""
import decimal
decimal.getcontext().prec = 100;
"""
class Planet:
    def __init__(self,name,m,r,v,F=[0,0],a=[0,0]):
        self.name = name
        self.r=r
        self.v=v
        self.F=F
        self.a=a
        self.m = m
        """
        for i in range(3):
            self.r[i] = decimal.Decimal(r[i])
            self.v[i] = decimal.Decimal(v[i])
            self.F[i] = decimal.Decimal(F[i])
            self.a[i] = decimal.Decimal(a[i])
            """
        self.pr = list()
        self.pv = list()
        self.pa = list()
        self.pr.append(r)
        self.pv.append(v)
        self.pa.append(a)
    def addPast(self,inp):
        if inp == 0:
            self.pr.append(self.r)
        elif inp == 1:
            self.pv.append(self.v)
        elif inp == 2:
            self.pa.append(self.a)

对于线：行星代码

from Planet import Planet as planet
#Mars  = Planet(m=6.471*10**23,r=(3396.2))
Mars = planet(name = "Mars",m = 6.4171*10**23,r = [1.093377665641286,-0.8345928928802313,-0.0451908725236022],v = [0.009032877133582754,0.01230339190527862,0.00003633981117282798])
Sun  = planet(name = "Sun",m= 1988500*10*24,r = [-0.005385370704196572,0.0068615262807763160,0.00006773925856402374],v = [-0.000007926400060373923,-0.000004128140422207171,0.000000229155339133599807])

planetList = [Mars,Sun]

主要代码：PE.py（物理引擎）

"""
import matplotlib.pyplot as plt
import decimal
"""
from Planets import planetList  as liste
from Graph import track

dt=0.1
G=6.674*(10**-11)
def calculateOrbit():
    for i in liste:
        for k in range(3):
            i.r[k] += i.v[k]*dt
            i.v[k] += i.a[k]*dt
            i.a[k] = 0
            i.F[k] = 0
            for j in liste:
                if(i!=j):
                    r = ((i.r[0] - j.r[0])**2 + (i.r[1] - j.r[1])**2 + (i.r[2] - j.r[2])**2)**(1/2)
                    i.F[k]+=-1*G*i.m*j.m*(i.r[k]-j.r[k])/r**3
                    i.a[k]+=i.F[k]/i.m
        i.addPast(0)
        i.addPast(1)
        i.addPast(2)
for s in range(1000):
    calculateOrbit()
track(liste[0])