syms x y a1=0;a2=50;b1=0;b2=0;r1=30;r2=40;%给出圆参数 [x,y]=solve((x-a1)^2+(y-b1)^2-r1^2,(x-a2)^2+(y-b2)^2-r2^2) %得到交点的解析解