j1 = 0.9
j3 = 0.9
j4 = 0.1
jz = 0.5
d1 = 0.1
d2 = 0.1
z = 4
r1 = 1/2 (Cos[x] + Cos[y])
r2 = Cos[x]*Cos[y]
a = 2*(j1 - j2 + j2*r2 - j1*r1)*z*sa + d1*(2*sa - 1) + jz*sb
b = 2*(j3 - j4 + j4*r2 - j3*r1)*z*sb + d2*(2*sb - 1) + jz*sa
c = jz*Sqrt[sa*sb]
ecl = (j2 - j1)*n*z*sa^2 - d1*n*sa^2 + (j4 - j3)*n*z*sb^2 - d2*n*sb^2 - jz*n*sa*sb
(*A+B>0*)
w1 = ((a - b)*((a + b)^2 - 4*c^2) + ((a + b)^2 + 4*c^2)*
Sqrt[(a + b)^2 - 4*c^2])/(2*((a + b)^2 - 4*c^2))
w2 = ((-a + b)*((a + b)^2 - 4*c^2) + ((a + b)^2 + 4*c^2)*
Sqrt[(a + b)^2 - 4*c^2])/(2*((a + b)^2 - 4*c^2))
q1 = w1/(Exp[w1/t] - 1)
q2 = w2/(Exp[w2/t] - 1)
sa = 0.5
sb = 0.5
j2 = 0.3
(*A+B>0*)
D[(1/(4*Pi^2))*NIntegrate[q1 + q2, {y, -Pi, Pi}, {x, -Pi, Pi}], t]
求助：1、最后一句求导中如果用NIntegrate，会导致出现以下错误并返回0值。N
Integrate::inumr: The integrand ((0. +4. (0.8 +Times[<<3>>]+Times[<<2>>])-4. (0.6 +Times[<<3>>]+Times[<<2>>])) ..... Plus[<<3>>]+4. Plus[<<3>>])^2)) has evaluated to non-numerical values for all sampling points in the region with boundaries {{0,3.14159},{0,3.14159}}.
怀疑应该是变量t没赋值导致的，但是我得对t求导，所以不能赋值啊，是不是就无法使用NIntegrate？
2.最后我用Integrate，但是一直无法运算出来，有什么便捷的方法么？
（PS：我最后只需要一个数值，不需要式子，所以才尝试用NIntegrate）
纠结了我好久了。
