% Displays the output of the wator program in two plots against time:
% RMS velocity (related to temperature), and time step size.

% Read in output
load graph.dat;
time = graph(:, 1);
mnsqvel = graph(:, 2);

%  plot mean square velocities versus time
      steps = length(time);
	  t = time(steps);
      subplot(121), 
	  axis([0,t,0,max(mnsqvel)]);
      plot(time,mnsqvel); title('Mean square velocity'); xlabel('time');
%  plot timesteps versus time
      subplot(122),
      timestep = time(2:steps)-time(1:steps-1);
      axis([0,steps-2,0,max(timestep)]);
      plot(timestep); title('time step'); xlabel('step number');