分子动力学模拟的基础是牛顿运动方程的积分：

$$\frac{F}{m} = \frac{d^2r}{dt^2}$$

我知道有多种基于泰勒展开的方法，可以及时整合并计算新位置。然而，这样做总是会在展开式中留下一个错误项。

为了准确地解决位置，我可以整合双方：

$$r = \frac{Ft^2}{m} + c_1 t + c_2$$

如果$r(t=0) = 0$然后$c_2 = 0$.

如果可以确定$c_1$也（？）我们有一个精确的解决方案，所以我可以计算下一个位置$t+dt$与泰勒展开式不同，具有完美的准确性。

谁能告诉我我的理解错在哪里？谢谢你。