Estoy leyendo este documento.
En este artículo después de la definición de fuerte derivado de Knuth va a calcular la derivada de $x^n$. Allí se utiliza la definición de los fuertes de derivados para expandir $(x+\epsilon)^{n+1}$ el siguiente,
$(x+\epsilon)^{n+1} = (x+\epsilon)[x^n + d_n(x)\epsilon + \mathcal{O}(\epsilon^2)]$
cuando me ampliar el lado derecho puedo conseguir
$x^{n+1} + (x d_n(x) + \ x^n)\epsilon + \color{red}{x \ \mathcal{O}(\epsilon^2) + \epsilon^2 \ d_n(x) + \epsilon \ \mathcal{O}(\epsilon^2)}$
Pero en Knuths los cálculos de la parte roja es sólo $\mathcal{O}(\epsilon^2)$.
La pregunta es ¿cómo? No sé cómo resolver esto en detalle.