功能: 多维函数高斯积分
格式:
IntegralFunByGauss(f, x1, a1, b1, x2, a2, b2, ..., xn, an, bn)
IntegralFunByGauss(f, x1, a1, b1, x2, a2, b2, ..., xn, an, bn, m)

f: n维函数符号对象 
x1,x2,...,xn: 积分变量名称 
a1,a2,...,an: 数值,积分变量下限 
b1,b2,...,bn: 数值,积分变量上限 
m : 数值,分割最大深度,默认为10000 

说明: 
1. 本函数使用3点高斯自适应方法进行积分
2. 本函数支持多维函数积分 
3. 本函数,在进行高纬度积分时,速度会很慢
4. 
$$
\begin{cases}
\mathbf{IntegralFunByGauss}(f, x, a, b) &=\mathbf{IntegralFun1D}(f, a, b)
\\
\\
\mathbf{IntegralFunByGauss}(f, x, ax, bx, y, ay, by)&=\mathbf{IntegralFun2D}(f, x, ax, bx, y, ay, by)
\end{cases}
$$

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

