我有用牛顿迭代法求非线性方组的过程中遇到问题,希望高手们帮忙看一下啊,附上我自己的程序   
    function [r,n]=mulNewton3(F,x0,eps)
if nargin==1
    eps=1.0e-4;
end
J=jacobian(sym(F));
fa=subs(sym(F),x0);
fb=subs(sym(J),x0);
r=x0-fa/fb;
n=1;
tol=1;
while tol>eps
    x0=r;
    r=x0-fa/fb;
    tol=norm(r-x0);
    n=n+1;
    if(n>100000)
        disp('******');
        return;
    end
end

下现的部分就是出错的地方的,谢谢
J=jacobian(sym(F));
fa=subs(sym(F),x0);
fb=subs(sym(J),x0);