gauss.asv

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

function x1=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(1,i);
        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;
            end
        else
            fprintf('此多项式方程组不能进行高斯迭代');
        end
    else
        fprintf('此多项式方程组的矩阵的主对角线上的元素不能进行高斯迭代');
    end
else
    fprintf('此多项式方程组不能进行高斯迭代');
end