错误的NumPy数组计算结果

Question

我一直在努力弄清楚为什么它不起作用。必须说我是python新手。 尝试计算何时a ^ d + b ^ d = c ^ d; （费马定理）使用np.arrays 以获得最佳性能。 问题：对于d> 2，数组B中的结果错误。 请指出正确的方向吗？ 谢谢！

Kubuntu 17.10 Python 3.6.4 | Anaconda自定义（64位）| （默认值，2018年1月16日，18:10:19） 输入“版权”，“信用”或“许可证”以获取更多信息。 IPython 6.4.0-增强的交互式Python。

这是问题显示时的代码示例：

import numpy as np

d = np.float64(4.)
A = np.float64([[99001457116.0, 99536254152.0],
                [99002341176.0, 99914925909.0],
                [99029463918.0, 99750866948.0]])

B = np.float64((np.sum(A**2, axis=1)**(1/d)).reshape((A.shape[0], 1)))
B