Gauss-Seidel迭代法的Matlab程序.doc

《Gauss-Seidel迭代法的Matlab程序.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、Gauss-Seidel迭代法的程序functionx=Gauss_Seidel(A,b)[m,n]=size(A);x0=zeros(n,1);x=zeros(n,1);k=0;fori=1:nx(i)=b(i);forj=1:i-1x(i)=x(i)-A(i,j)*x(j);endforj=i+1:nx(i)=x(i)-A(i,j)*x0(j);endx(i)=x(i)/A(i,i);endxwhilenorm(x-x0)>1e-10k=k+1;x0=x;fori=1:nx(i)=b(i);forj=1:i-1x(i)=x(i)-A(

2、i,j)*x(j);endforj=i+1:nx(i)=x(i)-A(i,j)*x0(j);endx(i)=x(i)/A(i,i);endxendk