Рассчитайте неявные интегральные функции контуров с питоном

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.

Что я делаю неправильно? Есть ли другой способ решить эту проблему?

Я должен добавить, что я работаю с meshgrid:

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

    X,Y = meshgrid(x, y) 

Спасибо заранее!

Answer 1

Чтобы избежать этой ошибки, z_func должны быть векторизация:

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