clc;
clear;

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

% creo array acc e forze
acc = dati(:,1);
forze = dati(:,2);

% controllo che non ci siano dati errati
i = 1;
while ((acc(i)>=0 || forze(i)>=0) && i<length(acc))
    if acc(i)<0 ||forze(i)<0
        disp('errore nei dati')
    end
    i = i+1;
end
 
% plot dati sperimentali
plot(acc,forze,'r.'),xlabel('acc (m/s^2)'),ylabel('F(Newton)');
hold on

% calcolo vettore delle masse e la sua media
masse = forze./acc;
media_masse = mean(masse);

%calcolo numero di valori maggiori della media
num = 0;
for i = 1:1:length(masse)
    if masse(i) > media_masse
        num = num + 1;
    end
end
disp('numero di medie maggiori della media: ')
disp(num);
    
% plot legge di Newton
acc_teo = [min(acc):0.1:max(acc)];
f_teo=media_masse .* acc_teo;
plot(acc_teo,f_teo);