function PlotGz()

close all;
load gr.dat
load peak.dat

semilogy(gr(:, 1), abs(gr(:, 2) - 1.0));
xlim([0 15]);
ylim([1e-3 1]);
hold on;
semilogy(peak(:, 1), peak(:, 2), 'go', 'MarkerSize', 7, 'LineWidth', 1);
coeff = polyfit(peak(:, 1), log(peak(:, 2)), 1);
x = [0 15];
y = polyval(coeff, x);
plot(x, exp(y), 'k', 'LineWidth', 1);
set(gca, 'FontSize', 16);
xlabel('z/\sigma');
ylabel('|g(z)-1|');
axes('Position', [0.56 0.61 0.32 0.29]);
box on;
plot(gr(:, 1), gr(:, 2));
xlim([0 15]);
ylim([0.5 2]);
set(gca, 'FontSize', 16);
xlabel('z/\sigma');
ylabel('g(z)');
saveas(gca, 'gzSemilog.png');