计算隐积分功能的轮廓与Python

Question

我想我缺少这个东西代码：

integrand = lambda t,x,y: (1/(Tiempo-t))*np.exp(-((x-U*(Tiempo-t))**2+y**2)/(4*a*(Tiempo-t))) 

def z_func(x,y,Rate,Conductivity): 
    integral, err = integrate.quad(integrand,0,Tiempo,args=(x,y,)) 
    return ((Rate/(2*math.pi* Conductivity))*integral) 

Z = z_func(X, Y, Ql, k) 

cs = plt.contour(X, Y, Z,[IncT]) 

我有一个隐函数与积分，像f(x,y,t)=A*Integral，其中A是恒定的。它集成了t。我需要为特定的t值计算轮廓。但是我收到了一些错误，例如"Supplied function does not return a valid float"，这是评估z_func时的实际错误。

我在做什么错？有其他解决方法吗？提前

 x = arange(-1.0,10.0,0.1) 
    y = arange(-1.0,10.0,0.1) 

    X,Y = meshgrid(x, y) 

感谢：

我要补充，我用meshgrid工作！

Answer 1

为了避免这种错误，z_func必须向量化：

vz_func = np.vectorize(z_func) 
Z = vz_func(X, Y, Ql, k)