使用Python为三体问题创建越级算法的问题

Question

我正在尝试使用跳蛙算法为 3 体问题编写代码。我使用 Piet Hut 和 Jun Makino 的“Moving Stars Around”作为指南。指南中的代码是用 C 编写的，但在尝试使用 Python 作为开始之前，我试图遵循确切的工作流程。我已经定义了时间步长dt = 0.01和持续时间t_end = 10，而不是提示输入。在第 5.4 节中，结果应该是：|gravity> g++ -o leapfrog2 leapfrog2.C|gravity> leapfrog2 > leapfrog2_0.01_10.outPlease provide a value for the time step0.01and for the duration of the run10Initial total energy E_in = -0.866025Final total energy E_out = -0.866025absolute energy error: E_out - E_in = 2.72254e-10relative energy error: (E_out - E_in) / E_in = -3.14372e-10连同一个圆形图。但是，我的代码的结果有所不同：Initial total energy E_in =  -0.8660254037844386Final total energy E_out =  -0.39922101519288833absolute energy error: E_out - E_in =  0.46680438859155027relative energy error: (E_out - E_in)/E_in =  -0.5390192788244604当然，在我绘制我的结果之后，它们不会绕圈子。我想知道我在翻译代码时是否犯了错误。任何帮助，将不胜感激！