e1=sym('(0.015*.85*35+0.985*1.15*25)/(0.985*1.15^2*25+0.003*0.85^2)=M/T');
e2=sym('(0.015*.85*35+0.985*1.15*25)/(0.985*1.15^2*25+0.003*0.85^2)=F');
e3=sym('(25*(a+0.85*b*F)(a*T+0.85*F*T)+45*(a+1.15*b*F))/((a+0.85*b*F)*45+Q)=400');
e4= sym('a*T+1.15*b*F*T=Z');
e5= sym('(1/(20*pi)^0.5*exp*(-(Z-M)^2/(2*(0.05*M)^2)=Q');
[a b T M F Z]=solve(e1,e2,e3,e4)