clc;
clear;

% carico in dati tutti i dati del foglio dati_richter
dati = xlsread('dati_richter');
n = size(dati);

% creo array 
magnitudo = dati(:,1);
energia = dati(:,2);
identificativo = dati(:,3);

% controllo input
i = 1;
while (magnitudo(i)<8.6 && i<length(magnitudo))
    if magnitudo(i)>8.6 
        disp('terremoto mai visto!')
    end
    i = i+1;
end

% magnitudo
num = 0;
vett_identificativo=[];
for i = 1:1:length(magnitudo)
    if magnitudo(i) > 7
        num = num + 1;
        vett_identificativo=[vett_identificativo,identificativo(i)];
    end
end
disp('numero terremoti gravi: ');
disp(num);
disp('identificativi');
disp(vett_identificativo);

% plot dati sperimentali
plot(magnitudo,energia,'r.'),xlabel('magnitudo (Richter)'),ylabel('energia (J)');
hold on

% plot legge di Richter
magnitudo_teo = [min(magnitudo):0.1:max(magnitudo)];
energia_teo=10^4.4 * 10 .^ (1.5 *  magnitudo_teo);
plot(magnitudo_teo,energia_teo);