jacobi.m

用于求积分

function Jacobi(A,b,n,x0,E,N)
k=1
x=zeros(n,1);
 fprintf('\n**********start**************\n',k);
for k=1:N
    fprintf('%d\t',k);
    for i=1:n
        x(i)=b(i)/A(i,i);
        for j=1:n
            if j~=i
                x(i)=x(i)-A(i,j)*x0(j)/A(i,i);
            end
        end
        fprintf('%f\t',x(i));
    end
    if norm(x-x0)<E
            break;
    else
        k=k+1;
        x0=x;
    end
    fprintf('\n');
end
 fprintf('\n***********end*************\n',k);
x0
x