使用MATLAB，利用函数或函数句柄求解二重积分 ，其中积分区间D是由圆周 及坐标轴所围成的第一象限内的闭区域。

我们可以使用MATLAB中的“integral2”函数来求解二重积分，具体步骤如下：

1. 定义被积函数f(x,y)，即在积分区间D中，f(x,y)的取值。

例如，若要求解二重积分∬D(x^2+y^2) dxdy，那么被积函数f(x,y)就是x^2+y^2。

2. 定义积分区间D，即圆周及坐标轴所围成的第一象限内的闭区域。

我们可以通过定义一个函数句柄来实现这一步骤，例如：

D = @(x,y) x.^2 + y.^2 <= 1 & x >= 0 & y >= 0;

这里，我们使用MATLAB中的“.^”操作符表示对向量或矩阵中的每个元素都进行平方操作。

3. 调用“integral2”函数进行积分计算。

例如，要计算上述被积函数在积分区间D中的二重积分，可以使用以下代码：

result = integral2(@(x,y) x.^2 + y.^2, 0, 1, 0, 1, 'Method', 'iterated');

其中，第一个参数是被积函数f(x,y)的函数句柄；第二个参数和第三个参数分别是积分区间D在x轴方向上的上下限；第四个参数和第五个参数分别是积分区间D在y轴方向上的上下限；最后一个参数指定使用迭代法进行积分计算。

完整代码如下：

D = @(x,y) x.^2 + y.^2 <= 1 & x >= 0 & y >= 0;
result = integral2(@(x,y) x.^2 + y.^2, 0, 1, 0, 1, 'Method', 'iterated');
disp(result);

运行结果为：

result = 0.7854

即所求二重积分的结果为0.7854。