% 
% Written by:
% -- 
% John L. Weatherwax                2006-05-17
% 
% email: wax@alum.mit.edu
% 
% Please send comments and especially bug reports to the
% above email address.
% 
%-----

close all; 

% Get statistics on the size of the domain required to plot:
% 
files = dir( 'sequential_ts_*.dat' ); 
% Sort the files by index

if( 0 ) 
x_max = -Inf; x_min = +Inf; 
y_max = -Inf; y_min = +Inf; 
for fi=1:length(files),
  X = load( files(fi).name ); 
  x_min = min([X(:,1);x_min]); x_max = max([X(:,1);x_max]); 
  y_min = min([X(:,2);y_min]); y_max = max([X(:,2);y_max]); 
end
end

figure; carray = 'rgbcmyk';
for fi=1:length(files),
  X = load( files(fi).name ); 
  %quiver( X(:,1), X(:,2), X(:,3), X(:,4) ); grid on; hold on; 
  h(fi)=plot( X(:,1), X(:,2), ['.',carray(mod(fi,length(carray))+1)] ); grid on; hold on; 
  %axis( [ 0.9*x_min, 1.1*x_max, 0.9*y_min, 1.1*y_max ] ); 
  axis( [ -1 +1 -1 +1 ] );
  title( ['File: ',files(fi).name], 'interpreter', 'none' ); 
  if( fi > 5 ) 
    delete(h(fi-5));
  end
  drawnow; 
  pause; 
end