四元方程求解，vpasolve函数。

大佬们帮忙看看这个四元方程为什么解不出来：（图为数学表达式）

eq1=sym('a/(pi/8*(a-b)^3*b+pi/8*(a-2*b-c)^3*c+pi/8*(a-2*b-2*c-0.02)^3*0.02))-645454.5455/(285551.776+50948.722*(1-a/5.8)+74796.986*(1-(a/5.8)^2))==0');

eq2=sym('a/100+0.00635-b==0');

eq3=sym('(a-2*b-((a-2*b)^4-64*d/pi)^0.25)/2-c==0');

eq4=sym('6.087*(4.772-pi/8*(a-b)^3*b-pi/8*(a-2*b-2*c-0.02)^3*0.02)-d==0');

S=vpasolve(eq1,eq2,eq3,eq4);

a=S.a

b=S.b

c=S.c

d=S.d