The equations for the Vacamole model in the vignette are as follows: $$dH = \eta I - \gamma H - \eta_2 H - \omega H$$ However, the equations in the header are: ```cpp dxdt.col(7) = iToH - hToR - hToD; // ηI - γH - ωH ```