%CASES   Compute the steady state for alternative parameterizations of model 3
%        with data for 2007.
%

%        Ellen McGrattan, 4-11-11
%        Revised, ERM, 6-20-11

disp('Alternative parameterizations: ')
disp('   1. Lump sum taxes')
disp('   2. Barro-Sahasakul tax rate on labor; ')
disp('      NIPA corporate income tax on profits ')
disp('   3. McGrattan tax rates')
disp(' ')
disp('Note: other parameters set to match NIPA variables and hours')
disp(' ')
disp(' ')
icase   = input('Which case?  ');

%
% NIPA data and hours
%

gdp     = 1;
c       = 0.7819;
g       = 0.0000;
yn      = 0.3244;
xt      = 0.1020;
xn      = 0.1161;
h       = 0.2777;
hn      = 0.0443;
kt      = 1.4920;
NIPAbw  = 0.4456;


%
% Other parameters
%
gn      = 0.01;
gz      = 0.019;
phi     = 0.00001;
q       = 1;
chi     = 0.5;
deli    = 0.03917652;
beta    = 0.979;
isymm   = 1;


%
% Tax rates used in different cases 
%

if icase==1;  % Lump sum taxes
  
  tauc    = 0.;
  tauh    = 0.;
  taux    = 0.;
  tauk    = 0.;
  taup    = 0.;
  tauu    = 0.;
  taud    = 0.;

end;

if icase==2; % Barro-Sahasakul tax rate (on labor income)
  
  tauc    = 0.;
  tauh    = 0.3530;
  taux    = 0.;
  tauk    = 0.;
  taup    = 0.2921;
  tauu    = 0.;
  taud    = 0.;

end;

if icase==3; % McGrattan's tax rates
  
  tauc    = 0.0620;
  tauh    = 0.3580;
  taux    = 0.;
  tauk    = 0.0136;
  taup    = 0.2862;
  tauu    = 0.;
  taud    = 0.078;

end;

param  = [gdp-yn,yn,c,xt,xn,g,kt,h,hn,NIPAbw,tauc,tauh,tauh,taux,tauk,taup, ...
          taud,gn,gz,phi,q,chi,deli,beta,isymm];
if (isymm~=1);
  param = [param',alph2,thet2]';
end;

target  = mod3ss(param);
