function run_large_scale_org(hemi)

back = pwd;

%=====================
% DEFAULTS
%=====================
% Use different values of seed and target branches to explore affinity
% matrices
seed_branches   = 8;
target_branches = 7;

affinity_matrix = make_affinity_matrix(seed_branches, target_branches, hemi);

if upper(hemi)=='R'
    plot_affinity_matrix(affinity_matrix,'R frontal to R parietal: affinity matrix');
elseif upper(hemi)=='L'
    plot_affinity_matrix(affinity_matrix,'L frontal to L parietal: affinity matrix');
end

cd(back);

cd(['results/largescale_organization/' hemi '_hemisphere/'])
savefig affinity_matrix

cd(back);

end