calculate only updated force of infection at present, this is the only state-dependent per capita rate maybe more efficient than modifying & returning the whole matrix

update_foi(state, params, beta)

Arguments

state

named vector of states

params

named vector of parameters

beta

vector or matrix of transmission rates, where (length(beta) | ncol(beta)) == length(state)