clc;
clear;

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

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

% controllo che non ci siano dati errati
i = 1;
while ((volt(i)>=0 || corr(i)>=0) && i<length(corr))
    if volt(i)<0 ||corr(i)<0
        disp('errore di input')
    end
    i = i+1;
end
 
% plot dati sperimentali
plot(corr,volt,'r.'),xlabel('I(A)'),ylabel('volt(V)');
hold on

% calcolo vettore delle resistenze e la sua media
r = volt./corr;
media = mean(r);

% plot legge di Ohm
corr_teo = [min(corr):0.1:max(corr)];
volt_teo=media .* corr_teo;
plot(corr_teo,volt_teo);

%calcolo numero di valori maggiori della media
scarto = sum(r-media)^2;
disp('scarto: ')
disp(scarto);