%
% dpgq3.m is a matlab file that solves for 
% the 'steady-state' for a simplified 
% ramsey model
%

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

clear all
diary dpgq3.mout
diary off
delete dpgq3.mout
diary dpgq3.mout

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

%format long
date
time0 = clock;

% ----- define the parameters ----- %

thetabar = 1.05;
alpha    = 0.30;
beta     = 0.99;
psi      = 1;
sigma    = 2;
delta    = 0.025 ;

% ----- dfine fsolve options ----- %

options = optimset('Display','iter','MaxIter',5000,'MaxFunEvals',5000);

% ----- solve for steady states ----- %

sol1 = fsolve(@dpgq3F,[1,1,.3],options,thetabar,alpha,beta,psi,sigma,delta)

% ----- exit ----- %

comptime = etime(clock, time0)
diary off