算法步骤
例题
程序代码
function x=GaussSeidel(A,b,x0,eps,N)
n=length(b);
if nargin<5
N=500;
end
if nargin<4
eps=1e-6;
end
if nargin<3
x0=zeros(n,1);
end
k=0;
x=zeros(n,1);
while k<N
x(1)=(b(1)-A(1
function x=GaussSeidel(A,b,x0,eps,N)
n=length(b);
if nargin<5
N=500;
end
if nargin<4
eps=1e-6;
end
if nargin<3
x0=zeros(n,1);
end
k=0;
x=zeros(n,1);
while k<N
x(1)=(b(1)-A(1