% This is the M-Script to run corona.slx
% Simulation time is T days
T = 100;
N = 1000; % Total Population is assumed to be constant, no death, no birth and no migration
Io = 2; % Initial number of people infected
So = 500; % 50 percents of the people are initially susceptible
beta = 1.5; % Every 2 people infected 3 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') 