This is inteded to be a simultaneous non linear solver with units. The solver would be written in C++ and the front end in Python.