改为有四个积分节点

将二重四点高斯求积公式改为四重四点高斯求积公式，可以得到以下 Matlab 代码实现：

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

% 定义积分区间和高斯积分系数
a = 0; b = 1;
c = 0; d = 1;
A = 1; B = 1;
x1 = -1/sqrt(3); w1 = 1;
x2 = 1/sqrt(3); w2 = 1;
x3 = -1/sqrt(3); w3 = 1;
x4 = 1/sqrt(3); w4 = 1;

% 计算积分
I = 0;
for i = 1:4
    for j = 1:4
        I = I + w1*w2*w3*w4*f((b+a)/2+x1*(b-a)/2,(d+c)/2+x2*(d-c)/2);
        x2 = -x2; % 对称性
    end
    x1 = -x1; % 对称性
    x2 = sqrt(3)/3; % 恢复
    x3 = -x3; % 对称性
    x4 = sqrt(3)/3; % 恢复
end
I = A*B*(b-a)*(d-c)/16*I;

在上述代码中，我们增加了两个高斯积分点，并相应地修改了计算过程。需要注意的是，四重高斯积分公式需要计算 $4 \times 4 = 16$ 个积分点，因此计算量比二重高斯积分公式要大得多。