I recently read over Nathaniel Martin's thesis, which has the following theorem on Page 19
Theorem 3.1. Let $a$ be any point in $\mathbf{R}$ and let $\delta$ be any real number with $0 \leq \delta \leq 1$. Then there exists a set $G$ whose density exists at $a$ and has the value $\delta$.
Here density defined to be the
$$ D_E(x) = \lim_{\epsilon\to 0} \frac{m(E \cap (x-\epsilon, x+\epsilon))}{2\epsilon}. $$
I can mostly follow his proof, but I was wondering if the problem is meaningfully simplified if instead of letting $a$ be any point in $\mathbf{R}$, it is instead set to be $0$? I am looking for an application of Vitali's covering lemma in the proof, in a similar way to the proof for the Lebesgue density theorem in the line, but am having a hard time progressing.
I saw elsewhere a suggestion for a constructive approach, using $$E_n = \pm (1/(n+1), (n+\alpha)/(n(n+1))]$$
I am struggling to find an application for these pairs of intervals, though.