% 一橋大学「金融経済論I」（２０１０年度、担当：塩路)
%
% ファイル名＝rbc2_model.m
% バージョン2
% 2010年11月16日
%
% 目的：rbc2.m(バージョン2)に付随するコード
% このファイルでモデルを定式化する

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% カナダ中央銀行（→現内閣府）の西山慎一氏には、バージョン1のチェックをしていただき、
% その誤りを指摘、修正していただいた。
% ここに記して感謝したい。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% 2012年11月27日、誤りを訂正


%% Mofification Note:
%% line 7: Modified to following la=x(1) ;ca=x(2) ; ka=x(3) ;  va=x(4) ; lt=x(5) ; ct=x(6) ; kt=x(7) ; vt=x(8) ;

function optcon=f1(x);
% Within "x, endogenous vars are in the order of ('a' denotes t+1 values):
       %  la ca   ka  va  lt ct   kt  vt   
la=x(1) ;ca=x(2) ; ka=x(3) ;  va=x(4) ; lt=x(5) ; ct=x(6) ; kt=x(7) ; vt=x(8) ;

global sigma alpha myu beta delta lamda phi ;


% Eliminate Price  
       ra = (va*alpha*(ka/la)^(alpha-1)) ;
       wt = (1-alpha)*vt*(kt/lt)^alpha;
   % Optimal Conditions  & state transition
       labor   = lt^lamda-wt/(myu*ct^sigma) ;       % LS = LD 
       euler   = ct^(-sigma) -(ca^(-sigma))*beta*(1+ra-delta) 	  ;     % C-Euler
       capital = ka - (1-delta)*kt-vt*(kt^alpha)*(lt^(1-alpha))+ct ; % K-trans
       tech    = (va-1) - phi*(vt-1);        
       
optcon  = [labor ; euler  ; capital  ; tech];    
