Text this: Lyapunov Functions for a Class of Discrete SIRS Epidemic Models with Nonlinear Incidence Rate and Varying Population Sizes