x,y))#不定积分print(sympy.integrate(f,(x,0,1),(y,0,1)))#二重积分定积分 ...
没试过，不过可以推荐尝试下这个scipy.integrate.dblquad
φ:np.sin(θ)*np.cos(θ)# 要计算的表达式# 二重积分value,error=dblquad(f,# 函数表达式0,# ...