syms x1 x2 u1 u2f1=0.01*u1*(1-x1)-7.2*1e4*x1*exp(-1e4/x2); f2=0.01*u1*(350-x2)+1.44*1e13*x1+0.01*u2*(1-exp(-700/u2))*(350-x2);[x1;x2]=solve([f1;f2],u1,u2);
