%
% Table 3B: Balanced Growth Aggregate Statistics
%   Permanent Change in Net Tax Function T^w(y)
%      

nv           = 13;
tab2a        = zeros(nv,4);
!rm -f trantf.dat

for icase=1:4;
  if (icase>=1 & icase<=2);
    GNP0    = 0.7450;
    C0      = 0.7437*GNP0;
    XT0     = 0.2123*GNP0;
    XI0     = 0.1346*GNP0;
    LINC0   = 0.5848*GNP0;
    CINC0   = 0.4152*GNP0;
    L0      = 0.2936;
    KTeoy0  = 4.1396*GNP0;
    KIeoy0  = 1.7236*GNP0;
    W0      = 1.4841;
    REV0    = 0.3921*GNP0;
    NWhh0   = 5.5743*GNP0;
  end;
  if (icase>=3 & icase<=4);
    GNP0    = 0.7465;
    C0      = 0.7392*GNP0;
    XT0     = 0.2135*GNP0;
    XI0     = 0.1346*GNP0;
    G0      = 0.0436*GNP0;
    LINC0   = 0.5848*GNP0;
    CINC0   = 0.4152*GNP0;
    L0      = 0.2903;
    KTeoy0  = 4.1649*GNP0;
    KIeoy0  = 1.7234*GNP0;
    W0      = 1.5039;
    REV0    = 0.3931*GNP0;
    NWhh0   = 5.5992*GNP0;

  end;

  if (icase==1);
    !cp ../benchmark/continue.dat trantf.dat
  end;
  if (icase==2);
    !cp ../benchmark/nomedss.dat trantf.dat
  end;
  if (icase==3);
    !cp ../standin/continue.dat trantf.dat
  end;
  if (icase==4);
    !cp ../standin/nomedss.dat trantf.dat
  end;
  clear AVG
  load trantf.dat
  plottrantf
  tab2b(:,icase) = ...
             [[C(240)/C0-1, ...
               (XT1(240)+XT2(240))/XT0-1, ...
               (XI1(240)+XI2(240))/XI0-1, ...
               w(240)*(L1(240)+L2(240))/LINC0-1, ...
               (cprof(240)+uprof(240))/CINC0-1, ...
               GDP(240)/GNP0-1, ...
               (L1(240)+L2(240))/L0-1, ...
               (KT1eoy(240)+KT2eoy(240))/KTeoy0-1, ...
               (KI1eoy(240)+KI2eoy(240))/KIeoy0-1, ...
               irate(240)-1, ...
               w(240)/W0-1, ...
               Rev(240)/REV0-1, ...
               NWhh(240)/NWhh0-1]*100];
  !rm -f trantf.dat
end;
tab2b(:,6) = tab2b(:,4)-tab2b(:,3);
tab2b(:,5) = tab2b(:,4);
tab2b(:,4) = tab2b(:,3);
tab2b(:,3) = tab2b(:,2)-tab2b(:,1);

clc
disp('                 TABLE 3B. CHANGES IN BALANCED GROWTH AGGREGATE STATISTICS              ')
disp('                       Workers'' Net Tax Function Changed Permanently')
disp('----------------------------------------------------------------------------------------')
disp('                             Heterogeneous Within Cohort      Homogeneous Within Cohort ')
disp('                             Continue  Reform  Difference    Continue  Reform Difference')
disp(' Statistic/Economy:            (1)      (2)     (2)-(1)        (4)      (5)    (5)-(4)  ')
disp('----------------------------------------------------------------------------------------')
disp(sprintf('  GNP                       %8.1f %8.1f %8.1f    %8.1f %8.1f %8.1f',tab2b(6,:)));
disp(sprintf('  Consumption               %8.1f %8.1f %8.1f    %8.1f %8.1f %8.1f',tab2b(1,:)));
disp(sprintf('  Tangible investment       %8.1f %8.1f %8.1f    %8.1f %8.1f %8.1f',tab2b(2,:)));
disp(sprintf('  Intangible investment     %8.1f %8.1f %8.1f    %8.1f %8.1f %8.1f',tab2b(3,:)));
disp(sprintf('  Labor income              %8.1f %8.1f %8.1f    %8.1f %8.1f %8.1f',tab2b(4,:)));
disp(sprintf('  Capital income            %8.1f %8.1f %8.1f    %8.1f %8.1f %8.1f',tab2b(5,:)));

disp(sprintf('  Tangible capital          %8.1f %8.1f %8.1f    %8.1f %8.1f %8.1f',tab2b(8,:)));
disp(sprintf('  Intangible capital        %8.1f %8.1f %8.1f    %8.1f %8.1f %8.1f',tab2b(9,:)));
disp(sprintf('  Labor input               %8.1f %8.1f %8.1f    %8.1f %8.1f %8.1f',tab2b(7,:)));

disp(sprintf('  Wage rate                 %8.1f %8.1f %8.1f    %8.1f %8.1f %8.1f',tab2b(11,:)));
disp(sprintf('  Interest rate, level (%%)  %8.1f %8.1f %8.1f    %8.1f %8.1f %8.1f',tab2b(10,:)));
disp(sprintf('  Tax revenues              %8.1f %8.1f %8.1f    %8.1f %8.1f %8.1f',tab2b(12,:)));

disp(sprintf('  Household net worth       %8.1f %8.1f %8.1f    %8.1f %8.1f %8.1f',tab2b(13,:)));
disp('----------------------------------------------------------------------------------------')
disp(' Note: Statistics are percent changes between the final and initial balanced growth paths,')
disp('       with the exception of the interest rate which is the level on the final growth path.')

