pythonscipy.optimize.newton说它不收敛，但它看起来会

Question

我使用python和scipy编写一些代码，将笛卡尔坐标收敛到开普勒元素，反之亦然。在

为了将笛卡尔坐标转换为开普勒，我使用以下表达式：

E = scopt.newton(self.f, self.M, self.df, args=(), tol=10^(-10), maxiter=10000)

与

^{pr2}$

运行整个代码时，我得到一个错误：

RuntimeError: Failed to converge after 10000 iterations, value is 5.25182613825

如果我运行它的迭代次数更少（50次），我会得到：

RuntimeError: Failed to converge after 50 iterations, value is 5.25182613825

比较这两个值，它明显收敛。即使我将容差减少到10^（-2），我仍然会得到相同的运行时错误。在

有人知道为什么会发生这种错误吗？在

Answer 1

Python中的求幂运算符是**。使用tol=10**(-10)或{}。在

^是位异或。在