def Cl(y):
Td = 420
tt = Td/y
fx = lambda x: x**4*np.exp(x)/(np.exp(x)-1)**2
[cl,err] = 9*N*kb*(y/Td)**3*quad(fx,0,tt)
return cl
здесь N, kb - глобально определенные константы.Здесь идея состоит в том, чтобы вычислить cl для каждого значения y, которое передается внутри Cl (y).Может кто-нибудь сказать, является ли это правильным способом использования 'quad' внутри функции?если да, не могли бы вы объяснить эту ошибку?