plotode.m

global beta delta gamma N omega q;

%str=['beta=', num2str(beta); 'delta=',  num2str(delta) ; 'gamma=', num2str(gamma)];
%newstr=join(str);

load 'Covid19infecAUT.mat'
covid19infection= table2array(covid19(:,2));
covid19recovered= table2array(covid19(:,4));
covid19deceased= table2array(covid19(:,3));
covid19active= table2array(covid19(:,5));
y=[1:size(covid19infection)]';


options = odeset('RelTol', 1e-4, 'NonNegative', [1 2]);
[t,x] = ode45('SEIRmodel', [0:1:size(covid19infection)], [17300 1 1 0 0]);

figure(1)
clf
plot(t,x,'-');
hold on
plot(y,covid19infection,'o');
hold on
plot(y,covid19recovered,'o');
hold on
plot(y,covid19deceased,'o');
hold on
plot(y,covid19active,'o');
legend('[S]', '[E]', '[I]', '[R]','[D]','Total Infections','Recovered','Deceased','Active','FontSize');
title('Covid-19 SEIR model Austria 26.02-now')
dim = [.2 .5 .5 .3];
%annotation('textbox',dim,'String',newstr ,'FitBoxToText','on');

xlabel('Days since 26.02.2020, t','FontSize', 18);
ylabel('Number of people','FontSize', 18);
ylim([0.0 20000]);