Before revision:  Third, by introducing randomness to some of the system coefficients, we study the model's sensitivity to these parameters and provide bounds and intervals of confidence for the model simulations .
After revision:  Existence, uniqueness and positivity results are proved using a nonstandard convex analysis-based approach .