matlab-solve函数CAT arguments dimensions not consistent.

这是我写的matlab代码，目的是用solve求解方程组

syms x0 y0
x1 = x1(1,1);y1 = y1(1,1);r1=r1(1,1);
x2 = x2(1,1);y3 = y2(1,1);r2=r2(1,1);
x3 = x3(1,1);y3 = y3(1,1);r3=r3(1,1);
f1=(x1-x0).^2+(y1-y0).^2-r1.^2;
f2=(x2-x0).^2+(y2-y0).^2-r2.^2;
f3=(x3-x0).^2+(y3-y0).^2-r3.^2;
[x0,y0] = solve([f1,f2,f3],[x0,y0])

运行出现以下问题：

错误使用 sym/cat>checkDimensions
CAT arguments dimensions not consistent.

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

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

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

请问我是哪里出问题了吗？
x1-3，y1-3，r1-3已确定值了，这是在matlab中的显示

x1 =
 
50.000000000000005176533632341287
 26.60444431189780282997260625709
 
 
y1 =
 
18.198511713310109693187965581336
46.080249255343817272810869349024
 
 
r1 =
 
53.208888623795609040627663555646
53.208888623795609040627663555646
 
 
x2 =
 
 26.60444431189777994361986611157
49.999999999999999641176942082799
 
 
y2 =
 
 -46.08024925534383048625278464076
-18.198511713310124901455480693607
 
 
r2 =
 
53.208888623795609040627663555646
53.208888623795609040627663555646
 
 
x3 =
 
              153.20888862379557733604333390191
0.000000000000010255259712890838942362539450747
 
 
y3 =
 
-0.000000000000031146691875878353378283308081177
0.0000000000000027249824454743459434383373780523
 
 
r3 =
 
99.999999999999990186904956574664
99.999999999999990186904956574664
 

报错只有这几行吗？没有get到你的程序中的报错行数，都是内部函数的报错，完整程序也只有这些吗？符号变量声明都不完整额，

clc,clear,close all;
x1 = [50.000000000000005176533632341287,26.60444431189780282997260625709];

y1 = [18.198511713310109693187965581336,46.080249255343817272810869349024];
 
 
r1 = [53.208888623795609040627663555646,53.208888623795609040627663555646];
 
 
x2 = [26.60444431189777994361986611157,49.999999999999999641176942082799];
 
 
y2 =[-46.08024925534383048625278464076,-18.198511713310124901455480693607];
 
 
r2 =[53.208888623795609040627663555646,53.208888623795609040627663555646];
 
 
x3 =[153.20888862379557733604333390191,0.000000000000010255259712890838942362539450747];
 
 
y3 =[-0.000000000000031146691875878353378283308081177,0.0000000000000027249824454743459434383373780523];
 
 
r3 =[99.999999999999990186904956574664,99.999999999999990186904956574664];
 
 
%syms x1 y1 x2 y2 x3 y3 x0 y0 r1 r2 r3
syms x0 y0
x1 = x1(1,1);y1 = y1(1,1);r1=r1(1,1);
x2 = x2(1,1);y3 = y2(1,1);r2=r2(1,1);
x3 = x3(1,1);y3 = y3(1,1);r3=r3(1,1);
f1=(x1-x0).^2+(y1-y0).^2-r1.^2;
f2=(x2-x0).^2+(y2-y0).^2-r2.^2;
f3=(x3-x0).^2+(y3-y0).^2-r3.^2;
[x0,y0]=vpasolve([f1,f2,f3],[x0,y0])

并无报错，但显示方程无解