闭曲面积分的高斯公式中闭曲面,对常数进行二重积分

from scipy.special import roots_legendre def f(x): return x**4-2*x+1N = 3 #取3个样本点a,b = 0,2 #积分上下限 x,w = roots_legendre(N)xp = x*(b-a)/2+(b+a)/2wp = w*(b-a)/2s = 0for i in range(N): s+=wp[i]*f(xp[i])print(s)