用matlab编写高斯-赛德尔迭代法的代码

以下是一个简单的 Matlab 高斯-赛德尔迭代法的实现：

function [x, iter] = gauss_seidel(A, b, x0, tol, maxiter)
% 高斯-赛德尔迭代法求解线性方程组 Ax = b
% A：系数矩阵，b：常数矩阵，x0：迭代初值，tol：容差，maxiter：最大迭代次数
% x：方程组的解，iter：实际迭代次数

n = length(b);
x = x0;
iter = 0;

while iter < maxiter
    x_old = x;
    for i = 1:n
        sigma = 0;
        for j = 1:n
            if j ~= i
                sigma = sigma + A(i, j) * x(j);
            end
        end
        x(i) = (b(i) - sigma) / A(i, i);
    end
    iter = iter + 1;
    if norm(x - x_old) < tol
        break;
    end
end

end

其中，A、b、x0、tol 和 maxiter 分别代表系数矩阵、常数矩阵、迭代初值、容差和最大迭代次数。函数返回方程组的解 x 和实际迭代次数 iter。