Get rid of xreplace, simplify, etc

Unnecessary complexity: Replacements should be done
directly in the objects producing constraint expressions instad
Add pyomo support

Minimal working example in test.py. Needs cbc on path to work.