function WR = script_v13WR_steadystate_helper(psi,LR,varpi,LAMBDAR,tau_l,PO,WR0)
    
        options = optimset('Display','Final','TolX',1e-10,'TolFun',1e-10);
        WR = fsolve(@(WR) psi*(LR^varpi)- (LAMBDAR*(1-tau_l)*WR^2)/PO,WR0, options);
   
end

