function [vecss, SteadyState] = deterministic_ss(parameters)
%parameters
beta= parameters(1);
phi= parameters(2);
beta_E = parameters(3);
alpha= parameters(4);
delta = parameters(5);
theta = parameters(6);
eps = parameters(7);
N = parameters(8);
R_b_s = parameters(9);
Rho_r = parameters(10);
K_pi = parameters(11);
K_y = parameters(12);
a = parameters(13);
b = parameters(14);
pi_meta = parameters(15);
y_pot = parameters(16);
m_A = parameters(17);
m_B = parameters(18);
P = parameters(19);
q = parameters(20);
f3 = parameters(21);


%start to solve the steady state with z = q = pi = f3 = 1


%Households
R = 1/beta;
l = 1 - phi/(c * w *(1 - tau));

%Entrepreneur's A
b_A = q * k_A - beta_E * nw_A;
nw_A = alpha * (y_A/x_A) + q * (1-delta) * k_A - R_b_s * b_A;
y_A = 1*((k_A)^alpha) * ((l_A)^(1-alpha));
w = ((1-alpha)* y_A)/(x_A * l_A);
x_A = P/P_A;
lamb_2_A = (1-(R_b_s * beta_E))/c_A;
c_A = (1- beta_E)*nw_A;
b_A = (m_A * q * k_A * (1-delta))/R_b_s;
q = beta_E * ((alpha * y_A)/k_A * x_A) + c_A * lamb_2_A * m_A * (q * (1-delta)/R_b_s);


%Entrepreneur's B
b_B = q * k_B - beta_E * nw_B;
nw_B = alpha * (y_B/x_B) + q * (1-delta) * k_B - R_b_B * b_B;
y_B = 1*((k_B)^alpha) * ((l_B)^(1-alpha));
w = ((1-alpha)* y_B)/(x_B * l_B);
x_B = P/P_B;
lamb_2_B = (1-(R_b_B * beta_E))/c_B;
c_B = (1- beta_E)*nw_B;
b_B = (m_B * q * k_B * (1-delta))/R_b_B;
q = beta_E * ((alpha * y_B)/k_B * x_B) + c_B * lamb_2_B * m_B * (q * (1-delta)/R_b_B);

%Retailers

MgC = (p_A/(a + b * (p_A * b/p_B * a)^(-rho/rho-1)) + p_B/(a * ( a * p_B/b * p_A)^(rho/1-rho) + b))^(-eps);
mc = MgC/P;
f1 = y * mc/(c * (1 - beta * theta));
f2 = y/( c * (1 - beta * theta));
% f2 = f1 * eps/(1-eps);                            As equações de cima já me permitem encontrar f1 e f2
P_otimo = P;


%Capital producer
K = i/delta;

%Banking Sector - Cournot (DELTA = beta)
R_b_B = R/(1 - (1/N) * (1/PED));
PED = 1 + (1/(1-alpha))*((m_B * (1 - delta))/(R_b_B * beta * MPK));
MPK = alpha * 1 * ((l_B)^(1-alpha)) * ((k_B)^(alpha-1)) * (1/x_B);

%State-owned bank
(R - R_b_B) * b_A = tau * w * (l_A + l_B);

%Central Bank
%R = Rho_r * R + (1 - Rho_r)*(R + K_pi * (1 - pi_meta) + K_y * (y - y_pot));     Não preciso pois R já sai do household        

% Market Clearing
c + c_A + c_B + i = y;
y = ((a * ((y_A)^rho) + b * ((y_B)^rho))^(1/rho));
l = l_A + l_B;
K = k_A + k_B;



%steady state vector
vecss = [R, c, w, b_A, k_A, nw_A, p_A, l_A, c_A, y_A, x_A, lamb_2_A, b_B, k_B, nw_B, p_B, l_B, c_B, y_B, x_B, lamb_2_B, MgC, mc, f1, f2, P_otimo, K, i, tau, PED, MPK, R_b_B, l, y];

%struct steady vector
SteadyState = struct();

SteadyState.R = R;
SteadyState.c = c;
SteadyState.w = w;
SteadyState.b_A = b_A;
SteadyState.k_A = k_A;
SteadyState.nw_A = nw_A;
SteadyState.p_A = p_A;
SteadyState.l_A = l_A;
SteadyState.c_A = c_A;
SteadyState.y_A= y_A;
SteadyState.x_A= x_A;
SteadyState.lamb_2_A= lamb_2_A;
SteadyState.b_B = b_B;
SteadyState.k_B = k_B;
SteadyState.nw_B = nw_B;
SteadyState.p_B = p_B;
SteadyState.l_B = l_B;
SteadyState.c_B = c_B;
SteadyState.y_B = y_B;
SteadyState.x_B= x_B;
SteadyState.lamb_2_B= lamb_2_B;
SteadyState.MgC = MgC;
SteadyState.mc = mc;
SteadyState.f1 = f1;
SteadyState.f2 = f2;
SteadyState.P_otimo = P_otimo;
SteadyState.K = K;
SteadyState.i = i;
SteadyState.tau = tau;
SteadyState.PED = PED;
SteadyState.MPK = MPK;
SteadyState.R_b_B = R_b_B;
SteadyState.l = l;
SteadyState.y = y;