我必须模拟由于外力（无阻尼）而沿横向移动的圆柱体。运动方程为：$M\frac{d^2Y}{dt^2}+kY = F_Y$

我的困惑是关于而是在任意位置时如何获得位移。圆柱体应围绕 ( ) 点自由摆动。$kY$$(0,0)$$x_c=0.5, y_c=0.5$$0.5,0.5$

我想用欧拉法求解这个方程，在这种情况下我的版本是：$Y$

Y = Y0 + dt * V0;
V = V0 + dt * (-K * Y0 + FY) / M;

V0并且Y0是旧/前一个时间步的值。