The sequence of real numbers $\{a_n\}$, $n \in \mathbb{N}$ satisfies the following condition: $a_{n+1}=a_n(a_n+2)$ for any $n \in \mathbb{N}$. Find all possible values for $a_{2004}$.
1. Given the sequence $\{a_n\}$, we know that $a_{n+1} = a_n(a_n + 2)$ for any $n \in \mathbb{N}$.
2. To simplify this recurrence relation, we add 1 to both sides:
   \[
   a_{n+1} + 1 = a_n(a_n + 2) + 1
   \]
3. Expanding the right-hand side:
   \[
   a_{n+1} + 1 = a_n^2 + 2a_n + 1
   \]
4. Notice that the right-hand side can be written as a perfect square:
   \[
   a_{n+1} + 1 = (a_n + 1)^2
   \]
5. Let us define a new sequence $\{b_n\}$ such that $b_n = a_n + 1$. Then the recurrence relation becomes:
   \[
   b_{n+1} = (a_n + 1)^2 = b_n^2
   \]
6. This implies that $b_{n+1} = b_n^2$. Starting from some initial value $b_1 = a_1 + 1$, we can see that:
   \[
   b_2 = b_1^2, \quad b_3 = b_2^2 = (b_1^2)^2 = b_1^4, \quad \text{and so on}
   \]
7. Generally, we have:
   \[
   b_n = b_1^{2^{n-1}}
   \]
8. For $n = 2004$, we get:
   \[
   b_{2004} = b_1^{2^{2003}}
   \]
9. Since $b_1 = a_1 + 1$, we have:
   \[
   b_{2004} = (a_1 + 1)^{2^{2003}}
   \]
10. The value of $b_{2004}$ can be any non-negative real number because $b_1$ can be any non-negative real number. Therefore, $a_{2004} = b_{2004} - 1$ can be any real number in the interval $[-1, \infty)$.

The final answer is $\boxed{[-1, \infty)}$