%GENTS      Print output with the time series generated from Fortran
%           code AEA3.F90.  

%           Ellen McGrattan, 4-6-02

load aea3.dat
aea    = aea3;
s      = aea(:,1);
k      = aea(:,2);
cp     = aea(:,3);
hr     = aea(:,4);
ip     = aea(:,5);
y      = aea(:,6);
w      = aea(:,7);
r      = aea(:,8);
g      = aea(:,9);
taul   = aea(:,10);
tauk   = aea(:,11);
tauc   = aea(:,12);
taux   = aea(:,13);
z      = aea(:,14);
delta  = aea(1,15);
gn     = aea(1,16);
gz     = aea(1,17);
t      = [1:length(s)];

disp('1.  Gross Domestic Product')
disp('2.  Private Consumption')
disp('3.  Public Consumption')
disp('4.  Gross Fixed Capital Formation')
disp('5.  Beginning-of-period Capital Stock')
disp('6.  Depreciation')
disp('7.  Total Hours')
disp('8.  Wages and salaries')
disp('9.  Propreitors'' income')
disp('10. Corporate Profits')
disp('11. Indirect Business Tax')
disp('12. Rental income ')
disp(' ')

comp = w.*hr;
pwinc= .125* comp;
ws   = comp-pwinc;

rinc = r.*k;
dep  = .06*k;
princ= .25*pwinc;
cprf = rinc/4;
rent = rinc-princ-cprf-dep;

gdp  = cp.*(1+tauc)+ip+g;

tr   = ((1+gn)*(1+gz)).^[-50:(-50+length(t)-1)]';
ntr  = ((1+gn)).^[-50:(-50+length(t)-1)]';
X    = [gdp.*tr,cp.*tr,g.*tr,ip.*tr,k.*tr,.06*k.*tr,hr.*ntr, ...
        ws.*tr,(pwinc+princ).*tr,cprf.*tr,(tauc.*cp).*tr,rent.*tr];

X    = X(51:150,:);
disp(sprintf([' %6.3f %6.3f %6.3f %6.3f %6.3f %6.3f %6.3f %6.3f %6.3f', ...
              ' %6.3f %6.3f %6.3f\n'],X'))
