Question: Non-negativity constraints in solve or fsolve command in Maple

Hi Mapleprimes,

I have tried to search solutions online without luck. My problem is this. Suppose I have three equations:

f(x,y,z)

g(x,y,z)

h(x,y,z)

I would solve using solve({f,g,h},{x,y,z}) which would give me solutions x*, y* and z*. I need to assign the solutions to use them in subsequent computations. I would like to impose non-negativity constraints, however. Hence, the solution of x is max(x*,0). Obviously if a non-negativity constraint binds, this affects the solutions of y* and z*. I would like to assign the solutions taking this into account. How would you propose I do this? Bear in mind that the solution of x* might be a-b where the relative magnitudes of paramenters a and b are not given. Further down in my code, I would like to assign values to parameters and then let Maple give me the solutions with the non-negativity constraints. I hope I am making sense.

Thanks very much in advance.

Best,

Christian