%
% ex6.m is an example matlab program 
% to solve the rbc model by muck
%

% ----- save output to file ----- %

clear all
diary ex6.mout
diary off
delete ex6.mout
diary ex6.mout

% ----- display date and time of computation ----- %

%format long
date
time0 = clock;


% ----- paramters ----- %

alpha = 1/3 ;
delta = .10;
beta  = 0.99;


% ----- steady states ----- %

cbar = 1.2132;
kbar = 5.2678;


% ----- set up sytem as required by solvem.m ----- %

AA = -kbar;
BB = alpha*kbar^alpha + (1-delta)*kbar ;
CC = -1;
DD = 0;


FF = 0;

GG = -1/cbar*beta*(alpha*(alpha-1)*kbar^(alpha-1)) ;
HH = 0;
JJ = 1/cbar*beta*(alpha*kbar^(alpha-1)+(1-delta)) ;
KK = -1/cbar;
LL = 0;
MM = 0;


NN = 0;


[PP,QQ,RR,SS] = solvem(AA,BB,CC,DD,FF,GG,HH,JJ,KK,LL,MM,NN)


comptime = etime(clock, time0)
diary off