在解决圆形受限三体问题时，我遇到了一个僵硬的方程式。[考虑到由固定在二维空间中的两个引力源引起的引力效应，物体正在移动。]

方程是这些：

$x''=-\frac{GM_1 (x-x_1)}{\sqrt{(x-x_1)^2+y^2}^3}--\frac{GM_2 (x-x_2)}{\sqrt{(x-x_2)^2+y^2}^3}$

$y''=-\frac{GM_1 y}{\sqrt{(x-x_1)^2+y^2}^3}--\frac{GM_2 y}{\sqrt{(x-x_2)^2+y^2}^3}$

Euler Method 或 Runge Kutta 都不起作用，因为或附近的属性不好。导数变化太快。模拟无法正确解决。该物体太容易撞击引力源。$(x_1, 0)$$(x_1, 0)$

我怎样才能解决这个问题？

谢谢！