more off;
load struc.out;
name = {
't', 'ps'
'T', 'K'
'Hamiltonian / atom', 'J/mol'
'structure factor', ''
'Mean squared displacement', 'A^2'
'H_{1x}', 'A'
'H_{2y}', 'A'
'H_{3z}', 'A'
'H_{2z}', 'A'
'H_{1z}', 'A'
'H_{1y}', 'A'
'\sigma_{xx}', 'MPa'
'\sigma_{yy}', 'MPa'
'\sigma_{zz}', 'MPa'
'\sigma_{yz}', 'MPa'
'\sigma_{xz}', 'MPa'
'\sigma_{xy}', 'MPa'
'C_{xxxx}', 'MPa'
'C_{xxyy}', 'MPa'
'C_{xxzz}', 'MPa'
'C_{xxyz}', 'MPa'
'C_{xxxz}', 'MPa'
'C_{xxxy}', 'MPa'
'C_{yyyy}', 'MPa'
'C_{yyzz}', 'MPa'
'C_{yyyz}', 'MPa'
'C_{yyxz}', 'MPa'
'C_{yyxy}', 'MPa'
'C_{zzzz}', 'MPa'
'C_{zzyz}', 'MPa'
'C_{zzxz}', 'MPa'
'C_{zzxy}', 'MPa'
'C_{yzyz}', 'MPa'
'C_{yzxz}', 'MPa'
'C_{yzxy}', 'MPa'
'C_{xzxz}', 'MPa'
'C_{xzxy}', 'MPa'
'C_{xyxy}', 'MPa'
};

fprintf (1, '\n');
n = size (name,1);
for i = 1:n,
  if (size(name{i,1},2) > 21) buf = '\t';
  elseif (size(name{i,1},2) > 13) buf = '\t\t';
  elseif (size(name{i,1},2) > 1) buf = '\t\t\t';
  else buf = '\t\t\t\t';
  end;
  fprintf (1, ['%d. %s' buf], i, name{i,1});
  if (mod(i,2) == 0) fprintf (1, '\n'); end;
end;
fprintf (1, '\n');

c = 1;
d = input (['Select a quantity to plot with ' name{c,1} ...
      ' [' name{c,2} ']: ']);

plot (struc(:,c), struc(:,d));
xlabel([name{c,1} ' [' name{c,2} ']']);
ylabel([name{d,1} ' [' name{d,2} ']']);
hold on;
