  In this paper, we propose and analyze a new fractional order malaria model
with temporary immunity. We extend the existing malaria model to include
fractional derivatives. Basic reproduction number is computed using the next
generation matrix method. Model equilibria are determined and their local
asymptotic stability analysis are considered. An efficient Adams-type
predictor-corrector method is used to solve the initial value problem. Finally,
phase plane portraits are generated to describe the dynamics between
susceptible fractions of human hosts and infected fractions of mosquito
vectors.
