Construct the infinitesimal generator matrix for (individual) adult female infection dynamics. Adult females follow SEI (Susceptible-Exposed-Infectious) style dynamics with a Gamma distributed EIP, with a mean duration 1/q and variance 1/nq^2 (following shape-scale parameterization, EIP ~ Gamma(n,1/nq)). This function only constructs the rate matrix for either a single mosquito or cohort that all emerged at the same time (the rate matrix for a population with emergence is infinite in dimension).
make_Q_SEI(q, n, mu, c, a, x)
rate matrix for a single (emergence) cohort of SEI mosquito
related to scale parameter of Gamma distributed EIP (1/q is mean length of EIP)
shape parameter of Gamma distributed EIP
mosquito mortality rate
human to mosquito transmission efficiency
human biting rate
prevalence of disease in humans