求高手，高手速进

15223053365

Parameters
Variable
Parameter x, y ,z, dx1, dx2 ,dy1, dy2 ,dz1 ,dz2 ;
    Constant A1=[76 89 110 125 130 136 142 163 172];
             A=A1.*3,14/180;
             B1=[21.33 35.82 60.73 79.24 85.46 92.91 100.32 125.58 135.93];
             B=B1.*3.14/180;
             C1=[69.25 54.63 29.81 11.99 6.95 5.99 11.58 36.04 46.34];
             C=C1.*3.14/180;
             l=35800+6367;
Function
               For(i=1:9,A,B)     ((l*cosA-x)*dx1+dy1*(l*sinA-y)+dz1*(-z))/((l*cosA-x)^2+(l*sinA-y)^2+z^2)^0.5=cosB;
               For(i=1:9,A,C)     ((l*cosA-x)*dx2+dy2*(l*sinA-y)+dz2*(-z))/((l*cosA-x)^2+(l*sinA-y)^2+z^2)^0.5=cosC;
                    dx1*dx2+dy1*dy2+dz1+dz2=0
用MATLAB做出来时空集，怎么用牛顿RAPHSON做啊，或者用泰勒公式做，或者其他高招