Esta es una forma fácil de hacerlo: $$\frac{1}{x^3(x-1)}=\frac{Ax^2+Bx+C}{x^3}+\frac{D}{x-1}$$
Elegí empezar en $x^2$ para el primero ya que el grado del denominador es $3$ y por el otro $x^0$ porque el denominador tiene grado $1$ (por lo que hay que ir hasta abajo empezando por $x^{\text{denominator degree}-1}$ ).
Combina las dos fracciones en una fracción con el denominador original: $$\frac{(Ax^2+Bx+C)(x-1)+D(x^3)}{x^3(x-1)}$$
Y recuerda que el numerador era igual a $1$ : $$(Ax^2+Bx+C)(x-1)+D(x^3)=1$$ Distribuir: $$Ax^3+Bx^2+Cx-Ax^2-Bx-C+Dx^3=1$$ Combina términos similares: $$(A+D)x^3+(B-A)x^2+(C-B)x-C=1$$
Y adivina qué $$(A+D)x^3+(B-A)x^2+(C-B)x-C=0x^3+0x^3+0x+1$$
Por lo tanto, tenemos: $$\begin{cases}A+D=0\\B-A=0\\C-B=0\\-C=1\end{cases}$$
Lo que te da: $$\begin{cases}A=-1\\B=-1\\C=-1\\D=1\end{cases}$$
Y así es: $$\frac{-x^2-x-1}{x^3}+\frac1{x-1}=-\frac1{x}-\frac1{x^2}-\frac1{x^3}+\frac1{x-1}$$
Que es fácil de integrar: $$\begin{align}&\int\left(-\frac1{x}-\frac1{x^2}-\frac1{x^3}+\frac1{x-1}\right)dx\\=&-\ln x+\frac1{x}+\frac1{2x^2}+\ln(x-1)\end{align}$$