高斯赛德尔迭代法

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

if nargin<4 it_max=100;end电脑爱好者网:Bn*yBwNr-[0[e

if nargin<3 ep=1e-5;end电脑爱好者网PG x2x+H*d `;p

n=length(A);k=0;

x=zeros(n,1);y=zeros(n,1);index=1;电脑爱好者网N|0T IMN}zZ

  while 1

            y=x;

            for i=1:n

            z=b(i);

            for j=1:n

                if j~=i

                    z=z-A(i,j)*x(j);

                end

            end

            if abs(A(i,j))<1e-10|k==it_max

                index=0;return;电脑爱好者网$y_\x dcY

            end

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

        end电脑爱好者网'GD\^VC

        if norm(y-x,inf)<ep电脑爱好者网4?_$pPmqd

            break;

        end电脑爱好者网L*qz3n3P:Q

        k=k+1;电脑爱好者网l(qW5UU
An"s

end电脑爱好者网-DS7]1de3p_:z

>> 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)

x =

 电脑爱好者网H#yz;} ` \` i}m

   -5.0769

   12.0000

    5.0769

 电脑爱好者网#fSH
v5f"y

 电脑爱好者网R4`}:sK

k =电脑爱好者网'e1s7? T0k

 电脑爱好者网0v/Y9|#Ie F#k.U3^E

    34电脑爱好者网d8?Z(YX:FqPh

index =电脑爱好者网dhb8B;]

     1电脑爱好者网/^q^!IYvU x^

 电脑爱好者网gA-mv%v7?/qD'HD)\r