solving ordinary differential equations programming the numerical methodthe special case the computation of ( )-( can be readily made in computer program sir pyfrom numpy import zeroslinspace import matplotlib pyplot as plt time unit beta /( * * gamma /( * dt min simulate for days n_t int( * /dtcorresponding no of hours linspace( n_t*dtn_t+ zeros(n_t+ zeros(n_t+ zeros(n_t+ initial condition [ [ [ step equations forward in time for in range(n_t) [ + [ndt*beta* [ ]* [ni[ + [ndt*beta* [ ]* [ndt*gamma* [nr[ + [ndt*gamma* [nfig plt figure( plt plot(tstitrfig legend(( )(' '' '' ')'upper left'plt xlabel('hours'plt show(plt savefig('tmp pdf')plt savefig('tmp png'this program was written to investigate the spreading of flu at the mentioned boarding schooland the reasoning for the specific choices and goes as follows at some other school where the disease has already spreadit was observed that in the beginning of day there were susceptibles and infectedwhile the numbers were and respectively hours later using as time unitwe then have from ( that among infectedit was observed that recovered during daygiving applying these parameters to new case where there is one infected initially and susceptiblesgives the graphs in fig these graphs are just straight lines between the values at times ti it as computed by the program we observe that reduces as and grows after about days everyone has become ill and recovered again we can experiment with and to see whether we get an outbreak of the disease or not imagine that "wash your handscampaign was successful and that the other school in this case experienced reduction of by factor of with this