$X_1, X_2,...,X_n$ is a random sample from $U(0,\theta)$. Find $E(X_{(n)}-X_{(1)})$.
I attempted this question by first finding the CDF of $X_{(n)}-X_{(1)}$ using the formula: $$F_{U}(u)= n\int_0^\theta f(x)[F(u+x)-F(x)]^{n-1}dx$$ Where $U=X_{(n)}-X_{(1)}$, $f(x)$ is the PDF of $U(0,\theta)$ and $F(x)$ is the CDF of $U(0,\theta)$.
Using this formula I obtained: $$F_{U}(u)=n\frac{u^{n-1}}{\theta^{n-1}}$$
Now, using the following formula for non-negative continuous random variables: $$E(U) = \int_0^\theta (1-F_{U}(u))du$$
I obtained: $$E(U) = \int_0^\theta \Big(1-n \Big(\frac{u}{\theta}\Big)^{n-1}\Big)du=\theta - \frac{n}{\theta^{n-1}}\Big[\frac{u^n}{n}\Big]_{0}^{\theta}=\theta-\theta=0$$
However, if I break the expectation and calculate individually, $$E(U) = E(X_{(n)})-E(X_{(1)})$$ I get the answer as: $$E(U) = \frac{n-1}{n+1}\theta$$ which I believe is the right answer. Can someone please explain why is the former method giving an incorrect answer?
EDIT: Proof of CDF of $X_{(n)}-X_{(1)}$. (General case when $X_{i}$'s are defined over the range $(-\infty,\infty))$
$X_{(n)}-X_{(1)}$. (General case when $X_{i}$'s are defined over the range $(-\infty,\infty))$" />