A slightly different approach which can be very useful e.g. when
applying the Egorychev method is to use the fact that for a rational
function all residues including the residue at infinity sum to
zero. (I don't have a reference ready at this time.) The residue at
infinity is given by
$$\mathrm{Res}_{z=\infty} f(z)
= - \mathrm{Res}_{z=0} \frac{1}{z^2} f(1/z).$$
This means the sum of the residues at the finite poles is given by
$$\mathrm{Res}_{z=0} \frac{1}{z^2} f(1/z)$$
which in the present case yields
$$\mathrm{Res}_{z=0} \frac{1}{z^2}
\frac{a_0+a_1/z+\cdots+a_{n-1}/z^{n-1}}{b_0+b_1/z+\cdots+b_n/z^n}
\\ = \mathrm{Res}_{z=0} \frac{1}{z^2}
\frac{a_0z^n+a_1z^{n-1}+\cdots+a_{n-1}z}{b_0z^n+b_1z^{n-1}+\cdots+b_n}
\\ = \mathrm{Res}_{z=0} \frac{1}{z}
\frac{a_0z^{n-1}+a_1z^{n-2}+\cdots+a_{n-1}}
{b_0z^n+b_1z^{n-1}+\cdots+b_n}.$$
This is
$$[z^0] \frac{a_0z^{n-1}+a_1z^{n-2}+\cdots+a_{n-1}}
{b_0z^n+b_1z^{n-1}+\cdots+b_n}$$
because the denominator is not zero at $z=0$ as per the problem
definition. With no pole at zero the constant coefficient can be
obtained from the ordinary Taylor series and is equal to
$$\left.\frac{a_0z^{n-1}+a_1z^{n-2}+\cdots+a_{n-1}}
{b_0z^n+b_1z^{n-1}+\cdots+b_n}\right|_{z=0} = \frac{a_{n-1}}{b_n}.$$
The rule that residues sum to zero may then be used in a variety of
creative ways, for example at this MSE link
I where it is
applied twice, with the residue at infinity vanishing in the second
application as per the integral formula. The same technique is used at
this MSE link II.
The Egorychev method with a non-zero residue at infinity is applied at
this MSE limk III
(streamlined proof).