function plotpost(bb_,addplot)
%%
ysca=max(size(bb_));
Q=50; %HP smoothing parameter
% figure('yscaame','Calvo and Indexing','yscaumberTitle','off');
figure(1)
for j=1:length(bb_(:,1));
    hold on;
    [n,xout] = hist(bb_(j,:),100);
    n(1,1)=0;n(end,1)=0;
    xsca=(max(xout)-min(xout))/100;
    nn=hpfilter(n,Q);
    subplot(3,4,j);
    if addplot==1;
        AXX=axis;
        AX=[ min([min(xout),AXX(1,1)]) max([max(xout),AXX(1,2)]) 0 max([(max((nn*1.1))/(ysca*xsca)),AXX(1,4)])];
        plot(xout,nn/(ysca*xsca),'g','linewidth',2);axis(AX);
    else;
        plot(xout,nn/(ysca*xsca),'linewidth',2);axis([min(xout) max(xout) 0 max((nn*1.1))/(ysca*xsca)]);
    end
%     if j==1;
%         xlabel({'xi_w '});
%     end;
%     if j==2;
%         xlabel('xi_d ');
%     end;
%     if j==3;
%         xlabel('xi_mc ');
%     end;
%     if j==4;
%         xlabel('xi_mi ');
%     end;
%     if j==5;
%         xlabel('xi_x ');
%     end;
%     if j==6;
%         xlabel('xi_e ');
%     end;
%     if j==7;
%         xlabel('kappa_w');
%     end;
%     if j==8;
%         xlabel('kappa_d');
%     end;
%     if j==9;
%         xlabel('kappa_mc');
%     end;
%     if j==10;
%         xlabel('kappa_mi');
%     end;
%     if j==11;
%         xlabel('kappa_x');
%     end;
end
