function pii=firm_vol2_2_steady_state_helper(cai_0, cai_1, beta, delta, eta, alpha, psi, gamma, lambda)
options=optimset('Display','Final','TolX',1e-4,'TolFun',1e-4);
 pii=fsolve(@(pii)(((cai_0*(1-cai_1)*((1/beta)-1)/(1-delta))*(1/((delta/cai_0)^(1/(1-cai_1))-(delta/cai_0)^(cai_1/(1-cai_1)))))^(1-cai_1))*((psi/(1-alpha))^(1/(1+eta))*((1-beta*(1-delta))/(alpha*beta^2))^(alpha/(eta*(1-alpha))))*((alpha*(1-beta)*(lambda-1)+(1-alpha)*((beta-gamma)+lambda*(1-beta*(1-delta))))*(pii*(1-beta*(1-delta)-alpha*beta)+((1-beta*(1-delta)-delta*alpha*beta^2))*(alpha*beta*lambda-pii*(beta-gamma+lambda*(1-beta*(1-delta)))))/(alpha*beta*(lambda*(1-beta*(1-delta))-beta*(beta-gamma+lambda*(1-beta*(1-delta))))))^(1/(1+eta)),3.6,options);