% This is the M-Script to run corona.slx
% Simulation time is T days
T = 53;
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
So = 0.5*2700000; % 50 percents of the people are initially susceptible
beta = 0.39; % Every 1 people infected 1 persons
gamma = 0.1; % 10% of infected people recovered each day
sim('corona.slx');
plot(time, susceptible, 'y', time, infected, 'r', time, recovered, 'g'), grid on, xlabel('days'), ylabel('S,R,I') 