帮忙看看程序的问题。三元非线性方程组求解 程序一直报错求帮助

有没有什么求解三元非线性方程组的方法，求
syms x y z vx1 vx2 vy1 vy2 vz1 vz2 x1 x2 y1 y2 z1 z2 T F
lamda=6.372*10^9;
c=299792458;
R=6377830;
eq1=(vx1*(x-x1)+vy1*(y-y1)+vz1*(z-z1))/sqrt((x-x1)^2+(y-y1)^2+(z-z1)^2)-(vx2*(x-x2)+vy2*(y-y2)+vz2*(z-z2))/sqrt((x-x2)^2+(y-y2)^2+(z-z2)^2)-F*lamda;
eq2=sqrt((x-x1)^2+(y-y1)^2+(z-z1)^2)-sqrt((x-x2)^2+(y-y2)^2+(z-z2)^2)-c*T;
eq3=x^2+y^2+z^2-R^2;
[x,y,z]=solve(eq1,eq2,eq3,x,y,z)