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，…）；