function Pi = Hansen2020_NotLinearized_steadystate_helper(Pi0, A, bet, thetaP, psiP, alph, chi, phi, phiPi, phiY, Tp)
options = optimset('Display', 'off', 'TolX', 1e-10, 'TolFun', 1e-10);
Pi = fsolve(@(Pi) Pi*(Pi-1) - bet*Pi*(Pi-1) - (thetaP/psiP)*(1/(1-alph) * chi * A^(-(1+chi)/(1-alph)) * bet^(-(phi+1)/(phiY*(1-alph))) * Pi^((1+phi)*(1-phiPi)/(phiY*(1-alph))) - (1+Tp)*(thetaP-1)/thetaP), Pi0, options);
