我试图用RK45实现一个关于地球和太阳的两体问题,但一直得到一个被零除的结果,这是我不明白的。从加速度函数看,除法似乎是正常的,但我不知道怎么做,也不知道怎么解决它。下面是代码: from scipy import optimize
from numpy import linalg as LA
import matplotlib.pyplot as plt
from scipy.optimize import fsolve
import numpy as np
AU=1.5e11
a=AU
e=0.5
mss=2E30
ms = 2E30
me = 5.98E24
mv=4.867E24