% This is the M-Script to run corona_modified_SRI.slx
% Simulation time is T days
T = 365;
N = 0.01*270000000; % 1% of the Total Population is assumed to be "related" to the pandemi in one way or the other
Io = 96; % Initial number of people infected on 14/03/2020
Ro = 3; % Initial number of prople recovered on 14/03/2020
Do = 0; % Initial DEATH number
% OBTAINED by SEARCHING METHOD: 
alpha = 0.005660580000000; % Death Rate = 0.56% Infected
beta = 0.108552899999996; % Every 1000 persons infected 109 persons each day
gamma = 0.018073000000000; % 1.8% infected people recovered each day
sim('corona_modified_SRI_revised.slx');
plot(time, susceptible, 'y', time, infected, 'r', time, recovered, 'g',time,death,'*'), grid on, xlabel('days after 20 Mar 2020'), ylabel('S(yellow),R(green),I(red),D(black)')
J = length(death);
Final_DEATH = round(death(J))
Max_INFECTED = round(max(infected))