这是我编写的M文件，但是不知道哪里错了，还请帮助！
%jacobi.m
function s=jacobi (A, b, X0, eps)
if nargin==3
    eps=1.0e-6
elseif nargin<3
    error
    return
end
% 求对角阵D和
D=diag(diag(A));
D=inv(D);
L=tril(A, -1);
U=triu(A, 1);
B=-D*(L+U);
f=D*b;
%进行迭代求解
while norm (s-X0) >=eps
    X0=s;
    s=B*X0+f;
end
return
当在 matlab中输入
a=[10 -1 2;2 5 -1;-1 2 5];
b=[5 10 15];
x0=[0 0 0];
eps=0.0001;
s=jacobi(a, b, x0, eps)
就会出现错误提示。