mature_survivor_rate = 0.888 INIT tainted_mosquitos = 0 TRANSIT TIME = 12 INFLOW LIMIT = ° CAPACITY = ° prob_bit = 0.3 INIT vectors = 0 INIT uninfected_mosquitoes = 300 mosquitoes = init(uninfected_mosquitoes) + init(tainted_mosquitos) + init(vectors) prob_vector = vectors / mosquitoes INIT uninfected_humans = 300 INIT human_hosts = 1 recovery_rate = .3 recovered = human_hosts * recovery_rate birth_rate = 0.0003 INIT immune = 0 humans = human_hosts + immune + uninfected_humans births = birth_rate * humans flow_to_host = prob_bit * prob_vector * uninfected_humans human_death_rate = 0.0003 uninfected_deaths = human_death_rate * uninfected_humans immunity_rate = .01 flow_to_immune = human_hosts * immunity_rate malaria_induced_death_rate = .005 human_host_deaths = (human_death_rate + malaria_induced_death_rate) * human_hosts prob_bite_human = .3 prob_host = human_hosts / humans mosquito_birth_rate = .03 mosquito_births = mosquito_birth_rate * mosquitoes infected = prob_bite_human * prob_host * uninfected_mosquitoes mosquito_death_rate = 0.01 uninfected_mosquito_deaths = mosquito_death_rate * uninfected_mosquitoes immune_deaths = human_death_rate * immune mature_survivor_rate = 0.888 vector_deaths = mosquito_death_rate * vectors vectors(t) = vectors(t - dt) + (tainted_to_vector - vector_deaths) * dt uninfected_mosquitoes(t) = uninfected_mosquitoes(t - dt) + (mosquito_births - infected - uninfected_mosquito_deaths) * dt uninfected_humans(t) = uninfected_humans(t - dt) + (recovered + births - flow_to_host - uninfected_deaths) * dt human_hosts(t) = human_hosts(t - dt) + (flow_to_host - flow_to_immune - recovered - human_host_deaths) * dt immune(t) = immune(t - dt) + (flow_to_immune - immune_deaths) * dt tainted_mosquitos(t) = tainted_mosquitos(t - dt) + (infected - tainted_to_vector) * dt mosquitoes = init(uninfected_mosquitoes) + init(tainted_mosquitos) + init(vectors) prob_vector = vectors / mosquitoes recovered = human_hosts * recovery_rate humans = human_hosts + immune + uninfected_humans births = birth_rate * humans flow_to_host = prob_bit * prob_vector * uninfected_humans uninfected_deaths = human_death_rate * uninfected_humans flow_to_immune = human_hosts * immunity_rate human_host_deaths = (human_death_rate + malaria_induced_death_rate) * human_hosts prob_host = human_hosts / humans mosquito_births = mosquito_birth_rate * mosquitoes infected = prob_bite_human * prob_host * uninfected_mosquitoes uninfected_mosquito_deaths = mosquito_death_rate * uninfected_mosquitoes immune_deaths = human_death_rate * immune immune_deaths = human_death_rate * immune vector_deaths = mosquito_death_rate * vectors