高斯赛德尔迭代法

 电脑爱好者网8SD7j&[4U$W5^

function[x,k,index]=Gau_Seid(A,b,ep,it_max)电脑爱好者网"m2`_&w1L.fv

if nargin<4 it_max=100;end

if nargin<3 ep=1e-5;end电脑爱好者网![ ]!Df!pff

n=length(A);k=0;

x=zeros(n,1);y=zeros(n,1);index=1;

  while 1

            y=x;电脑爱好者网.@;pr^C,^ u

            for i=1:n

            z=b(i);电脑爱好者网)Xe%N-sVF

            for j=1:n

                if j~=i

                    z=z-A(i,j)*x(j);电脑爱好者网b#y8G ?RnmC

                end

            end

            if abs(A(i,j))<1e-10|k==it_max电脑爱好者网/C,|/`-M\JM8bs

                index=0;return;

            end

            z=z/A(i,j);x(i)=z;

        end电脑爱好者网U,MV2n;U"j

        if norm(y-x,inf)<ep

            break;电脑爱好者网0W$I1~
d:U*M1V

        end电脑爱好者网3s?U3^0qDk

        k=k+1;电脑爱好者网wDT:JU\
M2V

end

>> A=[20 2 3;1 8 1;2 -3 15];b=[24 12 30]';ep=1e-5;

>> [x,k,index]=Gau_Seid(A,b,ep)

 电脑爱好者网ZiFat\

x =

   -5.0769电脑爱好者网 jYfU9s)b'])_-X

   12.0000

    5.0769电脑爱好者网1F,D#N^b#Lp

 电脑爱好者网FV;i7za4WTY

k =

 电脑爱好者网3D fD*VFT~4r

    34电脑爱好者网C)E2CDV\AA

index =

 电脑爱好者网d W6{1D#x"Io

     1电脑爱好者网zb)C%v @"Vp7i:MM