其中,f是被积函数,可以是匿名函数,句柄,内联函数等等。a,b是最外重积分的常数项。c,d可以使常数也可以是匿名函数,代表内层积分的上下限。
举个例子说明:被积函数f(x,y) = sqrt(10^4-x^2) 在x^2+y^2<=10^4区域的积分
大家可以比较下面两行代码的运算时间 在我的电脑上,上面两种方法速度相差500多倍。
可见,quad2d才是真正有效求解一般区域二重积分的函数。
可惜现在没见到一般区域三重积分的函数,求解一般区域三重积分可能还得将其外推到长方体,然后用triplequad。希望未来的版本看到求解一般区域三重积分的函数。
举个例子说明:被积函数f(x,y) = sqrt(10^4-x^2) 在x^2+y^2<=10^4区域的积分
大家可以比较下面两行代码的运算时间 在我的电脑上,上面两种方法速度相差500多倍。
可见,quad2d才是真正有效求解一般区域二重积分的函数。
可惜现在没见到一般区域三重积分的函数,求解一般区域三重积分可能还得将其外推到长方体,然后用triplequad。希望未来的版本看到求解一般区域三重积分的函数。