function PlotFig10Densest

close all;
load densityconvergence.dat

loglog(densityconvergence(:, 1), densityconvergence(:, 2), 'o');
coeff = polyfit(log(densityconvergence(:, 1)), log(densityconvergence(:, 2)), 1);
x = [0 13.8155];
y = polyval(coeff, x);
hold on;
loglog(exp(x), exp(y), 'k--');
xlim([1 1e6]);
ylim([1e-6 5]);
set(gca, 'FontSize', 16);
xlabel('n');
ylabel('\rho_n-\rho_\infty');
saveas(gca, 'Fig10.png');