system: only number equations when need to solve systems

-> allow to define options after solver creations
-> constraints must be alive until system has been solved !
-> mostly necessary for python binding, but also usefull in c++