Hello everyone, I am trying to develop a biological propagation model in Matlab.
The model is composed by 8 coupled ordinary differential equations, which describe how certain species evolve through the time. I have been using the ode45 solver. You can find the implementation at this link However, I am facing some difficulties about adding constraints to these equations. I would like to add some upper/lower bounds to the values, but apparently the solver is ignoring these constraints.
Anybody can help me?