与scipy整合

Question

我想将x**(-2/3)*(x-1)**(1/3)从x**(-2/3)*(x-1)**(1/3)进行积分，我曾经更改过变量并使用Simpson规则来解决它，现在我想尝试一些更快的方法。

我的代码：

x1 = lambda x: x**(-2/3)*(x-1)**(1/3)
integrate.quad(x1,0,1)

错误：

Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
  File "/usr/lib/python3/dist-packages/scipy/integrate/quadpack.py", line 281, in quad
    retval = _quad(func,a,b,args,full_output,epsabs,epsrel,limit,points)
  File "/usr/lib/python3/dist-packages/scipy/integrate/quadpack.py", line 345, in _quad
    return _quadpack._qagse(func,a,b,args,full_output,epsabs,epsrel,limit)
quadpack.error: Supplied function does not return a valid float.

我该如何解决？

Answer 1

即使标题说与numpy集成，我想您的意思是scipy。您的问题是您的函数中有复数。 例如，对于[0，1 (x-1)**(1/3) ，部分(x-1)**(1/3)变得复数，但是scipy.integrate.quad不处理复数。 有关更多信息，请参见例如， 使用scipy.integrate.quad来对复数进行积分 。