dS = (-beta * I ) * S
dE = (beta * I) * S - (a * E)
dI = (a * E) - (g * I)
dZ = pZ * g * I - gZ * Z
dH = pH * g * I - gH * H
dC = pH * pC * g * I - gC * ICU
dDs = pfD * pH * pC * g * I - gCD * Ds
dDa = pA * g * I - gA * Da
dRZ = gZ * Z
dRH = gH * H
dRC = gC * C
dRDs = gCD * Ds
dRDa = gA * Da