用Matlab 计算二重积分，其中D是由直线及双曲线 所围成的闭区域。

假设要计算的二重积分为：

∬Df(x,y)dxdy

其中D为由直线及双曲线所围成的闭区域，f(x,y)为被积函数。

可以通过将D分成若干个小矩形，再利用双重嵌套循环对小矩形进行遍历和面积计算，从而计算出二重积分的近似值。

具体步骤如下：

1. 绘制出D的图像，并确定D的极限值。

2. 在MATLAB中定义被积函数f(x,y)。

3. 定义二维区域D。

4. 计算每个小矩形的面积和中心坐标。

5. 在双重嵌套循环中，计算每个小矩形的面积和中心坐标，并将其加到积分结果中。

6. 显示结果。

以下是一个示例代码，实现对由直线及双曲线所围成的闭区域的二重积分计算：

% 定义被积函数
f = @(x,y) x.^2 + y.^2;

% 定义二维区域D
x = linspace(-2,2,201);
y1 = sqrt(x.^2 - 1);
y2 = -sqrt(x.^2 - 1);
D = [x', y1'; x', y2'];

% 计算每个小矩形的面积和中心坐标
n = size(D,1);
delta = (x(2)-x(1));
S = delta*sqrt(1+(2*x./(x.^2-1)).^2);
x_center = D(:,1);
y_center = D(:,2);

% 计算二重积分的近似值
result = 0;
for i = 1:n
    result = result + f(x_center(i),y_center(i))*S(i);
end

% 显示结果
disp(result);

在上述代码中，我们首先定义了被积函数f(x,y)，然后通过linspace函数生成了x坐标的取值范围，并计算出了y坐标的取值范围。接着，我们通过meshgrid函数生成了D中每个小矩形的顶点坐标，并计算出了每个小矩形的面积和中心坐标。最后，在双重嵌套循环中，我们计算了每个小矩形的面积和中心坐标，并将其加到积分结果中，最终得到了二重积分的近似值。