function[f]=foc_abierta(x)
global alpha beta gamma omega deltak deltag sigma psi d_b r_e
c_ss= x(1);
n_ss = x(2);
lambda_ss= x(3);
k_ss = x(4);
g_ss = x(5);
i_ss=x(6);
ig_ss=x(7);
y_ss=x(8);
q_ss=x(9);
miu_ss=x(10);


 
f(1)=(c_ss-omega^(-1)*(n_ss)^(omega))^(-sigma) - lambda_ss;
f(2)=-(c_ss-omega^(-1)*n_ss^(omega))^(-sigma)*n_ss^(omega-1) - (1-alpha-gamma)*lambda_ss*k_ss^alpha*g_ss^gamma*n_ss^(-alpha-gamma);
f(3)=-miu_ss+beta*miu_ss*((1-deltak)+(psi/2)*deltak^2)+beta*lambda_ss*(alpha*k_ss^(alpha-1)*g_ss^(gamma)*n_ss^(1-alpha-gamma));
f(4)=-miu_ss+beta*miu_ss*((1-deltag)+(psi/2)*deltag^2)+beta*lambda_ss*(gamma*k_ss^(alpha)*g_ss^(gamma-1)*n_ss^(1-alpha-gamma));
f(5)=y_ss-k_ss^alpha*g_ss^gamma*n_ss^(1-alpha-gamma);
f(6)=-i_ss+k_ss-(1-deltak)*k_ss;
f(7)=-ig_ss+g_ss-(1-deltag)*g_ss;
f(8)=-y_ss+c_ss+i_ss+ig_ss+r_e*d_b;
f(9)=-lambda_ss+miu_ss*(1+psi*deltak);
f(10)=(lambda_ss/miu_ss)-q_ss;

end
