jocabi_1.txt

来自「本代码是本人在做作业时所编」· 文本 代码 · 共 31 行

function x=jacobi_1(a,b)
tol=1e-6;
max=100;
[n,m]=size(a);
x0=[0;0;0];
xold=x0;
B=-a;
for i=1:n
B(i,i)=0;
B(i,:)=B(i,:)/a(i,i);
g(i,1)=b(i)/a(i,i);
end
i=1;
while i<max;
xnew=B*xold+g;
if norm(xnew-xold)<=tol
   x=xnew;
   disp('迭代收敛');
   return;
else
   xold=xnew;
  end
  disp([xnew']);
  i=i+1;
end
disp('不收敛');
disp('jiguo')
x=xnew;