function print_exploration_subject(t,relative1,relative2,banda)
%ch_names ={'F left';'F right';'C left';'C right'; 'P left'; 'P right'};
lgd = {'Familiar ACQ';...
           'Familiar STM';...
           'Familiar LTM';...
           'Familiar ACQ';...
           'Novel STM';...
           'Novel LTM';};

% colors1=[250 18 18;...
%                233 13 33;... %220 54 54;...
%                134 16 16]/256;
%                 
% 
% colors2 = [54 111 200;...
%                    13 69 155;...
%                    38 68 115]/256;

colors1=[233 162 13;...
               233 86 13;... %220 54 54;...
               233 33 13]/256;
                

colors2 = [13 203 233;...
                   13 74 233;...
                   114 13 233]/256;


f=figure;
f.Position = [5 1 1859 949];

[n_epochs1, n_muestras, n_channels] = size(relative1);
[n_epochs2, n_muestras, n_channels] = size(relative2);



for i=1:6    
%         for r=1:3
%             ax(i) = subplot(3,2,i);
%             media1 = squeeze(relative1(r,:,i));
%             media2 = squeeze(relative2(r,:,i));
%             plot(t,media1,'Color',colors1(r,:),'Linewidth',4);
%             hold on
%             plot(t,media2,'Color',colors2(r,:),'Linewidth',4);
%         end

        ax(i) = subplot(3,2,i);

        hold on
        for r=1:3
            media1 = squeeze(relative1(r,:,i));
            plot(t,media1,'Color',colors1(r,:),'Linewidth',4-r+1);
        end

        for r=1:3
            media2 = squeeze(relative2(r,:,i));
            plot(t,media2,'Color',colors2(r,:),'Linewidth',4-r+1);
        end

        hold off
        
            axis([-500 3000, -70,-0])
            %xlim([-500 3000])
            grid on
        
        
%             ax(i) = subplot(3,2,i);
%             plot(t,squeeze(relative1(:,:,i)),'Color',colors1,'Linewidth',4);
%             hold on
%             plot(t,squeeze(relative2(:,:,i)),'Color',colors2,'Linewidth',4);

            if i==1
                legend(lgd,'Location','Best')
            end

            p = get(gca, 'Position');
            %h = axes('Parent', gcf, 'Position', [p(1)+1.07*p(3) p(2)+p(4)*.55 0.15*p(3) 0.3*p(3)]);
            h = axes('Parent', gcf, 'Position', [p(1)+1.07*p(3) p(2)+p(4)*.15 0.15*p(3) 0.5*p(3)]);
            box on
            boxplot(h,[squeeze(relative1(:,:,i))',squeeze(relative2(:,:,i))'], 'boxstyle','filled','labels',lgd,'colors',[colors1;colors2])
            xticks([])
            hold off
    
end

% for i = 1:2
%         %subplot(3,2,i)
%         title (ax(i),ch_names{i});
%         grid on;
% end

subplot(3,2,1)
title (ax(1),'Left')

subplot(3,2,2)
title (ax(2),'Right')

canal={'Frontal';'Central';'Parietal'};

j=1;
for i = 1:2:6
    %subplot(3,2,i)
    
    text(ax(i),3500,10,canal{j},'HorizontalAlignment','Center', 'FontSize',14,'FontWeight','bold')
    j = j+1;
    
    if banda==1
        ylabel(ax(i),'\alpha relative power (AdB)')
    else 
        ylabel(ax(i),'\beta relative power (AdB)')
    end
end

for i = 5:6
    %subplot(3,2,i)
    xlabel(ax(i),'time (ms)')
end
    
