Untitled

function [var_perc, is_virtuoso, is_canaglia] = analisi(E20, E21)

    var_perc = ((E21 - E20) ./ E20) * 100;

    is_virtuoso = var_perc <= -5;

    is_canaglia = var_perc >= 2;
end

load('emissioni.mat');

[var_perc, is_virtuoso, is_canaglia] = analisi(E20, E21);

disp('Paesi virtuosi: ');
disp(find(is_virtuoso));
disp('Paesi canaglia:');
disp(find(is_canaglia));

var_globale = mean(var_perc);
disp(var_globale);

save('risultati.mat', var_perc, find(is_canaglia));


[max_incr, idx_max_incr] = max(var_perc);
[max_decr, idx_max_decr] = min(var_perc);

figure;
plot(1:length(var_perc), var_perc, 'b-');
hold on;
plot(idx_max_incr, max_incr, 'r*');
plot(idx_max_decr, max_decr, 'yo');
grid on;
xlabel('Paese');
ylabel('Variazione percentuale CO2');
title('Variazioni delle emissioni di CO2 per paese');