%===============================================================================
% Variables
%===============================================================================

% Growth Factors
var n g;

% Trend
trend_var(growth_factor=n) N;
trend_var(growth_factor=g) A;

% Endogenous
var(deflator=N*A) Y K C I;
var(deflator=N) L;
var(deflator=A) w;
var r;

% Predetermined
predetermined_variables K;

%===============================================================================
% Parameters
%===============================================================================

parameters hbar delta alpha beta gamma;

hbar  = 52*100;  % page 430
delta = 0.0556;  % page 433
alpha = 0.3590;  % page 434
beta  = 0.9752;  % page 440
gamma = 0.2846;  % page 440

%===============================================================================
% Model
%===============================================================================

model;
    Y = A*(K^alpha)*(L^(1-alpha));            % Equation (3)
    w = (1-alpha)*A*(K^alpha)*(L^(-alpha));   % Equation (4)
    r = alpha*A*(K^(alpha-1))*(L^(1-alpha));  % Equation (5)
    C + I = Y;                                % Equation (7)
    K(+1) = (1-delta)*K + I;                  % Equation (8)
    w*(hbar*N - L) = ((1-gamma)/gamma)*C;     % Equation (17)
    C(+1)/C = beta*(1 - delta + r(+1));       % Equation (18)
    n = N/N(-1);
    g = A/A(-1);
end;

%===============================================================================
% Initial Values
%===============================================================================

initval;
    %N = 3235000;
    %A = 364.98555;
    Y = 80810000000;
    K = 246272933253;
    L = 4354000;
    C = 56403546606;
    I = 24406453394;
    r = alpha*(Y/K);
    w = (1-alpha)*(Y/L);
end;

%===============================================================================
% End Values
%===============================================================================

endval;
    n = 1.0028662;
    g = 1.0238426;
    Y = 80810000000;
    K = 246272933253;
    L = 4354000;
    C = 56403546606;
    I = 24406453394;
    r = alpha*(Y/K);
    w = (1-alpha)*(Y/L);
end;

steady;
check;