clf;
load gr.out;
plot (gr(:,1), gr(:,2));
xlabel('r [Angstrom]');
ylabel('g(r)');
CHGX;

if exist('T'), 
 a = input(['T = ' T ' K? '], 's');
 if ~isempty(a), T = a; end;
else
 T = input('T = ', 's');
end;

if exist('rho'), 
 a = input(['density = ' rho ' mol/m^3? '], 's');
 if ~isempty(a), rho = a; end;
else
  rho = input('density(mol/m^3) = ', 's');
end;

title (['RDF of Argon at T = ' T 'K, \rho = ' rho ' mol/m^3']); 

CPRINT;

