Convergence of double series over $\mathbb{Z}$ I am attempting to show that the following double series converges
\begin{equation}
\sum\limits_{a,b \in \mathbb{Z}} \frac{1}{\left(\sqrt{a^2 +b^2} - M\right)^3}
\end{equation}
where $M \in \mathbb{R}^{+}$ is a constant. I've never dealt with double series before, but this is what I have attempted so far:
First let's define $$u_{ab} := \frac{1}{\left(\sqrt{a^2 +b^2} - M\right)^3}  \hspace{10mm}\left(\text{$\textbf{Notice:}$ $u_{a,b} = u_{-a,b}=u_{a,-b}=u_{-a,-b}$}\right)$$
Then our double series becomes
\begin{align*}
\sum\limits_{a=-\infty}^{\infty}\left(\sum\limits_{b=-\infty}^{\infty} u_{ab}\right) &= \sum\limits_{a=-\infty}^{\infty}\left(u_{a0} + 2\sum\limits_{b=1}^{\infty} u_{ab}\right) &&\text{(by $\textbf{Notice}$)}\\[10pt] 
&= \sum\limits_{a=-\infty}^{\infty}u_{a0} + 2\sum\limits_{a=-\infty}^{\infty}\left(\sum\limits_{b=1}^{\infty} u_{ab}\right)\\[10pt]
&=u_{00} +2 \left(\sum\limits_{a=1}^{\infty}u_{a0}\right) + 2\left[\sum\limits_{b=1}^{\infty}u_{0b} + 2\sum\limits_{a=1}^{\infty}\sum\limits_{b=1}^{\infty}u_{ab}\right]\\[10pt]
&=u_{00} + 2\left(\sum\limits_{a=1}^{\infty} u_{a0} + \sum\limits_{b=1}^{\infty}u_{0b}\right) + 2\sum\limits_{a,b=1}^{\infty}u_{ab}\\[10pt]
&=u_{00} + 4\sum\limits_{a=1}^{\infty}u_{a0} + 2\sum\limits_{a,b=1}^{\infty}u_{ab} &&\left(\text{since $u_{ab}=u_{ba}$}\right)
\end{align*}
Firstly: is my analysis above correct?

Secondly: I know that $\sum\limits_{a=1}^{\infty}u_{a0}$ converges, however the latter series is a double series, I have no experience/knowledge on how to show convergence for this. Any help? Intuition points towards some sort of $2d$ integral test? But I'm very unsure about that.

Your analysis so far looks good to me. I think it is worth stating explicitly in your question that $M$ is not equal to $\sqrt{a^2+b^2}$ for any integers $a,b$.
Now, with that out of the way: we have $a^2+b^2\ge 2ab$ for all $a,b$. So for $a\ge \frac{M^2}{(\sqrt 2-1)^2}$ or $b\ge \frac{M^2}{(\sqrt 2-1)^2}$, $|\sqrt{a^2+b^2}-M|\ge \sqrt{ab}$. So the sum over this region is bounded by
$$\sum_{a,b\ge 1} \left(\frac{1}{ab}\right)^\frac32=\sum_{a\ge 1} \left(\frac{1}{a}\right)^\frac32 \sum_{b\ge 1} \left(\frac{1}{b}\right)^\frac32$$
which converges.
