x = linspace(m-3*s,m+3*s).';
g = 1/(s*sqrt(2*pi))*exp(-(x-m).^2/(2*s^2));
histogram(Data,20,'Normalization','pdf','FaceColor', 'c');
plot(x,g,'c', 'LineWidth', 2);
xregion(1, m+3*s,'FaceColor', [1 0 0])
histogram(Data,20,'Normalization','pdf','FaceColor', 'c');
plot(x,g,'c', 'LineWidth', 2);
xregion(1, m+3*s,'FaceColor', [1 0 0],'FaceAlpha',1)
title('ConstantRegion FaceAlpha=1')
histogram(Data,20,'Normalization','pdf','FaceColor', 'c');
plot(x,g,'c', 'LineWidth', 2);
xregion(1, m+3*s,'FaceColor', [1 0 0])
title('Default with black axes')
histogram(Data,20,'Normalization','pdf','FaceColor', 'c');
plot(x,g,'c', 'LineWidth', 2);
xregion(1, m+3*s,'FaceColor', [1 0 0],'FaceALpha',1)
title('ConstantRegion FaceAlpha=1')
histogram(Data,20,'Normalization','pdf','FaceColor', 'c','FaceAlpha',1);
plot(x,g,'c', 'LineWidth', 2);
xregion(1, m+3*s,'FaceColor', [1 0 0])
title('Histogram FaceAlpha=1')
histogram(Data,20,'Normalization','pdf','FaceColor', 'c','FaceAlpha',1);
plot(x,g,'c', 'LineWidth', 2);
xregion(1, m+3*s,'FaceColor', [1 0 0],'FaceAlpha',1)
title('Histogram & ConstantRgn FaceAlpha=1')