% This code is to calibrate mug, mui, mux and delta1 by using steady state of endogenous variable %
function [MUG, MUI, MUX, DELTA1] = other_params

global BETTA SIG ALPA DELTA0 DELTA1 MUA MUY MUI;
global DELTA2 THETA_C THETA_L RHO_Z RHO_G RHO_XG RHO_X5 RHO_A KAPA CHI;

ys = 0;
junk = 0;

% [ys, check] = cycle_steadystate(junk,ys);
[ys, check] = cycle_uniform_3_steadystate(junk,ys);
% ys = [C ex1 ex2 ex3 ea1 ea2 ea3 ez1 ez2 ez3 eg1 eg2 eg3 G ga_obs gai_obs gc_obs gg_obs gy_obs H h_obs I K L LAMDA MU_A MU_K MU_X MU_Y Q Q_A U X_G Y Z]'
% ys(17) = G, ys(19) = H, ys(22) = K, ys(26) = MU_K, ys(27) = MU_X, ys(29) = Q, ys(30)= U, ys(34) = Z
MUG = ys(17);
MUI = ys(26);
MUX = ys(27);
DELTA1 = ys(34)*ALPA*(ys(30)/ys(26))^(ALPA-1)*(ys(22)/ys(19))^(ALPA-1)/ys(29);