We describe and analyze a periodically-forced difference equation model for malaria in mosquitoes that captures the effects of seasonality and allows the mosquitoes to feed on a heterogeneous population of hosts. We numerically show the existence of a unique globally asymptotically stable periodic orbit and calculate periodic orbits of field-measurable quantities that measure malaria transmission. We integrate this model with an individual-based stochastic simulation model for malaria in humans to compare the effects of insecticide-treated nets (ITNs) and indoor residual spraying (IRS) in reducing malaria transmission, prevalence, and incidence. We show that ITNs are more effective than IRS in reducing transmission and prevalence though IRS would achieve its maximal effects within 2 years while ITNs would need two mass distribution campaigns over several years to do so. Furthermore, the combination of both interventions is more effective than either intervention alone. However, although these interventions reduce transmission and prevalence, they can lead to increased clinical malaria; and all three malaria indicators return to preintervention levels within 3 years after the interventions are withdraw
