I am trying to solve systems of nonlinear equations with fsolve.
Since the variables are the concentrations with physical meaning, all of them should be always positive or zero.
However some values during the iteration and final results are negative. Is there any trick to keep all the values positive or zero during iteration and final results?
Your comment would be appreciated.