f = @(x,y)y*sin(x)+x*cos(y);
% 編寫匿名函數，將句柄賦給f
S = dblquad(f,pi,2*pi,0,pi)
% 計算二重積分

% 符號計算能得到高精度的解，語法請查閱下章
syms x y;
% 定義符號變量
S1 = int(int(y*sin(x)+x*cos(y),x,pi,2*pi),y,0,pi)
% 計算解析解

S2 = vpa(S1)
% 得高精度解