第一，使用MATLAB計算下圖中的二重積分。

第二，啟動MATLAB，新建指令碼（Ctrl+N），輸入如下程式碼：close all; clear all; clcfigure("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的影象。

第三，儲存和執行上述指令碼，在命令列視窗（Command Window）得到如下結果：Q1 = -118.4351也就是說，該二重數值積分的結果為-118.4351.

第四，同時得到被積分函式x.*sin(y)-cos(x)+y-3的影象。

第五，在Q1=dblquad(f,0,2*pi,-pi,pi,1.0e-3)後面再新增一行命令，Q2=dblquad(f,0,2*pi,-pi,pi,1.0e-3,"quadl")。也就是在絕對計算精度都是1.0e-3的情況下，使用quadl方法計算二重數值積分。

第六，儲存和執行上述指令碼，在命令列視窗（Command Window）得到採用quad方法（預設）和quadl方法計算的二重數值積分結果。Q1 = -118.4351Q2 = -118.4355也就是說，採用採用quad方法（預設）計算的二重數值積分結果為-118.4351，而採用quadl方法計算的二重數值積分結果為-118.4355，兩者在小數點第四位略有差別。

在MATLAB軟體中輸入二重積分的程式碼即可求二重積分，具體操作請參照以下步驟，演示軟體版本為MATLAB 2014版。

1、將要使用MATLAB計算下圖中的二重積分，首先在電腦上開啟MATLAB軟體。

2、新建指令碼（Ctrl+N），輸入圖中框住的程式碼內容。其中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的影象。

3、儲存和執行上述指令碼，在命令列視窗（Command Window）得到如下結果：Q1 =-118.4351。也就是說，該二重數值積分的結果為－118.4351。

4、同時得到被積分函式x.*sin(y)-cos(x)+y-3的影象。

5、也可採用quadl法計算二重積分，在Q1命令後面再新增一行命令Q2=dblquad(f,0,2*pi,-pi,pi,1.0e-3,"quadl")即可。