% 一橋大学「金融経済論I」（2012年度、担当：塩路)
%
% ファイル名＝kiyotakimoore1_model.m
% 2013年1月22日
%
% 目的：kiyotakimoore1.mに付随するコード
% このファイルでモデルを定式化する

function eqcon=f1(x);
% ベクトルxは内生変数を次の順に並べたもの
       %  qa ka ba aa qt kt bt at
qa=x(1) ; ka=x(2) ; ba=x(3) ;  aa=x(4) ; qt=x(5) ; kt=x(6) ; bt=x(7) ; at=x(8) ;

global R v pai lamda phi astar psi  ;

% 均衡条件は次の4本の式からなる  
eq5 = qt-qa/R-ka+v ; %講義資料の（5）式

eq8 = ka-(1-pai)*lamda*kt - pai*((at+qt+lamda*phi)*kt-R*bt)/(phi+qt-(1/R)*qa) ;
    %講義資料の(8)式
    
eq9 = ba-R*bt -qt*(ka-kt) - phi*(ka-lamda*kt) + at*kt ;
    %講義資料の（9）式
    
eqx = (aa-astar) - psi*(at-astar) ;
    %生産性aはAR1過程に従う
       
eqcon  = [eq5 ; eq8 ; eq9 ; eqx] ;
