xErr = randn(1,1000);
[mu, sig] = normfit(xErr);
h = histogram(xErr, 100, 'Normalization','pdf');
set(h_xErr, {'DisplayName'}, {['Standard deviation $sigma_{x} = $ ', num2str(sigX)]; ['Mean $mu_x = $ ', num2str(muX)]});
hl = legend('Location', 'NorthWest');

h = [h; histogram(xErr, 100, 'Normalization','pdf')];