mkdir('./result');
clear; clc;
dynare YYY.mod; % no
AA = oo_;
save('./result./AA.mat', 'AA');
save('./result./M_.mat', 'M_');

clear; clc;
dynare YYY2.mod; % Qhat
BB = oo_;
save('./result./BB.mat', 'BB');

clear; clc;
dynare YYY3.mod; % Qhat-Qhat(-1)
CC = oo_;
save('./result./CC.mat', 'CC');

clear; clc;
dynare YYY1.mod; % Bhat
DD = oo_;
save('./result./DD.mat', 'DD');
% --------------
clear; clc;
load('./result./AA.mat');
load('./result./BB.mat');
load('./result./CC.mat');
load('./result./DD.mat');
load('./result./M_.mat');
ypos     = strmatch('Ygap', M_.endo_names ,'exact');
pipos    = strmatch('Pigap', M_.endo_names ,'exact');
npos     = strmatch('LHgap', M_.endo_names ,'exact');
rpos     = strmatch('Rgap', M_.endo_names ,'exact');
bpos     = strmatch('Bgap', M_.endo_names ,'exact');
vol = round(sqrt([AA.var(ypos, ypos)  AA.var(pipos, pipos) AA.var(npos, npos) AA.var(rpos, rpos) AA.var(bpos,bpos);
                  BB.var(ypos, ypos)  BB.var(pipos, pipos) BB.var(npos, npos) BB.var(rpos, rpos) BB.var(bpos,bpos);
                  CC.var(ypos, ypos)  CC.var(pipos, pipos) CC.var(npos, npos) CC.var(rpos, rpos) CC.var(bpos,bpos);
                  DD.var(ypos, ypos)  DD.var(pipos, pipos) DD.var(npos, npos) DD.var(rpos, rpos) DD.var(bpos,bpos);
                  
                  ])*100, 3);
wel = round([
    AA.var(pipos,pipos)+1*AA.var(ypos,ypos)+0.1*AA.var(bpos,bpos) ;
    BB.var(pipos,pipos)+1*BB.var(ypos,ypos)+0.1*BB.var(bpos,bpos) ;
    CC.var(pipos,pipos)+1*CC.var(ypos,ypos)+0.1*CC.var(bpos,bpos) ;
    DD.var(pipos,pipos)+1*DD.var(ypos,ypos)+0.1*DD.var(bpos,bpos) ;
    ]*10000, 3);
vol_wel = [vol wel];
save('./result./vol_wel.mat', 'vol_wel');