function F=steady_helper1(x)
load('_param.mat','beta','omega','theta','Rs','tau');
Phis=x;
F=tau*(1-beta*(1-omega*Phis))*Phis-(1-theta)*beta*((1-omega*Phis)/theta-Rs)-(1-theta);
end