
function demo_covid

% load covid.dat;
covid=dati_covid;

% grado polinomio minimi quadrati
grado=5;

% -------------- PARTE INFETTI --------------

% vettore colonna "nuovi_infetti" di nuovi infetti, in cui la prima 
% componente e' la piu' recente
nuovi_infetti=covid(:,1);

% ricavo vettore colonna "nuovi_infetti" di nuovi infetti, in cui la prima 
% componente e' la meno recente
nuovi_infetti=flipud(nuovi_infetti);

% elimina contagio pari a 0
i=0; contagio0=1;
while contagio0 == 1
    i=i+1;
    if nuovi_infetti(i) > 5
        contagio0=0;
    end
end

nuovi_infetti_purged=nuovi_infetti(i:end);
totale_giorni=length(nuovi_infetti_purged);
giorni=(1:totale_giorni)';

% minimi quadrati
p=polyfit(giorni,nuovi_infetti_purged,grado);
pvals=polyval(p,giorni);

% grafico dati e minimi quadrati
figure(1);
plot(giorni,pvals,'r-',giorni,nuovi_infetti_purged,'ko');

% previsione tra 1 giorno
previsione=max([polyval(p,totale_giorni+1) 0]);
fprintf('\n \t previsione infetti tra 1 giorno: %5.0f \n',previsione);






% -------------- PARTE DECEDUTI --------------

% vettore colonna "nuovi_deceduti" di nuovi deceduti, in cui la prima 
% componente e' la piu' recente
nuovi_deceduti=covid(:,2);

% ricavo vettore colonna "nuovi_deceduti" di nuovi deceduti, in cui la prima 
% componente e' la meno recente
nuovi_deceduti=flipud(nuovi_deceduti);

% elimina contagio pari a 0
i=0; contagio0=1;
while contagio0
    i=i+1;
    if nuovi_deceduti(i) > 6
        contagio0=0;
    end
end

nuovi_deceduti_purged=nuovi_deceduti(i:end);
totale_giorni=length(nuovi_deceduti_purged);
giorni=(1:totale_giorni)';

% minimi quadrati
p=polyfit(giorni,nuovi_deceduti_purged,grado);
pvals=polyval(p,giorni);

% grafico dati e minimi quadrati
figure(2);
plot(giorni,pvals,'r-',giorni,nuovi_deceduti_purged,'ko');

% previsione tra 1 giorno
previsione=max([polyval(p,totale_giorni+1) 0]);
fprintf('\n \t previsione deceduti tra 1 giorno: %5.0f \n',previsione);


fprintf('\n \n');


%------------------------------------------------------------
% >> demo_covid
% 
% 	 previsione infetti tra 1 giorno:  2102 
%
% 	 previsione deceduti tra 1 giorno:   154 
%  
% >>
%------------------------------------------------------------