数值分析 jacobi迭代法求解线性方程组 MATLAB程序实现

jacobi迭代法 参考数值分析第四版 颜庆津著 P36 

数值分析 jacobi迭代法MATLAB程序实现   运行输入为：

运行结果为：

以下是函数内容（保存为jacobi.m文件，在MATLAB中运行）：

%function [G,d,x,N]=jacobi(A,b)

%Jacobi迭代法求解线性方程组。迭代公式x(k+1)=Gx(k)+d,(k=0,1,…).
%返回迭代矩阵G，d和方程组的解x，以及迭代次数N


%2015.11.6  密密编写  (*^__^*) 嘻嘻……
function [G,d,x2,N]=jacobi(A,b)
n=size(A);
if n(1)~=n(2)
    error('矩阵A不是方阵');
end

%初始化
N=1;%迭代次数
%X=zeros(N,n)%存储每次迭代解出的x，迭代次数N
B=zeros(n)