% This is the M-Script to run corona_modified_SRI.slx
% Simulation time is Takhir days
Takhir = 100;
N = 0.01*19000000; % 1% of the Total Population is assumed to be "related" to the pandemi in one way or the other
N = 0.01*19000000; % Assume: 1% of the total population is affected by the virus
Io = 1; % Initial number of people infected on 14/03/2020
Ro = 0; % Initial number of prople recovered on 14/03/2020
Do = 1; % Initial number of death on 14/03/2020
To = Io + Ro + Do;  % Initial number of tested on 14/03/2020% OBTAINED by SEARCHING METHOD: 
alpha = 0.006347080000000; % Death Rate = 0.56% Infected
beta = 0.229505899999997; % Every 1000 persons infected 109 persons each day
gamma = 0.044089000000000; % 1.8% infected people recovered each day
delta = 0.127692011000000;
sim('corona_modified_SRI_delta_mod.slx');
plot(time, tested,'o', time, infected, 'r', time, recovered, 'g',time,death,'*'), grid on, xlabel('days after 14 Mar 2020'), ylabel('T(blue),R(green),I(red),D(black)')
J = length(death);
for I = 1:J
    PERSEN(I) = (infected(I)/tested(I))*100;
end
pause
clf
plot(time,PERSEN), grid on, xlabel('days after 20/03/20'), ylabel('INFECTED/TESTED [%]')
Final_DEATH = round(death(J))
Final_TESTED = round(tested(J))
Max_INFECTED = round(max(infected))