%
% Example 2 from Sec 5.6
%
A = [.8 0; 0 .64]

% Invoke plot screen
axis([-100,100,-100,100])
grid
hold on

% Choose an initial vector
x = [input('x-axis = '); input('y-axis = ')]

% Choose the running time
n = input('n = ')

% Graph the numerical solution.
for i = 1:n
   xa = A * x;
   arrow(x,xa-x,'r',2);
   x = xa;
end

