第二，启动MATLAB，新建脚本(Ctrl+N)，输入如下代码：

close all; clear all; clc

figure('Position',[50,50,600,500],'Name','double integration')

f=@(x,y) x*sin(y)-cos(x)+y-3;

Q1=dblquad(f,0,2*pi,-pi,pi,1.0e-3)

x=linspace(0,2*pi);y=linspace(-pi,pi);

ff=x.*sin(y)-cos(x)+y-3;

plot3(x,y,ff,'r','LineWidth',3)

grid on; view(35,35)

xlabel('x','FontSize',15); ylabel('y','FontSize',15)

其中Q1=dblquad(f,0,2*pi,-pi,pi,1.0e-3)采用默认方法quad计算二重积分，绝对计算精度设为1.0e-3。

plot3(x,y,ff,'r','LineWidth',3)是绘制被积分函数ff=x.*sin(y)-cos(x)+y-3的图像。