符号计算
syms x y;                    %定义两个符号变量
a=int(int(x^y,x,0,1),y,1,2)  %积分x,0,1   ,y,1,2
b=simple(a)                  %化简
c=vpa(b,4)                   %得到4位近似解，也可以任意N位解
数值计算
%%二重积分f= @(x,y)exp(sin(x))*ln(y)，y从5*x积分到x^2，x从10积分到20
1 (7.X后版本才有此函数quad2d)
   y1=quad2d(@(x,y) exp(sin(x)).*log(y),10,20,@(x)5*x,@(x)x.^2)

2  y2 =quadl(@(x) arrayfun(@(x) quadl(@(y)exp(sin(x)).*log(y),5*x,x.^2),x),10,20)

3  y3 = dblquad(@(x,y)exp(sin(x)).*log(y).*(y>=5*x & y<=x.^2),10,20,50,400)