function table = moments(economy)

eval(['load ' economy '_run.mat'])


%concatenate controls and states 
variables_cu = [controlvar_cu statevar_cu];

%Var-Cov matrix of innovation to state vector
varshock = ETAMATRIX*ETAMATRIX';

%Compute var-cov matrix of x_t
[VARY,VARX]=mom([gx;eye(size(hx))],hx,varshock,0);

ntaud = find(variables_cu=='taud_cu');
nr = find(variables_cu=='r_cu');
npai = find(variables_cu=='pai_cu');
noutput = find(variables_cu=='output_cu');
nh = find(variables_cu=='h_cu');
nc = find(variables_cu=='c_cu');
nw = find(variables_cu=='w_cu');
nvariables = [ntaud;nr;npai;noutput;nh;nc;nw];

var_variables = diag(VARY)*100^2;
var_variables([nr npai]) = var_variables([nr npai])*16;
std_variables = var_variables.^(1/2);

table=zeros(7,3);
table(:,1) = std_variables(nvariables);

eval(['save ' economy '_moments.mat'])

disp('rows: taud r pai output h c w')