ykb.m

来自「雅可比迭代法解方程 高斯法解方程 查值法曲线拟和 2-9次方曲线拟和」· M 代码 · 共 28 行

%a=input('input array A:');
%b=input('input array b:');
a=[10 -1 -2;-1 10 -2;-1 -1 5];
b=[72 83 42];
n=length(b);
x=zeros(1,n);
N=20;   e=1e-5;    k=0;
while k<N
  y=x;
    for i=1:1:n   
      t=0;
      for j=1:1:n           
          t=t+a(i,j)*x(j);         
      end 
    x(i)=(b(i)-t+a(i,i)*x(i))/a(i,i);
    end
    if max(abs(x-y))<e   
      disp('done')                    %达到精度 输出done并显示结果和计算次数
      break
    else             
      k=k+1;                           %否则k加一返回继续
    end
end
disp(k+1)
disp(x)
disp(y)    %当计算完k次运算后仍未带到正确结果 则不执行y=x 保留y以观察