新人求教，望大神帮忙看一眼！求不出q和a

clear
l=input('AB跨度(/m)=')*1e3;
sigma_max=input('最大正应力(/MPa)=');
Wz=input('25a号工字钢截面系数(/cm3)=')*1e3;
%求A、B处支座反力
syms q a Ma Mc sigma_a sigma_c fun1 fun2 fun3 fun4
[Fa,Fb]=solve('Fa+Fb-q*(l+2*a)=0','Fa*a+Fb*(l+a)-q*(l+2*a)^2/2=0','Fa','Fb');
display(Fa);
display(Fb);
%求A、B、C处截面弯矩值
Ma=q*a^2/2;
Mb=Ma;
Mc=q*l^2/8-q*a^2/2;
%求A、B、C截面处正应力
sigma_a=Ma/Wz;
sigma_b=sigma_a;
sigma_c=Mc/Wz;
display(sigma_a);
display(sigma_b);
display(sigma_c);
%为使支座处截面A、B上及跨度中央截面C上的最大正应力均为140MPa
[q,a]=solve('sigma_a-sigma_max==0','sigma_c-sigma_max==0','q','a');