MATLAB求解非線性方程組

matlab中有專門的solve函數來解決方程組的
(a-x)^2+(b-y)^2=e^2
(C-x)^2+(D-y)^2=v^2
已知a,b,c,d,e,v 值
求解 X,Y 請問用 matlab 如何寫，就是求2個園的交點問題。
仿真程序為：
global a b c d e v;
>> a=1;b=0;c=-1;d=0;e=1.5;v=1.5;%設定你這幾個未知數的值
>> syms x y;%%%%%%x,y是變量
>> [x,y]=solve('x^2+y^2-2*a*x-2*b*y=e^2-a^2-b^2','x^2+y^2-2*c*x-2*d*y=v^2-c^2-d^2');%%%%把平方展開了
>> x=vpa(x,4);y=vpa(y,4);%%%%%%%%%%%取4位有效數字
%%%%%%%%%%
格式就是用solve(方程1，方程2，…求解變量1,變量2，…）；