Let $f(z)=\frac{1}{1-\cos z}$. I want to calculate the Laurent series of $f$ around $z=0$, up to the coefficient of $z^6$.
It's easy to see that $g(z)=z^2f(z)$ (Defined as its limit at $0$ is holomorphic since $z=0$ is a pole of order $2$. Furthermore, it's easy to see that $g(0)=2$. I can calculate the Maclaurin series of $g(z)$ until $z^8$ and divide by $z^2$ to get the result. Also, since $g(z)$ is an even function, the odd coefficients are all $0$. And since $g(z)$ is of class $C^\infty$ on a small enough disk, I think I can just (In theory) calculate $g^{(n)}(0)=\lim_{z\to 0} g^{(n)}(z)$ by L'Hospital rule.
The problem is that the expression $g^{(n)}(z)$ becomes more and more complicated when $n$ grows, and I must calculate all its even derivatives at $0$.
So I wanted to know if there is a more efficient way than calculating all the even derivatives $g^{(n)}(z)$ until $n=8$, and then the derivatives of the numerator and denominator of each expression until one of them isn't $0$ at $0$. For example, in the denominator there will be a $(1-\cos z)^{2^n}$, which has a zero of order $2^{n+1}$ which must cancel with some zero in the numerator. It's too big a calculation just to use L'Hospital until $n=8$.
There are some related questions about computing a general term of the series; I'm just thinking about a moderately fast way to compute some of the coefficients.