高斯赛德尔迭代法

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

if nargin<4 it_max=100;end

if nargin<3 ep=1e-5;end

n=length(A);k=0;

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

  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;

            end

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

        end

        if norm(y-x,inf)<ep

            break;

        end

        k=k+1;

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)



x =



   -5.0769

   12.0000

    5.0769





k =



    34





index =



     1

..呵呵 我们这学期的计算方法讲过这个东东。。迭代的事还是交给计算机吧。

几头雾水…

本人姓:常...名:有理..

让我们和财富一起成长!!!