%  mfile plotoutput
%  plot mean square velocities versus time
      steps = length(time);
      subplot(121), 
      axis([0,t,0,max(mnsqvel)]); 
      plot(time,mnsqvel,'k'); 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,'k'); title('time step'); xlabel('step number');