$\newcommand{\eps}{\varepsilon}$Sugerencia: Si $(x, y)$ $(x_{0}, y_{0})$ son puntos arbitrarios, entonces
$$
f(x, y) - f(x_{0}, y_{0})
= \bigl(f(x, y) - f(x, y_{0})\bigr) + \bigl(f(x, y_{0}) - f(x_{0}, y_{0})\bigr).
$$
En el primer término de la derecha, sólo $y$ varía; en el segundo, sólo $x$ varía.
Edit: Más de sugerencia, ya que el epsilontics son un poco delicados.
Suponer sin pérdida de generalidad que $f$ es no decreciente en $y$.
Fix$(x_{0}, y_{0})$$\eps > 0$. Uso la continuidad de $f$ $y$ a recoger $\delta_{2} > 0$ tal que
$$
|y - y_{0}| < \delta_{2} \text{ implica }
|f(x_{0}, y) - f(x_{0}, y_{0})| < \tfrac{1}{3}\eps.
$$
Ahora uso la continuidad en $x$ (y el "estándar mínimo de dos condiciones" modismo) para recoger $\delta_{1} > 0$ tal que
$$
|x - x_{0}| < \delta_{1} \text{ implica }
|f(x, y_{0} \pm \delta_{2}) - f(x_{0}, y_{0} \pm \delta_{2})| < \tfrac{1}{3}\eps.
$$
Si $|x - x_{0}| < \delta_{1}$$0 \leq y - y_{0} < \delta_{2}$, luego
$$
f(x, y) - f(x_{0}, y_{0})
= \bigl(f(x, y) - f(x_{0}, y)\bigr) + \bigl(f(x_{0}, y) - f(x_{0}, y_{0})\bigr).
\etiqueta{1}
$$
Ahora para la parte delicada: Desde $f$ es monótona en $y$,
\begin{align*}
\bigl(f(x, y) - f(x_{0}, y)\bigr)
&\leq \bigl(f(x, y_{0} + \delta_{2}) - f(x_{0}, y_{0})\bigr) \\
&= \bigl(f(x, y_{0} + \delta_{2}) - f(x_{0}, y_{0} + \delta_{2})\bigr)
+ \bigl(f(x_{0}, y_{0} + \delta_{2}) - f(x_{0}, y_{0})\bigr)
\tag{2a}
\end{align*}
(incremento $y$ en el primer término y la disminución en la segunda, a continuación, utilizar mi primera sugerencia), y
$$
\bigl(f(x_{0}, y) - f(x_{0}, y_{0})\bigr)
\leq \bigl(f(x_{0}, y+ \delta_{2}) - f(x_{0}, y_{0})\bigr).
\etiqueta{2b}
$$
Totalmente similar estimaciones de manejar el caso $-\delta < y - y_{0} < 0$.
Poner todo junto, si $|x - x_{0}| < \delta_{1}$$|y - y_{0}| < \delta_{2}$, luego
\begin{align*}
|f(x, y) - f(x_{0}, y_{0})|
&\leq \big|f(x, y) - f(x_{0}, y)\big| + \big|f(x_{0}, y) - f(x_{0}, y_{0})\big| \\
&\leq \big|f(x, y_{0} + \delta_{2}) - f(x_{0}, y_{0} + \delta_{2})\big|
+ \big|f(x_{0}, y_{0} + \delta_{2}) - f(x_{0}, y_{0})\big| \\
&\qquad + \big|f(x_{0}, y+ \delta_{2}) - f(x_{0}, y_{0})\big| \\
&< \eps.
\end{align*}