“…The ode45 function evaluates the differential equations using an explicit 4 th order Runge-Kutta method for solving ODEs. In these simulations, initial population sizes are chosen based on populations used in similar studies I product of human transmission probability (.5) and biting rate (.5), from [3] II product of recovery rate (.005) and modified recovery rate (8.04), from [3] III sum of treatment rate (.2) and recovery rate (.005), from [3] IV product of mosquito transmission probability (.83) and biting rate (.5), from [3] V division by 365 allows for the rate to have the correct units of days [4,3,5]. The initial conditions used are S h = 300, E h = 0, I h = 1, R h = 0, S m = 300, E m = 0, and I m = 0.…”