
NV=10.^(1:0.1:10);
dV=0:10;

cmap = hsv(20);

k=1;
for d=dV
    color = cmap(k,:);
    k=k+1;

    QMC_exp=(log(NV)).^d./NV;
    MC_exp=1./sqrt(NV);

    loglog(NV,QMC_exp,'-','Color',color,'LineWidth',4);
    hold on;
end
legendStrings = "d = " + string(dV);
legend(legendStrings);
loglog(NV,MC_exp,':','Color',color,'LineWidth',4);

hold off;
