I seek to prove that $$\sum_{n=1}^\infty \frac{1}{\lambda_{n}^2} = \frac{1} {10},$$ by applying Cauchy Theorem to $$ f(z) = \left(\frac{z\tan(z)}{z-\tan(z)}+\frac{3}{z}\right) \frac{1}{z^2},$$ where $\lambda_n$ are positive solutions to $\tan(x) = x$.
How can I solve this?
Let's see:
$f(z)$ has three types of poles $z_0=0$ , $z_{1,n}=n \pi/2$ and $z_{2,n}=\lambda_n$
It turns out that the poles at $z_{1,n}$ are removeable so the corresponding residue will yield zero (why is that? because $f(z)$ is exactly designed to do so). Now let us integrate $f(z)$ around a contour which consists of two lines just above and below the real axis (connected at infinity), with a sense of orientation which is clockwise. Due to the residue theorem we have
$$ \oint f(z)dz = -2\pi i [\text{res}(z_0)+\sum_{n=-\infty}^{\infty}\text{res}(z_{2,n})] $$
But on the other hand, because $tan(z)$ is a bounded function away from the real axis we can also close our contour of integration by adding big semicircles to the now disconnected lines. Because there will no now poles enclosed during this procedure we see that $\oint f(z)dz=0$ and therefore
$$ -\text{res}(z_0)=\sum_{n=-\infty}^{\infty}'\text{res}(z_{2,n}) \quad (1) $$
Here $\sum'$ means that we omit $0$ under the sum sign. The residue at $0$ is given by $\text{res}(z_0)=\frac{1}{5}$, which follows from a straightforward Laurent expansion.
The tricky part know is to calculate the residues at $\lambda_n$. First of all because the poles are simple we can use the
$$ \text{res}(\lambda_n)=\frac{g(\lambda_n)}{f'(\lambda_n)} \quad (2) $$
with $g(z)=\frac{z^2\tan(z)+3(z- \tan(z))}{z^3}$ and $f'(z)=1-\sec^2(z)=-\tan^2(z)$
using $\tan(\lambda_n)=\lambda_n$ we might rewrite (2) as $$ \text{res}(\lambda_n)=-\frac{\lambda_n^2 \lambda_n+ 3(0)}{\lambda_n^3\lambda_n^2}=-\frac{1}{\lambda_n^2} $$ Putting this into (1) we might observe that
$$ -\frac{1}{5}=-\sum_{n=-\infty}^{\infty}'\frac{1}{\lambda_n^2} $$
now we observe that the equation $\tan(x)=x$ is anti-symmetric under $x\rightarrow -x$ and all zeros come as pairs $\pm\lambda_n$. This allows us to rewrite
$$ \frac{1}{10}=\sum_{n=1}^{\infty}\frac{1}{\lambda_n^2} $$
Q.E.D