task_syn    =  [0.0274029;  0.0707994;  0.132135; 0.176905; 0.0320425];
task_full   =  [0.00700702; 0.00428036; 0.0254512; 0.0598569; 0.00436699];
task_hybrid =  [0.00436661; 0.0158798;  0.0744908; 0.0928952; 0.00233551];

act_syn =    [0.036104;  0.048403;  0.0467951; 0.0583683; 0.0447815];
act_full =   [0.124487;  0.10533;   0.0973407; 0.0814257; 0.0836098];
act_hybrid = [0.0695954; 0.063004;  0.0713405; 0.073988; 0.0689496];

total_syn =    [0.190716;  0.127497;  0.185994; 0.249279; 0.0982359];
total_full =   [0.186303;  0.11946;   0.164708; 0.179036; 0.113494];
total_hybrid = [0.0907391; 0.0913529; 0.157159; 0.17843; 0.0972878];

y = [mean(task_syn), mean(act_syn), mean(total_syn - task_syn - act_syn);
     mean(task_full), mean(act_full), mean(total_full  - task_full - act_full);  
     mean(task_hybrid), mean(act_hybrid), mean(total_hybrid - task_hybrid - act_hybrid);
];

width = 3;     
height = 3;    
alw = 2.75;    
fsz = 24;      
lw = 2.5;      
msz = 8;       

bar(y,'stacked');
legend('Task', 'Activation', 'Other');
ylabel('Cost Function Value');
set(gca, 'xticklabels', {'12 Syn', 'Full', 'Hybrid', }, 'FontSize', fsz, 'LineWidth', lw);
legend boxoff;
box off;

% Other Statistics in Section 3.4
display('Synergies vs Full');
other_statistics(task_syn, act_syn, task_full, act_full);
display('Hybrid vs Full');
other_statistics(task_hybrid, act_hybrid, task_full, act_full);