matlab里实现Gauss_Seidel迭代法的代码

以下是一个简单的 matlab 代码实现 Gauss-Seidel 迭代法的示例：

function [x, it] = gaussSeidel(A, b, x0, tol, maxit)
% Gauss-Seidel 迭代求解线性方程组 Ax = b
% A: 系数矩阵
% b: 常数向量
% x0: 初始解向量
% tol: 允许误差范围
% maxit: 最大迭代次数
% x: 解向量
% it: 实际迭代次数

n = length(b);
x = x0;
it = 0;
while it < maxit
it = it + 1;
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

提示：这里假设输入的线性方程组是无奇异的，并且 Gauss-Seidel 迭代法的收敛行为取决于系数矩阵 A 的条件数，因此在应用时需要格外注意收敛性的问题。