% This is the M-Script to run corona_modified_SRI.slx
% Simulation time is T days
T = 200;
N = 2700000; % Total Population is assumed to be constant, no death, no birth and no migration
Io = 96; % Initial number of people infected on 14/03/2020
Ro = 3; % Initial number of prople recovered on 14/03/2020
% OBTAINED by SEARCHING METHOD: 
alpha = 0.006767420000000; % Death Rate = 0.68% Infected
beta = 0.110091900000000; % Every 1000 persons infected 110 persons each day
gamma = 0.017003500000000; % 1.7% infected people recovered each day
sim('corona_modified_SRI.slx');
plot(time, susceptible, 'y', time, infected, 'r', time, recovered, 'g', time, death, '*'), grid on, xlabel('days after 14/03/2020'), ylabel('S,R,I') 