clc;
clear;

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

% creo array 
pressione = dati(:,1);
volume = dati(:,2);

% controllo che non ci siano dati errati
i = 1;
while ((pressione(i)>=0 || volume(i)>=0) && i<length(pressione))
    if pressione(i)<0 || volume(i)<0
        disp('errore di input')
    end
    i = i+1;
end
 
figure(1)
% plot dati sperimentali
plot(pressione,volume,'r.'),xlabel('pressione (atm)'),ylabel('volume (l)');
hold on

% plot legge dei gas
pressione_teo = [min(pressione):0.1:max(pressione)];
vol_teo=(R * 250)./pressione_teo;
plot(pressione_teo,vol_teo);
hold off;

figure(2)
%vettore temperature
temp = input('inserisci vettore delle temperature');
temp_K = temp + 273;
vol = [];
for i = 1:1:length(temp)
    vol=(R * temp_K(i))./pressione_teo;
    plot(pressione_teo,vol),xlabel('pressione (atm)'),ylabel('volume (l)');
    hold on;
end