function [z1,z2,z3] = hwlogS_SOEML_steady_state_helper(delta,omega,kappa,beta,phi,alphaN,alphaT,rho,psii,bbar,theta,rho_AT,rho_AN)

z = fsolve(@(x) SS_equations_R1(x,delta,omega,kappa,beta,theta,phi,alphaN,alphaT,rho,psii,bbar),...
    [1;1;1],optimset('Display','Off','TolFun',1.0e-08));
z1=z(1);
z2=z(2);
z3=z(3);