function [ys,check] = soe1_steadystate(ys, junk)

global M_;

growth_us_ss = get_param_by_name('growth_us_ss'); 

% because of that crapy (i am sorry) Koopman&durbin stuff, which doesn't
% take into account cointagration of unit-root variables, we initialize
% all unit root variables to the first observation

data;



Y_US = 0;
LGDP_US = LGDP_US(1);
LGDP_US_BAR = LGDP_US;
G_US = growth_us_ss;
E1_Y_US     = 0;
GROWTH_US = growth_us_ss;
GROWTH4_US = growth_us_ss;
GROWTH4_US_BAR = growth_us_ss;



check = 0;

ys = [
Y_US 
LGDP_US 
LGDP_US_BAR 
G_US 
E1_Y_US
GROWTH_US 
GROWTH4_US 
GROWTH4_US_BAR
     ];
