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

题目

用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*t

答案

你能不能用标准格式编写了给我,我用lingo帮你求解；