首先是对一元函数求积分，使用Scipy下的integrate函数： from scipy import integrate def g(x): return (1-x**2)**0.5 #用integrate下的quad函数可以同时求出积分结果和误差 res,err=integrate.quad(g,-1,1) #-1和1表示积分上下限，如果是正无穷用np.inf print(res,err) 得到的结果如下，前者是积分结果，后者是误差 1.5707963267948986 1.0002

问题1：求解如下二重积分 I=∫y=01/2∫x=01−2yxy dxdy. I = \int_{y=0}^{1/2}\int_{x=0}^{1-2y} xy~dxdy . I=∫y=01/2​∫x=01−2y​xy dxdy. 程序1，如下： from scipy.integrate import dblquad area = dblquad(lambda x, y: x*y, 0, 0.5, lambda x: 0, lambda x: 1-2*x) print