Python的运算符的一类问题，Z = np.cos((X**3+Y**3)**(1/3))

Question

import numpy as np
x = np.linspace(-4,4,num=50)
y = np.linspace(-4,4,num=50)
X,Y = np.meshgrid(x,y)
Z = np.cos((X**3+Y**3)**(1/3))
程序提示报错：
D:\Anaconda3\lib\site-packages\ipykernel_launcher.py:5: RuntimeWarning: invalid value encountered in power
"""

但是当我将程序中的Z = np.cos((X**3+Y**3)**(1/3))改成Z = np.cos((X**2+Y**2)**(1/3))，程序就能正常运行，请问各位这是怎么回事啊？

Answer 1

可以将这段代码改下
x = np.linspace(-4,4,num=50,dtype=np.complex)
y = np.linspace(-4,4,num=50,dtype=np.complex)
不改dtype的话，在出现负数的根时会出现问题。