The version of rigorous infinitesimals I've seen didn't let you have these identities, because e^2 != 0 (because it still maintains the field axiom xy = 0 => x=0 V y=0); but the argument still goes through by carrying the whole power series in place of these simplifications – you can write the higher order terms as e*f(x, e) and drop them when you evaluate st(...).