高斯塞尔德迭代法matlab

高斯塞尔迭代法（Gauss-Seidel Iteration）是数值计算中一种迭代法，用于求解线性方程组。Matlab代码如下：

function x = gauss_seidel(A, b, x0, tol, maxiter)
% Input:
%   A: coefficient matrix
%   b: right-hand side vector
%   x0: initial guess
%   tol: tolerance
%   maxiter: maximum number of iterations
% Output:
%   x: solution vector

n = size(A, 1);
x = x0;
for k = 1:maxiter
    for i = 1:n
        x(i) = (b(i) - A(i, 1:i-1) * x(1:i-1) - A(i, i+1:n) * x0(i+1:n)) / A(i, i);
    end
    if norm(x - x0) < tol
        return
    end
    x0 = x;
end
end

其中，A为系数矩阵，b为右侧向量，x0为初始猜测，tol为容差，maxiter为最大迭代次数。函数输出解向量x。