Fixes double-counting bug in constant_renumber

* dsolve((x+1)*f(x).diff(x) - f(x)- 1,f(x)) now returns correctly
* tests in test_constantsimp and test_ode updated
* doctest corrected
* minor readability editing