% This is the M-Script to run corona_modified_SRI.slx
% Simulation time is Takhir days
Takhir = 100;
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
To = 96 + 3; % Initial PCR-TCM-Tested number
% OBTAINED by SEARCHING METHOD: 
alpha = 0.005989380000000; % Death Rate = 0.56% Infected
beta = 0.152280899999996; % Every 1000 persons infected 109 persons each day
gamma = 0.023577000000000; % 1.8% infected people recovered each day
delta = 0.106677111000000;
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 14/03/20'), ylabel('INFECTED/TESTED [%]')
Final_DEATH = round(death(J))
Final_TESTED = round(tested(J))
Max_INFECTED = round(max(infected))