bar(Cvec)
hold on;
plot([0, 9], [mean(Cvec), mean(Cvec)]);
plot([0, 9], [mean(Cvec)+std(Cvec), mean(Cvec)+std(Cvec)], 'm--');
plot([0, 9], [mean(Cvec)+2*std(Cvec), mean(Cvec)+2*std(Cvec)], 'r:', 'LineWidth', 2);
plot([0, 9], [mean(Cvec)-std(Cvec), mean(Cvec)-std(Cvec)], 'm--');
plot([0, 9], [mean(Cvec)-2*std(Cvec), mean(Cvec)-2*std(Cvec)], 'r:', 'LineWidth', 2);
xlabel('Index of data', 'FontSize', 18);
ylabel('Estimation of C', 'FontSize', 18);
legend('C', 'mean', 'mean+sd', 'mean+2sd');
orient landscape;
print('-dpng', 'Csbarchart.png');
hold off;