close all;
y = randn(500,1);
% four groups
groups = randi(4,500,1);
% number of bins
nbins = 10;
[ng, hb] = histFS(y,nbins,groups);
title('Histogram with default color rotation (brcmykg)','interpreter','latex','FontSize',18);

figure;
y = randn(500,1);
groups = randi(4,500,1);
nbins = 10;
[ng, hb] = histFS(y,nbins,groups,{'BMW','FIAT','VOLVO','FERRARI'});
title('The same, with personalized legends','interpreter','latex','FontSize',18);

figure;
y = randn(500,1);
groups = randi(4,500,1);
nbins = 10;
[ng, hb] = histFS(y,nbins,groups,{'BMW','FIAT','VOLVO','FERRARI'},gca,'kgbr');
title('FERRARI must be red! Color sequence changed to kgbr','interpreter','latex','FontSize',18);

close all;
y = randn(500,1);
groups = randi(4,500,1);
nbins = 10;
hs = gscatter(1:numel(y),y,groups);
hfs = gcf;                              % get the handle of the scatterplot
has = get(hfs,'CurrentAxes');           % it is the same as has = gca
hlegends  = get(has,'Children');        % get the legend entries
getleg = get(hlegends,'DisplayName');   % get the names of the legend entries
getcol = get(hlegends,'Color');         % get the color of the legend entries
getcolm = cell2mat(getcol);             % arrange the RGB vectors into a matrix
figure;
[ng, hb] = histFS(y,nbins,groups,getleg,gca,getcolm);
title('Color sequence is taken from the scatterplot','interpreter','latex','FontSize',18);

y = randn(500,1);
groups = randi(4,500,1);
nbins = 10;
myleg = {'my group 1' 'my group 2' 'my group 3' 'my group 4' };
[ng, hb] = histFS(y,nbins,groups,myleg);
title('Example with personalised clickable legends','interpreter','latex','FontSize',18);

close all;clear all;
X = crosstab2datamatrix([10 20 30]); X=X(:,2);
X1 = 444*ones(60,1); X1(40:60)= 999; X1 = shuffling(X1);
X = X+5;
if verLessThan('matlab','8.5')
groups = num2str(X1);
else
groups = categorical(X1);
end
W = ones(size(X,1),1); W(12) = 100;
histFS(X,10,groups,[],[],[],W);
title({'Weighted histograms:' , 'please check with tabulateFS(W) and tabulateFS(X)'},'interpreter','latex','FontSize',18);
figure;
histFS(X,10,groups);
title('Standard histogram','interpreter','latex','FontSize',18);
cascade;
disp('tabulateFS(X)');
tabulateFS(X);
disp('tabulateFS(W)');
tabulateFS(W);