%FD1       plot results from fd.f.

%         Ellen McGrattan, 2-23-97


load fd.dat
n  = length(fd)/2;
x  = fd(1:n,1);
c  = fd(1:n,2);
c0 = fd(1:n,3);
v  = fd(n+1:2*n,2); 
v0 = fd(n+1:2*n,3); 
clear fd

figure(1)
plot(x,c,x,c0,'o'),title('Consumption'), xlabel('capital stock')
% legend(' Approximate',' True if omega=theta=.5'), 
% Because Matlab 4 does not recognize "legend"

figure(2) % Creates a second figure
plot(x,v,x,v0,'o'),title('Value Function'), xlabel('capital stock')
%  legend(' Approximate',' True if omega=theta=.5') 

