用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*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]*Cos[W*t1+g]+Sin[Ф]*Sin[W*t1+g]);
By1=Bx31*Sin[Ф]*Cos[Q]+By31*(Cos[Ф]*Cos[W*t1+g]-Sin[Ф]*Sin[Q]*Sin[W*t1+g])+Bz31*(-Cos[Ф]*Sin[W*t1+g]-Sin[Ф]*Sin[Q]*Cos[W*t1+g]);
Bz1=Bx31*Sin[Q]+By31*Cos[Q]*Sin[W*t1+g]+Bz31*Cos[Q]*Cos[W*t1+g];
Bx2=Bx32*Cos[Ф]*Cos[Q]+By32*(-Cos[Ф]*Sin[Q]*Sin[W*t2+g]-Sin[Ф]*Cos[W*t2+g])+Bz32*(-Cos[Ф]*Sin[Q]*Cos[W*t2+g]+Sin[Ф]*Sin[W*t2+g]);
By2=Bx32*Sin[Ф]*Cos[Q]+By32*(Cos[Ф]*Cos[W*t2+g]-Sin[Ф]*Sin[Q]*Sin[W*t2+g])+Bz32*(-Cos[Ф]*Sin[W*t2+g]-Sin[Ф]*Sin[Q]*Cos[W*t2+g]);
Bz2=Bx32*Sin[Q]+By32*Cos[Q]*Sin[W*t2+g]+Bz32*Cos[Q]*Cos[W*t2+g];
Bsx11=Bx1*Sin[a]*Cos[Q2]-By1*Cos[a]*Cos[Q2]+Bz1*Sin[Q2];
Bsy11=Bx1*(Cos[a]*Cos[Q2]-Sin[a]*Sin[Q2]*Sin[Q2])+By1*(Sin[a]*Cos[Q2]+Cos[a]*Sin[Q2]*Sin[Q2])+Bz1*Sin[Q2]*Cos[Q2];
Bsz11=Bx1*(-Sin[a]*Sin[Q2]*Cos[Q2]-Cos[a]*Sin[Q2])+By1*(Cos[a]*Sin[Q2]*Cos[Q2]-Sin[a]*Sin[Q2])+Bz1*Cos[Q2]*Cos[Q2];
Bsx22=Bx2*Sin[a]*Cos[Q2]-By2*Cos[a]*Cos[Q2]+Bz2*Sin[Q2];
Bsy22=Bx2*(Cos[a]*Cos[Q2]-Sin[a]*Sin[Q2]*Sin[Q2])+By2*(Sin[a]*Cos[Q2]+Cos[a]*Sin[Q2]*Sin[Q2])+Bz2*Sin[Q2]*Cos[Q2];
Bsz22=Bx2*(-Sin[a]*Sin[Q2]*Cos[Q2]-Cos[a]*Sin[Q2])+By2*(Cos[a]*Sin[Q2]*Cos[Q2]-Sin[a]*Sin[Q2])+Bz2*Cos[Q2]*Cos[Q2];
麻烦达人了,谢谢
麻烦达人将编程过程和结果一起弄出来,对了,上面的式子中-7是指-7次方,谢谢

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