def sir_conserved ( t, y ):

#*****************************************************************************80
#
## sir_conserved() returns a conserved quantity for the SIR ODE.
#
  import numpy as np
#
#  Each entry of H is the sum of a row (axis=1) of Y values.
#
  h = np.sum ( y, axis = 1 )

  return h