The equations are highly non-linear, and there could be multiple solutions. Is there a way to impose specific constraints on the solution?
Also, in many cases, it is very hard to converge to exact solution, so if MATLAB is used, how to set tolerances for solutions?