Let us stick to $\mathbb{R}$ and $\mathcal{B}$, the Borel Sigma field. Given that $\mathbb{F}$ is a sigma field of $\Omega$, consider a random variable $X$
$X:\begin{bmatrix} \Omega \\ \mathbb{F} \\ P \end{bmatrix} \to \begin{bmatrix} \mathbb{R}\\\mathcal{B}\\ P_X \end{bmatrix}$
Definition: $S_X:=$ smallest closed subset of $\mathbb{R}$ such that $P_X(S_X)=1$. Call this support.
Definition: $U_X:=\{x\in\mathbb{R} \ | \ \forall \ \varepsilon>0, \ P_X(V_{\varepsilon}(x))>0 \}$ (where $V_\varepsilon(x)$ is the $\varepsilon$ neighbourhood of $x$)
I was able to show that the definition $S_X$ is equivalent to the definition $U_X$.
Consider now, that $X$ is also a continuous random variable, which means: $$P_X(a,b)=\int_a^b p_d(t)dt$$ where $p_d(t)$ is the density function.
this answer, at the very end, defines the support $G_X$, for this case, as follows: $$G_X=\text{closure}\{x\in\mathbb{R}\ | \ p_d(x)>0 \}$$
Now I want to show that, for the case of continuous random variable, $G_X \equiv S_X\equiv U_X$.
It was rather easy to show that $ S_X \subseteq G_X$, but the reverse assertion is proving to be trouble some. I suspect some application of Lagrange Mean value theorem or something of that sort. Any ideas?