%% define the indices for states and controls

% s = [psi(-1), mu(-1), rn, log(xi), LOG q^u, LOG q^u_{-1}]
psilag_idx      = 1;
mulag_idx       = 2; 
rn_idx          = 3;  
u_idx           = 4;
logxi_idx       = 4+1;
logqu_idx       = 5+1;
loglagqu_idx    = 6+1;


% x = [ygap, pi, i, psi, lambda]
ygap_idx    = 1;
pi_idx      = 2;
i_idx       = 3;
psi_idx     = 4;
lambda_idx  = 5;

% e = [e_rn, e_xi]; 
e_rn_idx = 1;
e_u_idx = 2;
e_name = {'rn','u'}; % corresponding names