想要自己复现一下文献，但是下面这个方程组用mathematica一直运行不出来，想请教大家有什么技巧能够完成计算呢？需要求解出z11，z12，z21，z22四个量（其他均为已知量，用符号表示），或者求出（z21-z11）和（z22-z12）也行。是不是符号求解太困难而只能用把已知量数据代入数值求解呢？
eq1 = -r^2*z21 + a1*z11 + 3/4*a3*(z11^3 + z11*z12^2) +
5/8*a5*(z11^5 + 2*z11^3*z12^2 + z11*z12^4) + 2*xi1*r*z12
eq2 = -r^2*z22 + a1*z12 + 3/4*a3*(z12^3 + z11^2*z12) +
5/8*a5*(z12^5 + 2*z11^2*z12^3 + z11^4*z12) - 2*xi1*r*z11
eq3 = -r^2/mu*(z21 - z11) - a1*z11 - 3/4*a3*(z11^3 + z11*z12^2) -
5/8*a5*(z11^5 + 2*z11^3*z12^2 + z11*z12^4) - 2*xi1*r*z12 +
a*(z21 - z11) + 2*xi2*d/mu*(z22 - z12)
eq4 = -r^2/mu*(z22 - z12) - a1*z12 - 3/4*a3*(z12^3 + z11^2*z12) -
5/8*a5 (z12^5 + 2*z11^2*z12^3 + z11^4*z12) + 2*xi1*r*z11 +
a*(z22 - z12) + 2*xi2*d/mu*(z11 - z21) - X
eq1 = -r^2*z21 + a1*z11 + 3/4*a3*(z11^3 + z11*z12^2) +
5/8*a5*(z11^5 + 2*z11^3*z12^2 + z11*z12^4) + 2*xi1*r*z12
eq2 = -r^2*z22 + a1*z12 + 3/4*a3*(z12^3 + z11^2*z12) +
5/8*a5*(z12^5 + 2*z11^2*z12^3 + z11^4*z12) - 2*xi1*r*z11
eq3 = -r^2/mu*(z21 - z11) - a1*z11 - 3/4*a3*(z11^3 + z11*z12^2) -
5/8*a5*(z11^5 + 2*z11^3*z12^2 + z11*z12^4) - 2*xi1*r*z12 +
a*(z21 - z11) + 2*xi2*d/mu*(z22 - z12)
eq4 = -r^2/mu*(z22 - z12) - a1*z12 - 3/4*a3*(z12^3 + z11^2*z12) -
5/8*a5 (z12^5 + 2*z11^2*z12^3 + z11^4*z12) + 2*xi1*r*z11 +
a*(z22 - z12) + 2*xi2*d/mu*(z11 - z21) - X
