gauss.m

为配套王能超版本数值分析简明教程编写,方便学生教师演示!

function gauss(A,b,e)
[i,j]=size(A);
[o,q]=size(b);
I=eye(i);
[l,u,p]=lu1(A);
if i==j & i==o;
if diag(p)~=zeros(i,1);
    D=inv(p-l);
    B=D*u;
    a=abs(eig(B));
    a=a';   
    n=0;
    if a<ones(1,i);
        m=1;
        x1=zeros(i,1);
        x2=B*x1+D*b;
       while max(abs(x1-x2))>e;
            x1=x2;
            x2=B*x1+D*b;
            fprintf('%d',m);
            m=m+1;
            for k=1:i;
                fprintf('%13.8f',x1(k));
            end
            fprintf('\n');
        end
    else
        fprintf('此多项式方程组不能进行雅可比迭代');
    end
else
    fprintf('此多项式方程组的矩阵的主对角线上有零元,不能进行雅可比迭代');
end
else
    fprintf('此多项式方程组不能进行雅可比迭代');
end