我自定义函数[x,k,err]=Jacobi(A,b,n)后怎么才能输出k和err？

shuddkk · 发表于 2020-7-7 10:07

EDA365欢迎您登录！

您需要登录才可以下载或查看，没有帐号？注册

x

function [x,k,err]=Jacobi(A,b,n)
%Input A:系数矩阵
%Input b:右端向量
%Input n:系数矩阵的阶数
%Input epsilon:计算精度
%Input K:最大迭代次数
%Input x0:初始向量
%Output x:迭代近似解
epsilon=1e-6;
K=5000;
x0=rand(n,1);
k=1;
x=x0;
while k<K
for i=1:n
      s=b(i);
      for j=1:n
         s=s-A(i,j)*x(j);
      end
      s=(s+A(i,i)*x(i))/A(i,i);
      x(i)=s;
end
if norm(x-x0,inf)<epsilon
      err=norm(x-x0,inf);
      break;
end
k=k+1;
end
end

在命令窗口输入 Jacobi([10,5,0,0;5,10,-4,0;0,-4,8,-1;0,0,-1,5],[6;25;-11;-11],4)后
只得到了x的值，怎么输出k和err？