% delta

close all;

msize = 10;
psize = 18;
ColorOrders = [0    0.4470    0.7410
    0.8500    0.3250    0.0980
    0.8790    0.6640    0
    0.4940    0.1840    0.5560
    0.4660    0.6740    0.1880
    0.3010    0.7450    0.9330
    0.6350    0.0780    0.1840];
colororder(ColorOrders);

Zndata = dlmread('../data/pcZn.dat', '', 1, 0);
Dndata = dlmread('../data/pcDn.dat', '', 1, 0);
Endata = dlmread('../data/pcEn.dat', '', 1, 0);

plot(Zndata(:,1), Zndata(:, 4), 'ks-', 'markersize', msize, 'color', ColorOrders(1,:), 'linewidth', 1);
hold on;
plot(Dndata(:,1), Dndata(:, 4), 'kd-', 'markersize', msize, 'color', ColorOrders(2,:), 'linewidth', 1);
plot(Endata(:,1), Endata(:, 4), 'k^-', 'markersize', msize, 'color', ColorOrders(4,:), 'linewidth', 1);

plot(Zndata(:,1), Zndata(:, 6), 'k+--', 'markersize', msize, 'color', ColorOrders(1,:), 'linewidth', 1.25);
plot(Dndata(:,1), Dndata(:, 6), 'k*--', 'markersize', msize, 'color', ColorOrders(2,:), 'linewidth', 1);
plot(Endata(:,1), Endata(:, 6), 'kx--', 'markersize', msize+2, 'color', ColorOrders(4,:), 'linewidth', 1.25);

plot([2, 14], [1, 1], 'color', [0.5,0.5,0.5]);

legend({'$p_\mathrm{c}^\mathrm{site}({\bf Z}^n)$', '$p_\mathrm{c}^\mathrm{site}(D_n)$', '$p_\mathrm{c}^\mathrm{site}(E_n,\Lambda_9)$', '$p_\mathrm{c}^\mathrm{bond}({\bf Z}^n)$', '$p_\mathrm{c}^\mathrm{bond}(D_n)$', '$p_\mathrm{c}^\mathrm{bond}(E_n,\Lambda_9)$'}, 'interpreter', 'latex', 'location', 'southeast', 'numcolumns', 2);
legend boxoff;

xlim([1.9, 13.1]);
ylim([0.35, 1.02]);
xticks(3:3:12);
yticks([0.4:0.2:1]);
xlabel('$n$', 'interpreter', 'latex');
ylabel('$\delta$', 'interpreter', 'latex');
set(gca, 'fontname', 'times new roman', 'fontsize', psize);
print('../figures/Figure2.eps', '-depsc');