多元高次方程的求解，公式推导？

图片中给出了公式，最后要推导出a的表达式，我推导出来一个关于a的四次方程，不会解了，求助大佬们。
解决了发红包。

syms a h f b E d u r x y w q z D sigma_x sigma_y Mx M_y M_xy M_yx shear_xy shear_yx sigma_max y1 y2 t_max s ;

% h=a/(2*f);

% q=h*r;
% w=(4*q*(a^4)*(1+cos(2*pi*x/a))*(1+cos(2*pi*y/b)))/((pi^4)*D*(3+2*(a^2)/(b^2)+3*(a^4)/(b^4)));
E=36000;tm=1;b=10;r=3.8;u=0.3;f=0.24;d=0.5;
% D=(E*(d^3))/12*(1-u^2);
% Mx=-D*(diff(w,'x',2)+u*diff(w,'y',2));
% M_y=-D*(diff(w,'y',2)+u*diff(w,'x',2));
% % M=diff(w,'x',1);
% % M_xy=-D*(1-u)*(diff(M,'y',1));
% % M_yx=-D*(1-u)*(diff(M,'y',1));
% sigma_x=(12*(Mx)*z)/(d^3);
% sigma_y=(12*(M_y)*z)/(d^3);
% % shear_xy=12*(M_xy)*z/(d^3);
% % shear_yx=12*(M_xy)*z/(d^3);
% % sigma_max=(0.123*(a^2)*(b^2)*(pi^2)*q*((a^2)+u*(b^2)))/((d^2)*(3*(a^4)+2*(a^2)*(b^2)+3*(b^4)));
% x=0;y=0;z=-0.25;
% sigma_t=-(140737488355328*E*a^3*b^2*r*z*pi^2*(u^2 - 1)*(b^2*cos((2*pi*x)/a) + b^2*cos((2*pi*x)/a)*cos((2*pi*y)/b) + a^2*u*cos((2*pi*y)/b) + a^2*u*cos((2*pi*x)/a)*cos((2*pi*y)/b)))/(1713638851887625*D*f*(3*a^4 + 2*a^2*b^2 + 3*b^4));
E=36000;tm=1;b=10;r=3.8;u=0.3;f=0.24;d=0.5;D= 3.412500000000000e+02;s=0;
% sigma_t_max=(0.203*E*a^3*b^2*r*(2*u*a^2 + 2*b^2)*(u^2 - 1))/(tm*(D*f*(3*a^4+ 2*a^2*b^2 + 3*b^4)));

https://wenku.baidu.com/view/7714de35bf23482fb4daa58da0116c175e0e1e6a.html