$$ f(z) = \frac{z}{(\sin z)^2} $$ at $z_0 = 0$ (for the first four terms).
So I thought I knew what to do, but I don't. Since it appears to be an indeterminte form, could I by L'hopital turn it into:
$$\frac{1}{2(\sin z)(\cos z)} $$ but even here I am jammed
Why just the first four terms? Let's compute all of them. For first, since $z=0$ is a removable singularity for $\frac{\sin z}{z}$ we have that $z=0$ is a simple pole for $f(z)=\frac{z}{\sin^2 z}$, and the problem boils down to finding the Taylor series, in a neighbourhood of the origin, of: $$ g(z) = \left(\frac{z}{\sin z}\right)^2. \tag{1}$$ From the Weierstrass product: $$\frac{\sin z}{z}=\prod_{n\geq 1}\left(1-\frac{z^2}{n^2\pi^2}\right)\tag{2} $$ we have that any $z=m\pi$, for $m\in\mathbb{Z}\setminus\{0\}$, is a simple pole for $\frac{z}{\sin z}$, and by considering the logarithmic derivative of $(2)$ we get: $$ \cot z-\frac{1}{z} = \sum_{n\geq 1}\frac{2z}{z^2-n^2\pi^2},\tag{3} $$ a really useful identity, that leads to: $$ \frac{1-z\cot z}{2}=\sum_{n\geq 1}\frac{z^2}{n^2\pi^2-z^2}=\sum_{n\geq 1}\frac{\frac{z^2}{n^2\pi^2}}{1-\frac{z^2}{n^2\pi^2}}=\frac{\zeta(2)}{\pi^2}z^2+\frac{\zeta(4)}{\pi^4}z^4+\ldots\tag{4} $$ Rearranging we have: $$\cot z = \frac{1}{z}-2\sum_{n\geq 1}\frac{\zeta(2n)}{\pi^{2n}}z^{2n-1}\tag{5}$$ and now we just need to differentiate both terms with respect to $z$ and multiply them by $-z^2$ to get:
from which: $$ g(z)=\left(\frac{z}{\sin z}\right)^2 = 1+\frac{z^2}{3}+\frac{z^4}{15}+o(z^5) $$ readily follows.