load model3f.dat
nkt    = 10;
nki    = 10;
kt     = linspace(1.25,1.75,nkt)';
ki     = linspace(1.2,2,nki)';
ivr    = input('Which endogenous variable? (1=c,2=xt,3=xi,4=h1,5=h2,6=k1t,7=k2t) ');
ist    = input('Which exogenous state? (1,...,ns) ');
c     = zeros(nkt,nki);
xt    = zeros(nkt,nki);
xi    = zeros(nkt,nki);
h1    = zeros(nkt,nki);
h2    = zeros(nkt,nki);
k1t   = zeros(nkt,nki);
k2t   = zeros(nkt,nki);
c(:)  = model3f(1:nkt*nki,ist);
xt(:) = model3f(nkt*nki+1:2*nkt*nki,ist);
xi(:) = model3f(2*nkt*nki+1:3*nkt*nki,ist);
h1(:) = model3f(3*nkt*nki+1:4*nkt*nki,ist);
h2(:) = model3f(4*nkt*nki+1:5*nkt*nki,ist);
k1t(:)= model3f(5*nkt*nki+1:6*nkt*nki,ist);
k2t(:)= model3f(6*nkt*nki+1:7*nkt*nki,ist);
cmod  = model3f(1:nkt*nki,:);
xtmod = model3f(nkt*nki+1:2*nkt*nki,:);
ximod = model3f(2*nkt*nki+1:3*nkt*nki,:);
h1mod = model3f(3*nkt*nki+1:4*nkt*nki,:);
h2mod = model3f(4*nkt*nki+1:5*nkt*nki,:);
k1tmod= model3f(5*nkt*nki+1:6*nkt*nki,:);
k2tmod= model3f(6*nkt*nki+1:7*nkt*nki,:);

if ivr==1
  plot(kt,c)
elseif ivr==2
  plot(kt,xt)
elseif ivr==3
  plot(kt,xi)
elseif ivr==4
  plot(kt,h1)
elseif ivr==5
  plot(kt,h2)
elseif ivr==6
  plot(kt,k1t)
elseif ivr==7
  plot(kt,k2t)
end;
