betta=0.9;
omegga=0.2;
fi=0.75;
siggma=1;
etta=1;
xi=1;
epsion=10;
thetta=1;
rho_a=0.95;
rho_r=0.9;
fi_pi=1.5;
fi_y=0.7;
tau_n=0.05;
tau_c=0.15;
r=0.1;

pi=betta*(1+r)-1;

A=1;
x0=[1,1,1,1];
x = nk_commitment_ss(x0,pi,betta,omegga,fi,siggma,etta,xi,epsion,thetta,tau_n,tau_c);
pi_star=x(1);
w=x(2);
v=x(3);
N=x(4);
mc=w;
Y=A*N/v;
G=omegga*Y;
C=Y-G;
x1=1/(1-fi*betta)*C^(-siggma)*mc*Y;
x2=1/(1-fi*betta)*C^(-siggma)*Y;