m2str(bp2.community.N),'$'],...
    'FontSize',16, 'interpreter','latex');
xlabel('Optimizing by isolated components','FontSize', 20);
set(gcf,'position', [70, 223, 1224, 673]);
%%
% We can see that the default optimization gives a higher modularity value
% (we are optimizing at the entire matrix) but lower number of modules than
% the component optimization. By optimizing at the component level, we do
% not have any more the constraint of the entire network optimization and
% therefore it is possi