用Mathematica（或者别的软件）解方程组已知常量M=365;Ф=1.3Degree;Q=0.3Degree;a=60.4312Degree;Q2=3Degree;W=16.27715;t1=0.448567;t2=0.641098;Bsx11=5.088;Bsy11=5.78498;Bsz11=-17.1598;Bsx22=4.34157;Bsy22=3.6327;Bsz22=-21.4012;未知量为

　　用Mathematica（或者别的软件）解方程组

　　已知常量

　　M=365;

　　Ф=1.3Degree;

　　Q=0.3Degree;

　　a=60.4312Degree;

　　Q2=3Degree;

　　W=16.27715;

　　t1=0.448567;

　　t2=0.641098;

　　Bsx11=5.088;

　　Bsy11=5.78498;

　　Bsz11=-17.1598;

　　Bsx22=4.34157;

　　Bsy22=3.6327;

　　Bsz22=-21.4012;未知量为R,B,r,g,下列等式中Pi为圆周率π,未知量之间的关系见后,需要求解四个未知量R,B,r,g,其中R在0到60之间,B在-π到π之间,r在-π/2到π/2之间,g在0到2π之间,求解的四个未知量结果均在实数范围内.关系式如下（Mathematica软件式子）：

　　x0=R*Sin[Ф]*Cos[Ф-B]/Sqrt[Sin[Ф]^2+Sin[Ф-B]^2*Tan[Q-r]^2];

　　y0=R*Sin[Ф]*Sin[Ф-B]/Sqrt[Sin[Ф]^2+Sin[Ф-B]^2*Tan[Q-r]^2];

　　z0=R*Sin[Ф-B]*Tan[Q-r]/Sqrt[Sin[Ф]^2+Sin[Ф-B]^2*Tan[Q-r]^2];

　　x1=x0*Cos[Ф]*Cos[Q]+y0*Sin[Ф]*Cos[Q]+z0*Sin[Q];

　　y1=x0*(-Sin[Ф]*Cos[W*t1+g]-Cos[Ф]*Sin[Q]*Sin[W*t1+g])+y0*(Cos[Ф]*Cos[W*t1+g]-Sin[Ф]*Sin[Q]*Sin[W*t1+g])+z0*Cos[Q]*Sin[W*t1+g];

　　z1=x0*(-Cos[Ф]*Sin[Q]*Cos[W*t1+g]+Sin[Ф]*Sin[W*t1+g])+y0*(-Sin[Ф]*Sin[Q]*Cos[W*t1+g]-Cos[Ф]*Sin[W*t1+g])+z0*Cos[Q]*Cos[W*t1+g];

　　x2=x0*Cos[Ф]*Cos[Q]+y0*Sin[Ф]*Cos[Q]+z0*Sin[Q];

　　y2=x0*(-Sin[Ф]*Cos[W*t2+g]-Cos[Ф]*Sin[Q]*Sin[W*t2+g])+y0*(Cos[Ф]*Cos[W*t2+g]-Sin[Ф]*Sin[Q]*Sin[W*t2+g])+z0*Cos[Q]*Sin[W*t2+g];

　　z2=x0*(-Cos[Ф]*Sin[Q]*Cos[W*t2+g]+Sin[Ф]*Sin[W*t2+g])+y0*(-Sin[Ф]*Sin[Q]*Cos[W*t2+g]-Cos[Ф]*Sin[W*t2+g])+z0*Cos[Q]*Cos[W*t2+g];

　　Bx31=(4*Pi*10^-7*M*3*x1*z1/(4*Pi*(Sqrt[x1^2+y1^2+z1^3])

　　^5))*10^9;

　　By31=(4*Pi*10^-7*M*3*y1*z1/(4*Pi*(Sqrt[x1^2+y1^2+z1^3])

　　^5))*10^9;

　　Bz31=(4*Pi*10^-7*M*(2*z1^2-x1^2-y1^2)/(4*Pi*(Sqrt[x1^2+y1^2+z1^3])

　　^5))*10^9;

　　Bx32=(4*Pi*10^-7*M*3*x2*z2/(4*Pi*(Sqrt[x2^2+y2^2+z2^3])

　　^5))*10^9;

　　By32=(4*Pi*10^-7*M*3*y2*z2/(4*Pi*(Sqrt[x2^2+y2^2+z2^3])

　　^5))*10^9;

　　Bz32=(4*Pi*10^-7*M*(2*z2^2-x2^2-y2^2)/(4*Pi*(Sqrt[x2^2+y2^2+z2^3])

　　^5))*10^9;

　　Bx1=Bx31*Cos[Ф]*Cos[Q]+By31*(-Cos[Ф]*Sin[Q]*Sin[W*t1+g]-Sin[Ф]*Cos[W*t1+g])+Bz31*(-Cos[Ф]*Sin[Q]*C