


function [ F ]  = find_steady_env(Z,eta,erow,gamma1,gamma2,tau,theta1,theta2,d0,d1,d2,d3,a,epsilon,alpha,rk,delta,kappaL,sigma,g,phi)
y=Z(1);
x=Z(2);
h=Z(3);

e=x*(1-eta)-erow;
mu=(gamma1*tau/(theta1*theta2)*y^(-gamma2))^(1/(theta2-1));
D=d3*(d0+d1*x+d2*x^2);
A=(1-d3*(d0+d1*x+d2*x^2))*a;
mc=(epsilon-1)/(epsilon)-(1-gamma2)*gamma1*tau*y^(-gamma2)*(1-mu)-theta1*mu^(theta2);
k=alpha*mc*y/rk;             
i=delta*k;
w=(1-alpha)*mc*y/h;          
c=y-i-g-y*theta1*mu^(theta2);                    
                                     
F=[y-A*(k^(alpha)*h^(1-alpha))
    kappaL-w/(c^(sigma)*h^(phi))
    e-(1-mu)*gamma1*y^(1-gamma2)
    ];

end

