dblquadandquad2d

其中，f是被积函数，可以是匿名函数，句柄，内联函数等等。a,b是最外重积分的常数项。c,d可以使常数也可以是匿名函数，代表内层积分的上下限。
举个例子说明：被积函数f(x,y) = sqrt(10^4-x^2) 在x^2+y^2<=10^4区域的积分
大家可以比较下面两行代码的运算时间

1. tic,y1 = dblquad(@(x,y) sqrt(10^4-x.^2).*(x.^2+y.^2<=10^4),-100,100,-100,100 ),toc
   
2. tic,y2 = quad2d(@(x,y) sqrt(10^4-x.^2),-100,100,@(x)-sqrt(10^4-x.^2),@(x)sqrt(10^4-x.^2)),toc

复制代码

在我的电脑上，上面两种方法速度相差500多倍。

可见，quad2d才是真正有效求解一般区域二重积分的函数。

可惜现在没见到一般区域三重积分的函数，求解一般区域三重积分可能还得将其外推到长方体，然后用triplequad。希望未来的版本看到求解一般区域三重积分的函数。