% This is the M-Script to run corona_modified_SRI.slx
% Simulation time is T days
% SOUTH SULAWESI CASE
T = 365;
N = 0.01*19000000; % 1% of the population
Io = 2; % Initial number of people infected on 20/03/2020
Ro = 0; % Initial number of prople recovered on 20/03/2020
Do = 1; % Initial Death
% OBTAINED by SEARCHING METHOD: 
alpha = 0.006409420000000; % Death Rate = 0.64% Infected
beta = 0.154031899999999; % Every 1000 persons infected 148 persons each day
gamma = 0.038943500000000; % 3.9% 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))
