%% Australia covid cases

data = readtable("Australia_covid_cases.xlsx")

t = data.Date_reported(:);
cases = data.Cumulative_cases(:);
figure
plot(t,cases,'-','DisplayName','Cumulative cases')
legend, grid 

%

x = cases(1:end-10);
sys1 = ar(x,2)
compare(cases,sys1,10)