fp = fopen('scorr.out', 'r');
delta = fscanf(fp, '%e', 1);
[scorr numscorr] = fscanf(fp, '%e');
tau = (0:(numscorr-1)) * delta;

szeroval = scorr(1) / delta;
s = sprintf('<\\sigma_{xy}^2>V/kT = %.3e Pa', szeroval);
fprintf(1, '%s\n', s);

clf;
plot (tau, scorr/scorr(1));
xlabel ('\tau [s]');
ylabel ('<\sigma_{xy}(0)\sigma_{xy}(\tau)>/<\sigma_{xy}^2>');
title (s);
v = axis;
axis([0 max(tau)/10 0 1]);
CHGY;
CHGX;

CPRINT;
