没其他条件吗？那例子太多了。
给你一个刚才帮别人调的Jacobi迭代求解线性方程组的程序吧。
函数（保存成Jacobi.m）：
epsilon=0.001;
n=length(b);
x1=zeros(n,1);
M=100;
for k=1:1:M
for i=1:1:n
s=sum(A(i,:)*x0(:))-A(i,i)*x0(i);
x1(i)=(b(i)-s)/A(i,i);
end
fprintf('%.4f ',x1);
fprintf('\n ');% 按格式输出数据
if norm(x1-x0,'inf')
end
x0=x1;
end
x=x1;
调用：
-1., 11., -1., 3.;
2., -1., 10., -1.;
0.0, 3., -1., 8.];
b = [6., 25., -11., 15.]';
x0 = zeros(4,1);
x=Jacobi(A,b,x0)
