功能: 二维定积分

格式: IntegralFun2D(f, x, xa, xb, y, ya, yb)

f : 一个二维符号函数
x : 积分变量名称
xa: 数值, x的积分下限
xb: 数值, x的积分上限
y : 积分变量名称 
ya: 数值或者以x为变量的符号函数, 表示y的积分下限
yb: 数值或者以x为变量的符号函数, 表示y的积分上限

说明:
1. 本函数主要对函数f进行二重积分,且函数里的变量名称必须与x,y一致
2. xa,xb必须为数值
3. y的积分限可以为x的函数,也可以为数值,本函数原理如下
$$
\int_{xa}^{xb}\int_{ya}^{yb}f(x,y)dxdy
$$
4. 内部采用3点高斯积分

例子1:
这里要求对如下函数积分
$$
\int_0^\pi dx \int_0^1 x\cos(xy)dy
$$
f = @x * cos(x*y);
a = IntegralFun2D(f,x,0,pi,y,0,1)//执行得到如下结果,和理论解一致
a =
[ 2.00000000000000 ]

例子2:
这里要求对如下函数积分
$$
\int_{-1}^1dt\int_t^1 y\sqrt{1+t^2-y^2}dy
$$
ya =@t;
yb = 1;
xa = -1;
xb = 1;
f  = @ y * sqrt(1 + t^2 - y^2);
a = IntegralFun2D(f, t, xa, xb, y, ya, yb)//执行得到如下结果,和理论解一致
a =
[ 0.50000000000000 ]