这个错误是什么呀，像求解f点关于sita的函数，

sym x_f
sym y_f
l_ef=16;
x_e=-30;
y_e=-36;
l_bc=43;
l_cd=20;
y_d=0;x_d=0;

F=[x_f,y_f];
E=[x_e,y_e];
eny1=(y_f-y_e)*(y_f-y_c)/(x_f-x_e)/(x_f-x_c)+1;
eny2=(F-E)'*(F-E)-l_ef^2;
F_0=[-10,-15];
for sita=1:0.01:2*pi
    x_c=x_d+l_cd*cos(sita);
    y_c=y_d+l_cd*sin(sita);
    ENY=[eny1;eny2];
    n=1;
    while norm(ENY)>=0.01
        J(1,1)=diff(eny1,x_f);J(1,2)=diff(eny1,y_f);diff(eny2,x_f);diff(eny2,y_f);
        F=F_0-inv(J)*ENY;
        ENY=[eny1,eny2];
        n=n+1;
        F_0=F;
        if n>=1000
            n;
            break
        end
    end
end

ans =

x_f


ans =

y_f

错误使用 sym/cat>checkDimensions (第 68 行)
CAT arguments dimensions not consistent.

出错 sym/cat>catMany (第 33 行)
[resz, ranges] = checkDimensions(sz,dim);

出错 sym/cat (第 25 行)
    ySym = catMany(dim, args);

出错 sym/vertcat (第 19 行)
    ySym = cat(1,args{:});

>>