Gauss-Seidel迭代法(8)

else

MaxInteration=10; n=length(A(:,1)); D=zeros(n,n); L=zeros(n,n); U=zeros(n,n); for i=1:n for j=1:n if(i==j)

D(i,j)=A(i,j); else

if(i>j)

L(i,j)=-A(i,j); else

U(i,j)=-A(i,j); end end end end

G=inv(D-L)*U; f=inv(D-L)*b; x1=zeros(n,1); xhis=x1;

for i=1:MaxInteration x2=G*x1+f; x1=x2;

xhis=[xhis,x1]; end

xhis=xhis'; y=x1;

if(n==3)

plot3(xhis(:,1),xhis(:,2),xhis(:,3),'or');

title('解的迭代情况'); end end