sir_simulation_test

sir_simulation_test, an Octave code which calls sir_simulation(), which simulates the spread of a disease through a hospital room of M by N beds, using the SIR (Susceptible/Infected/Recovered) model.

Licensing:

The computer code and data files described and made available on this web page are distributed under the MIT license

Related Data and Programs:

sir_simulation, an Octave code which simulates the spread of a disease through a hospital room of M by N beds, using the Susceptible/Infected/Recovered (SIR) model.

Source Code:

sir_simulation_test.m, calls all the tests.
sir_simulation_test.sh, runs all the tests.
sir_simulation_test.txt, the output file.

As the SIR_SIMULATION program runs, it displays a plot of the hospital room, with susceptibles in green, infecteds in shades of red, and recovereds in gray. Here is an example:

sir_day13.png, a plot of the patient status array on day 13, with M = N = 10, a single infected patient at A(5,5), K = 4, TAU = 0.2, and T_MAX = 50,

Using the command

      plot ( 1:t_max, sir(1,:), 'g', 1:t_max, sir(2,:), 'r', 1:t_max, sir(3,:), 'k' )

we can make line plots of the S, I and R populations. For M = N = 10, a single infected patient at A(5,5), K = 4, TAU = 0.2, and T_MAX = 50, we did this several times:

Using the commands

      h = areas ( sir' );
      set ( h(1), 'FaceColor', [ 0, 1, 0 ] );
      set ( h(2), 'FaceColor', [ 1, 0, 0 ] );
      set ( h(3), 'FaceColor', [ 0.8, 0.8, 0.8 ] );

we can make an area plot of the S, I and R populations. For M = N = 10, a single infected patient at A(5,5), K = 4, TAU = 0.2, and T_MAX = 50, we did this:

sir_area.png

Last revised on 17 March 2019.