Context
Consider a simplified version of the ODE model presented in Marini et al. 2017 (see section Model calibration of Additional file 1 in that paper).
In this model there are two populations of mosquitoes, that go through their different life stages: egg -> larva -> pupae -> female adult. For simplicity the pupae stage is omitted in this example without loss of information (see attached file 1). When evaluating this system in Wolfram Mathematica and obtaining the Jacobian we arrive at the steady states (see attached file 2).
Problem
If we focus on the scenario where all populations are present (i.e. Ec > 0, Ea > 0, Lc > 0, La > 0, Ac > 0, Aa > 0) and substitute the parameters with any values at this steady state, the populations of larvae, Lc and La, exceed their carrying capacities, Kc and Ka respectively.
Question
Why is this happening?
Can populations exceed the carrying capacity if they are only limited by it on a single life stage?
Or is there a problem in the implementation in Wolfram Mathematica?