Ostrogradsky-método de Hermite. Como he escrito en más detalle aquí, podemos utilizar el Ostrogradski-Hermite método para integrar una fracción racional $P(x)/Q(x)$ sin descomposición en fracciones parciales y sin encontrar múltiples raíces del denominador. Usted también puede encontrar una descripción simple de este método aquí, inciso 4.5.2, así como algunos ejercicios (1301-1304).
Suponga que $\deg P<\deg Q$. Entonces, existen polinomios $$P_{1},\quad Q_{1}=\gcd \left\{ Q,Q^{\prime }\right\} ,\quad P_{2}\quad \text{and}\quad Q_{2}=Q/Q_{1},$$ with $\deg P_{1}<\gr Q_{1}$, $\gr P_{2}<\gr Q_{2}$, que
\begin{equation*}
\int \frac{P}{Q}\, dx=\frac{P_{1}}{Q_{1}}+\int \frac{P_{2}}{
Q_{2}}\, dx.\tag{1}
\end{ecuación*}
Puesto que el integrando es
\begin{equation*}
\frac{P}{Q}=\frac{\left( x^{2}-3x+1/3\right) ^{2}}{\left(
x^{3}-x+1\right) ^{2}},\tag{2}
\end{ecuación*}
y
\begin{eqnarray*}
Q &=&\left( x^{3}-x+1\right) ^{2} \\
Q^{\prime } &=&2\left( x^{3}-x+1\right) \left( 3x^{2}-1\right) \\
Q_{1} &=&\gcd \left\{ Q,Q^{\prime }\right\} =x^{3}-x+1,\\
Q_{2}&=&\frac{Q}{Q_{1}}=x^{3}-x+1,
\end{eqnarray*}
podemos ampliar la integral dada como
\begin{equation*}
\int \frac{\left( x^{2}-3x+1/3\right) ^{2}}{\left( x^{3}-x+1\right) ^{2}}\, dx=
\frac{Ax^{2}+Bx+C}{x^{3}-x+1}+\int \frac{Dx^{2}+Ex+F}{x^{3}-x+1}\, dx,\tag{3}
\end{ecuación*}
donde $A$, $B$, $C$, $D$, $E$, $F$ son constantes. Se pueden encontrar por
la diferenciación $(3)$, la reducción de ambos lados a un denominador común e igualando los coeficientes de los poderes de $x$ de los numeradores. Por lo tanto, tienen
\begin{eqnarray*}
\frac{\left( x^{2}-3x+1/3\right) ^{2}}{\left( x^{3}-x+1\right) ^{2}} &=&
\frac{\left( 2Ax+B\right) \left( x^{3}-x+1\right) -\left( 3x^{2}-1\right)
\left( Ax^{2}+Bx+C\right) }{\left( x^{3}-x+1\right) ^{2}}\\ &&+\frac{\left(Dx^{2}+Ex+F\right)\left(x^{3}-x+1\right)}{\left( x^{3}-x+1\right) ^{2}}.
\end{eqnarray*}
Desde $\left( x^{2}-3x+1/3\right) ^{2}=x^{4}-6x^{3}+(29/3)x^{2}-2x+1/9$, después de algunos álgebra obtenemos
\begin{align}
x^{4}-6x^{3}+(29/3)x^{2}-2x+1/9&= Dx^{5}+\left( -A+E\right) x^{4}+\left( -2B-D+F\right) x^{3}\\
&\qquad+\left(-A-3C+D-E\right) x^{2}+\left( 2A+E-F\right) x\\
&\qquad+\left( B+C+F\right),
\end{align}
lo que significa que las constantes satisfacen el siguiente sistema de ecuaciones
\begin{equation*}
\left\{
\begin{array}{l}
D =0,\\ -A+E=1,\\ -2B-D+F=-6, \\
-A-3C+D-E =\frac{29}{3},\\ 2A+E-F=-2, \\ B+C+F=\frac{1}{9},
\end{array}
\right.
\end{ecuación*}
cuya solución es
\begin{equation*}
A=-1,\, B=3,\, C=-\frac{26}{9},\, D=0,\, E=0,\, F=0.
\end{ecuación*}
Por lo tanto
\begin{equation*}
\int \frac{\left( x^{2}-3x+1/3\right) ^{2}}{\left( x^{3}-x+1\right) ^{2}}\, dx=
\frac{-x^{2}+3x- 26/9 }{x^{3}-x+1}+C.\tag{4}
\end{ecuación*}