Generic stable algorithm implemented for polynomials and tested up to cubics
Fixes for the new substitution engine, preparing for quartic tests