1.程序代码
funcon x=Gauss(A,b,x0,ep,N)
%用途:用高斯迭代法解线性方程组 Ax=b
%A 为系数矩阵,b 为右端向量,x0 为初始向量(默认零向量)
%ep 为精度(1e-6),N 为最大迭代次数(默认 500 次),x 返回近似解向量
n=length(b);
if nargin<5
N=500;
end
if nargin<4
ep=1e-6;
end
if nargin<3
x0=zeros(n,1);
end
x=zeros(n,1);
k=0;
while k<N
for i=1:n
if i==1
x(1)=(b(1)-A(1,2:n)*x0(2:n))/A(1,1);
elseif i==n