function [ys,check] = nkopenadjbankpruramsey_steadystate(ys,exo);

chi = -0.013;
sigma = 1;
varphi = 3;
alpha = 0.3;
varepsilon = 6;
thetaf = 0.75;
delta = 0.025;
omega = 0.4;
thetab = 0.97; 
xi = 0.4; 
upsilon = 0.00376299619999994; 
c = 0.005;

epsr = 0;
pih = 0;
pir = 0;
pi = 0;
Pch = 1;
DE = 1;
H = 1;
d = 1;
Q = 1;
RF = 1.01316908;
R = RF;
Rk = R + 0.01; 
RkF = Rk;
Pk = Rk - (1 - delta);
beta = 1/R;
C = (beta^(1/chi)) - 1;
mc = (varepsilon -1)/varepsilon;
W = (mc*((alpha^alpha)*((1 - alpha)^(1 - alpha)))*Pk^(-alpha))^(1/(1-alpha));
kl = alpha*W/(Pk*(1 - alpha));
L = (W*(C^-sigma))^(1/varphi);
K = kl*L;
SH = K;
SF = ((1 - xi)/xi)*SH;
I = delta*K;
Y = (K^alpha)*(L^(1 - alpha));
BF = (Y + (RkF - 1)*SF - C - I)/(RF - 1);
X = Y - (1 - omega)*C - I;
x1 = ((C^(-sigma))*Y*mc)/(1 - beta*thetaf);
x2 = ((C^(-sigma))*Y)/(1 - beta*thetaf);
V = (log(C) - (L^(1+varphi))/(1+varphi))/(1 - beta);
RkF = Rk;
varpi = 2;
nu = (beta*(1 - thetab)*(Rk - (1 + taus)*R - c))/(1 - beta*thetab*(R + varpi*(Rk - (1 + taus)*R - c)));
eta = nu*(R/(Rk - (1 + taus)*R - c));
No = thetab*((Rk - R - c)*varpi + R)*SF/(xi*varpi);
Nn = upsilon*SF/xi;
N = SF/(xi*varpi);
D = (1 + taus)*(SF/xi) - (1 + taun)*N;

check = 0;

ys(1) = R;
ys(2) = RF;
ys(3) = Pk;
ys(4) = Rk;
ys(5) = log(Q); 
ys(6) = log(Y);
ys(7) = log(C); 
ys(8) = log(I);
ys(9) = pir;
ys(10) = pih;
ys(11) = pi;
ys(12) = log(Pch);
ys(13) = log(x1);
ys(14) = log(x2);
ys(15) = log(mc);
ys(16) = log(K);
ys(17) = log(L);
ys(18) = log(W);
ys(19) = log(d);
ys(20) = log(DE);
ys(21) = log(H);
ys(22) = log(BF);
ys(23) = log(beta);
ys(24) = X;
ys(25) = log(SF);
ys(26) = log(varpi);
ys(27) = log(nu);
ys(28) = log(eta);
ys(29) = RkF;
ys(30) = log(No);
ys(31) = log(Nn);
ys(32) = log(N);
ys(33) = log(D);
ys(34) = V;
ys(37) = epsr;



 
 
 


 


