I am working on a more detailed version of a proof that has appeared on Math Stack before. Requests for a proof of the following theorem have been posted several times, but I am asking about the specific strategy suggested in (1) below.
Theorem.
Let $f$ be any real-valued, continuous function with domain $\mathbb{R}$. Then $\Delta(f) := \{x\in \mathbb{R}: f'(x) \in \mathbb{R}\}$ is Borel measurable.
In other words, the set of points of finite differentiability belongs to the sigma-algebra generated by the closed sets.
I am trying to flesh out the proof given in (1) Show that $\Delta(f)$ is a $F_{\sigma\delta}$-set, where it is suggested that $\Delta(f)$ is exactly $$S\,=\,\bigcap_{k=1}^{\infty} \;\bigcup_{n=1}^{\infty} \;\; \bigcap_{0 < |\eta| < \frac{1}{n}} \;\; \bigcap_{0 < |\delta| < \frac{1}{n}} \; \left\{x:\;\; \left| \frac{f(x + \eta) - f(x)}{\eta} \; - \; \frac{f(x + \delta) - f(x)}{\delta} \right| \; \leq \frac{1}{k}\right\}.$$
I have showed that $\Delta(f)$ is a subset of the above set; I outline the reasoning below. For any fixed and non-zero $\eta$ and $\delta$, the function $$x \mapsto \left| \frac{f(x + \eta) - f(x)}{\eta} \; - \; \frac{f(x + \delta) - f(x)}{\delta} \right|$$ is continuous, so that the set $\{x: \cdots\}$ appearing above is closed, being a continuous pre-image of a closed set. So if $S=\Delta(f)$ the theorem is proved.
What I am stuck on is showing $S\subset \Delta(f)$. Set $S$ doesn't "know" about the value of the derivative, which is why I think it is harder. It seems (from looking at other approaches posted) that we need to pass to the rationals (maybe something like this: https://math.stackexchange.com/a/1395360), or build a Cauchy sequence, or a sequence of nested compacts. It makes me wonder, more generally, about sufficient conditions for differentiability at $x_0$ that do not involve pre-knowledge of $f'(x_0)$. What is a good way to rigorously show the inclusion $S\subset \Delta(f)\,$?
Sketch of proof that $\Delta(f)\subset S$:
Let $x \in \Delta(f)$, and let $k\geqslant 1$ be fixed. Let $\epsilon = 1/k>0.$ Let $t=f'(x)$.
$\exists n\geqslant 1,$ $$\left| \frac{f(x + \eta) - f(x)}{\eta} -t \right| \; \leqslant \; \frac{\epsilon}{2} \;\text{ whenever }\; 0 \lt |\eta|\lt \frac1n.$$ From here it is just a matter of applying the triangle inequality: $$\left| \frac{f(x + \eta) - f(x)}{\eta} - \frac{f(x + \delta) - f(x)}{\delta} \right| $$ $$\leqslant \left| \frac{f(x + \eta) - f(x)}{\eta} - t\right| + \left| \frac{f(x + \delta) - f(x)}{\delta} -t \right| \leqslant \epsilon=\frac{1}{k}.$$
To finish the proof of the theorem, we need only show $S\subset \Delta(f)$, where $$S\,=\,\bigcap_{k=1}^{\infty} \,\bigcup_{n=1}^{\infty} \; \bigcap_{0 < |\eta| < \frac{1}{n}} \; \bigcap_{0 < |\delta| < \frac{1}{n}} \left\{x:\, \left| \frac{f(x + \eta) - f(x)}{\eta} \; - \; \frac{f(x + \delta) - f(x)}{\delta} \right| \, \leq \frac{1}{k}\right\},$$ $$\Delta(f)\,=\,\bigcup_{t\in \mathbb{R}}\,\bigcap_{\epsilon>0}\,\bigcup_{\zeta>0} \; \bigcap_{0 < |h| < \zeta} \left\{x:\, \left| \frac{f(x + h) - f(x)}{h} \; - \; t\right| \, \leq \epsilon\right\}.$$
(I acknowledge the suggestion of HallaSurvivor, who indicated that the Cauchy sequence idea is one way to succeed here.)
Let $x\in S$ be fixed. Let $g_x(\theta):\mathbb{R}\setminus\{0\}\to\mathbb{R}$ be $$g_x(\theta)=\frac{f(x+\theta)-f(x)}{\theta}.$$ The sequence $g_x(1),g_x(\frac12),g_x(\frac13),\dots$ is Cauchy, because $x\in S$. Let $t$ be the limit of this sequence.
Let $\epsilon$ be given. Let $k>2/\epsilon,$ and let $n\geqslant1$ correspond to $k$ in the definition of $S$. Let $m>n$ be such that $$\left|g_x\left(\frac1m\right)-t\right|<\frac\epsilon2.$$ Let $\zeta$ be such that $0<\zeta<\frac1m,$ and let $0<|h|<\zeta.$ Now estimate
$$\left| \frac{f(x + h) - f(x)}{h} -t \right| \; = \; \left|g_x(h) -t \right|$$ $$ \leqslant \left|g_x(h)- g_x\big(\frac1m\big)\right| + \left|g_x\big(\frac1m\big)-t\right| < \frac{1}{k}+\frac\epsilon2<\epsilon,$$ because $|\eta|=|h|<\zeta<\frac1n$ and $|\delta|=\frac1m<\frac1n.$