clear all
beta = 0.1;
%for k = 1:2
for k = 1:221
i = k - 11 ;
  tk = i / 100.0 ;
  t(k) = tk ;
  y(k) = 0.1 * (tk-1.0)^2 ;
z(k) = 0.1 - 0.2 * tk;
a(k) = 0.1 - 0.2 * beta * tk;
end

plot(t,y)
, axis([-0.1,2.0,-0.05,0.14])
hold on
plot(t,z,'--g')
, axis([-0.1,2.0,-0.05,0.14])
plot(t,a,'--r')
, axis([-0.1,2.0,-0.05,0.14])

alphamin = [ 0 0 ] ;
vertical = [ -0.05 ,0.14 ];
plot( alphamin, vertical, ':k' );

alphamax = [ 1.8 1.8 ] ;
vertical = [ -0.05 ,0.14 ];
plot( alphamax, vertical, ':k' );


hold off

print -depsc eps/l23
