I wish to calculate the integral $$\int_{|z|\mathop=2}\sqrt{z^4-z}\,dz$$
$$z^4-z=z(z-1)(z-e^{2i\pi/3})(z-e^{-2i\pi/3})$$ We must make branch cuts which go through the branch points. The branch points of $\sqrt{z^4-z}$ are $z=0,1,e^{2i\pi/3},e^{-2i\pi/3}$. All of these lie inside this contour, so I first thought that I must work out the residue at each of these points, and then use the residue theorem. However there is no nice way that I could think of to calculate these residues.
Then it occurred to me that I am not even sure if they exist - these are not poles, they are branch points with branch cuts through them. I then saw that the question has a hint saying to use Laurent Series. Most of the Laurent series I have ever seen have been functions with (sort of) Taylor series, but with some change of variables, or multiplied by some factor of $1/z^n$. The point is - I don't know any (efficient) way of finding the Laurent Series for this function.
My questions: How can I
Define these branch cuts,
Find the different Laurent Series, and find the value of the coefficient of $1/z$, thus finding value of the integral?
Also:
I thought Laurent Series existed within an annulus of the point about which they exist - how is this possible for a function with branch cuts? Is there something I am confusing here? These two ideas I had seem to contradict.
If the Laurent Series do exist, then is it possible to calculate the residues directly? (Also assuming these exist)
We can cut the plane such that the branch cuts join pairs of the branch points. For example, the cuts are from $0$ to $1$ and from $e^{i2\pi/3}$ and $e^{-i2\pi/3}$.
Note that with this choice of branch cuts, $f(z)=\sqrt{z^4-z}$ is analytic in any annulus $1<|z|<R$. Therefore, we can deform the contour $|z|=2$ to $|z|=R\to \infty$.
The integral of interest becomes
$$\begin{align} \oint_{|z|=2}f(z)\,dz&=\oint_{|z|=R}f(z)\,dz\\\\ &=\oint_{|z|=R} z^2\sqrt{1-\frac{1}{z^3}}\,dz\\\\ &=\oint_{|z|=R} z^2\left(1-\frac{1}{2z^3}+O\left(\frac{1}{z^6}\right)\right)\,dz\\\\ &=\int_0^{2\pi}R^2e^{i2\phi}\left(1-\frac{1}{2R^3e^{i3\phi}}+O\left(\frac{1}{R^6}\right)\right)\,iRe^{i\phi}\,d\phi\\\\ &=-i\pi+O\left(\frac1{R^3}\right)\\\\ &\to -i\pi \,\,\text{as}\,\,R\to \infty \end{align}$$