编写Gauss-Seidel迭代法的matlab程序，使其能够求解n阶线性方程组

可以回答这个问题。以下是一个示例的 Matlab 代码：

function [x, iter] = gauss_seidel(A, b, x0, tol, max_iter)
% Gauss-Seidel 迭代法求解线性方程组 Ax=b
% 输入参数：
% A: n x n 的系数矩阵
% b: n x 1 的常数向量
% x0: n x 1 的初值向量
% tol: 迭代精度
% max_iter: 最大迭代次数
% 输出参数：
% x: n x 1 的解向量
% iter: 实际迭代次数

n = length(b);
x = x0;
iter = 0;
while iter < max_iter
x_old = x;
for i = 1:n
x(i) = (b(i) - A(i,1:i-1)*x(1:i-1) - A(i,i+1:n)*x_old(i+1:n)) / A(i,i);
end
if norm(x - x_old) < tol
break;
end
iter = iter + 1;
end

if iter == max_iter
warning('达到最大迭代次数，可能未收敛！');
end

end

使用示例：

A = [4 -1 0 0; -1 4 -1 0; 0 -1 4 -1; 0 0 -1 3];
b = [15; 10; 10; 10];
x0 = [0; 0; 0; 0];
tol = 1e-6;
max_iter = 1000;
[x, iter] = gauss_seidel(A, b, x0, tol, max_iter);
disp(['解向量：']);
disp(x);
disp(['实际迭代次数：']);
disp(iter);