Hi, I have a nonlinear set of equations that I have to solve implicitly, it's a flow network where I can't escape the square law relationship between pressure and flow, and as a consequence an explicit solution is looking either very nasty indeed or infeasible. My simulation tools can handle a single tear variable (using Newton Raphson), but the situation has just got more complicated, and I now have to converge two separate tear variables at once to solve the equation set.

Can anyone point me to the name of a technique that can be used for this? I'd really prefer a codeable approach rather than a packaged tool, because the solution has to be embedded into existing C++ code.

TIA