clear all; close all;
dynare hmdy521s    

sdym       = oo_.endo_simul(46,1000:1230);
srm       = oo_.endo_simul(47,1000:1230);
sdrhpm       = oo_.endo_simul(48,1000:1230);
sdrspm       = oo_.endo_simul(49,1000:1230);
sdhcm       = oo_.endo_simul(50,1000:1230);

load hm.mat

figure(1)
plot(t1,sdym,'b-')
hold on
plot(t1,dym,'r')

figure(2)
plot(t1,srm,'b-')
hold on
plot(t1,rm,'r')

figure(3)
plot(t1,sdrspm,'b-')
hold on
plot(t1,drspm,'r')

figure(4)
plot(t1,sdrhpm,'b-')
hold on
plot(t1,drhpm,'r')

figure(5)
plot(t1,sdhcm,'b-')
hold on
plot(t1,dhcm,'r')
