jacobi.asv

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

function jacobi(A,b,e)
[i,j]=size(A);
a=ones(i);
a=a';
if (abs(eig(A))<a);
    [l,u,p]=lu1(A);
    m=1;
    x1=zeros(i,1);
    I=eye(i);
    x2=(I-inv(p)*A)*x1+inv(p)*b;
    while max(abs(x1-x2))>e;
       x1=x2;
       x2=(I-inv(p)*A)*x1+inv(p)*b;
       fprintf('%d',m);
       m=m+1;
       for k=1:i;
            fprintf('%13.8f',x1(k));
       end
       fprintf('\n');
    end
else
    fprintf('此多项式方程组不能进行雅可比迭代');
end