实现velocity verlet算法

Question

我使用了以下公式来实现velocity verlet算法(在这里用一种伪代码编写)：

time+=timeStep;
newPosition+=timeStep*(velocity+acceleration*timeStep*0.5);

newA = Force/mass, force is a function of newPosition

v+= 0.5*timeStep*(newA-acceleration);

这不管用！newA-acceleration非常小。我用错公式了吗？或者速度verlet算法对于小的timeSteps (我使用的是2秒)不能很好地工作？

(如果这个放错了地方，我很抱歉，我不知道还能把它放在哪里)

Answer 1

应该是newA +加速，而不是newA-加速。糟了！