There are many joint density functions of N randoms for which each of the randoms is uniformly distributed (on, for concreteness, the interval $[0,1)$), yet they variables are not independent. The formula will depend on whi8ch jdf you choose.
Here is an extreme example: The $N$ variables are labeled $X_k | k = 0 \ldots N-1$ and variable $X_0$ is chosen as a uniform variate on $(0,1)$. Then for all $k > 0$,
$$
X_k = \left( X_0 + \frac{k}{N} \right) \pmod 1
$$
(where $a \pmod 1$ means the fractional part of $a$; $1.6 \pmod 1 = 0.6$).
Here the jdf is
$$
\text{jdf }(X_0 \ldots X_{N-1}) = \prod_{k=1}^{N-1} \left[ \delta\left(X_k - \frac{k}{N} -X_0 \right) + \delta\left(X_k - \frac{k}{N} - X_0 -1 \right)\right]
$$