一道c语言或matlab编程题

function varargout=Jacobiliu(varargin)
A=[5 2 1;2 8 -3;1 -3 -6];
b=[8;21;1];
x0=[0;0;0];
x=Jacobi(A,b,x0)
function x=Jacobi(A,b,x0)
n=100;%最大迭代次数
ee=0.0001;%精度
n1=length(b);
x1=zeros(n1,1);
for i=1:n
for j=1:n1
s=0;
for k=1:n1
if k~=j
s=s+A(j,k)*x0(k);
end
end
x1(j)=(b(j)-s)/A(j,j);
end
if norm(x1-x0) break
end
x0=x1;
end
x=x1;

结果
x =

1.0000
2.0000
-1.0000

---

5x1+2x2+x3=8
2x1+8x2-3x3=21
X1-3x2-6x3=1
要求源代码及运算结果

---