I have just looked at about 10 different Laurent series in the past days, but I cannot solve $f(z)=\frac{1}{z^2}$.
I attempt, $|z|>0$, with u=$z^2+1$:
\begin{equation} \frac{1}{z^2}=\frac{1}{z^2-1+1}=\frac{1}{u-1}=-\frac{1}{1-u} \end{equation}
From here, assuming there exists substitution method at all for Laurent series, I get:
\begin{equation} -\frac{1}{1-u}=\sum_{n=-\infty}^\infty(-1)^nu^n=\sum_{n=-\infty}^\infty(-1)^n(z^2+1)^n \end{equation}
but this is not correct. The correct result is :
\begin{equation} -\frac{1}{1-u}=\sum_{n=0}^\infty(-1)^n(1+n)(-1+z)^n \end{equation}
It seems there is a rule of differentiation here, so I try this:
Put $g(z)=\frac{1}{z}$, and $f(z)=\frac{1}{z^2}$, we see that $g'(z)=-f(z)$
So
\begin{equation} g(z)=\frac{1}{z}=\frac{1}{1-z-1}=-\frac{1}{z+1}\frac{1}{-\frac{1}{z+1}+1}=\frac{1}{z+1}\frac{1}{1-(-\frac{1}{z+1})} \end{equation}
This gives
\begin{equation} \frac{1}{z+1}\sum_{n=0}^{\infty}(-1)^n\bigg(\frac{1}{z+1}\bigg)^n = \sum_{n=0}^{\infty}(-1)^{n-1}\frac{1}{(z+1)^n} \end{equation}
Since $g'(z)=-f(z)$, then $-g'(z)=f(z)$, so I differentiate the series:
\begin{equation} f(z)=\frac{1}{z^2}=\sum_{n=0}^{\infty}\bigg((-1)^{n-1}\frac{1}{(z+1)^n}\bigg)' =\sum_{n=0}^{\infty}(-1)^{n-1}(-n)(z+1)^{-n-1} \end{equation}
But this is not entirely correct, though we are getting closer. What is the wrong in this last procedure here?