jacobi2.m

来自「Jacobi求解线性方程组、possison方程」· M 代码 · 共 26 行

clear all
clc

format long e
A=[8,-3 2;4 11 -1;6 3 12];
b=[20;33;36];
[n,m]=size(A);
x=zeros(n,1);
y=zeros(n,1);

diedai=0;
jingdu=1;
while jingdu>1e-6
    diedai=diedai+1;
    for j=1:n
        for k=1:n
            if k~=j
                b(j)=b(j)-A(j,k)*y(k);
            end
        end
        x(j)=b(j)/A(j,j);
    end
    jingdu=norm(x-y,inf);
    y=x;
end
x