function x = nk_commitment_ss(x0,pi,betta,omegga,fi,siggma,etta,xi,epsion,thetta,tau_n,tau_c)
x = csolve(@nk_commitment,x0,[],1e-8,100,pi,betta,omegga,fi,siggma,etta,xi,epsion,thetta,tau_n,tau_c);
end

function F = nk_commitment(x,pi,betta,omegga,fi,siggma,etta,xi,epsion,thetta,tau_n,tau_c)

F=[(1+pi)^(1-epsion)-(1-fi)*(1+x(1))^(1-epsion)-fi;
 (1+x(1))*(1-fi*betta*(1+pi)^(epsion))-epsion/(epsion-1)*(1+pi)*x(2)*(1-fi*betta*(1+pi)^(epsion-1));
 x(3)*(1-(1+pi)^(epsion)*fi)-(1-fi)*(1+x(1))^(-epsion)*(1+pi)^(epsion);
 xi*x(4)^(etta)*x(3)^(-siggma)-(1-tau_n)/(1+tau_c)*(1-omegga)^(-siggma)*x(4)^(-siggma)*x(2);];

end