Edit: Explicit calculations now included at the bottom. Thanks go to Dylan Moreland for suggesting I use Sage to compute Groebner bases.
Edit 2: I have now included a lemma with proof towards the bottom that extends this method so that it can show that an ideal is primary, even if it is not prime.
Using any techniques I'm familiar with, proving an ideal prime (radical) is generally painful. The method I use below is tailored to your specific example but can be adapted to work more generally. If we want to show that $I$ is prime (radical), the easiest way IMHO is to show that $K[x_0,\ldots,x_3]/I$ is an integral domain (has only nilpotent zero-divisors). The methods that follow work to prove an ideal prime, but not to prove an ideal radical.
Note: I abuse notation, using $x_0,x_1$, etc. for both elements of $K[x_0,\ldots,x_3]$ and their images under various homomorphisms. I do not believe this will result in any confusion, but if it does I apologize and will attempt to clarify.
Outline of Method: We begin by determining that certain elements are not zero-divisors, hence localizing to invert these elements does not change whether or not the ring is an integral domain. Then we use substitution to reduce the number of generators and eventually get to just $1$ homogeneous polynomial, which we show is irreducible by elementary methods.
The key is judicious localization. First note that the image of $x_1$ in $R$ is not a zero-divisor. To see this, suppose $x_1f\in I$ for some $f\in K[x_0,\ldots,x_3]$, so we have $f\in (I : (x_1))$. We want to show that $f\in I$, or equivalently that $(I:(x_1)) = I$. This can be done using Groebner bases. Since $K[x_0,\ldots,x_3]$ is an integral domain,
$$(I : (x_1)) = \frac{1}{x_1}(I\cap (x_1)) = \frac{1}{x_1}\left((tI+(1-t)(x_1))\cap K[x_0,\ldots,x_3]\right)$$
which is generated by the elements of a Groebner basis for $tI+(1-t)(x_1)$ that are independent of $t$ divided by $x_1$, and in fact this basis is a Groebner basis. Once you have this, you can compute a Groebner basis for $I$ and verify that the two ideals are the same by performing polynomial division on each generator of each ideal by all the generators of the other ideal, which should return $0$. It so happens that my desktop, which has mathematica and a fast processor, is a thousand miles away at the moment, so I will not do these calculations at this time (you however are free to, and I may do so later).
The following requires a couple of theorems concerning localization and polynomial rings. Having established that $x_1$ is not a zero-divisor, we let $S_1 = \{1,x_1,x_1^2,\ldots\}$ and note that $R$ is an integral domain iff $R_{S_1}$ is. We want to show that $x_2$ is not a zero divisor in order to eliminate another variable. We can certainly use the same technique as before to verify that it is not a zero-divisor in $R$, and if $fx_2=0$ in $R_{S_1}$ then we can clear denominators by multiplying by $x_1^k$ to get $(fx_1^k)x_2=0$ in $R$ and $fx_1^k=0$ iff $f=0$, so this verifies that it is not a zero-divisor in $R_{S_1}$. Again, I will not do these calculations here but will assume that they check out. Since $x_2$ is not a zero-divisor in $R_{S_1}$, we let $S_2 = \{1,x_2,x_2^2,\ldots\}$ and note that $R_{S_1}$ is an integral domain iff $R' = R_{S_1S_2}\cong (R_{S_1})_{S_2}$ is. Because $x_1,x_2$ are invertible in $R'$, we can let $x_0' = x_0+x_2^{-2}x_1^2x_3$ to get
$$R' \cong K[x_0,x_1,x_2,x_3]_{S_1S_2}/I_{S_1S_2} = (K[x_1,x_2])_{S_1S_2}[x_0',x_3]/I_{S_1S_2}$$
and using this we eliminate variables mercilessly
$$\begin{eqnarray*}I_{S_1S_2} &=& ((x_0'-x_2^{-2}x_1^2x_3)x_3-x_1x_2,
(x_0'-x_2^{-2}x_1^2x_3)^2x_2-x_1^3,
x_1x_3^2-x_2^3,
(x_0'-x_2^{-2}x_1^2x_3)x_2^2-x_1^2x_3)\\
&=& ((x_0'-x_2^{-2}x_1^2x_3)x_3-x_1x_2,
(x_0'-x_2^{-2}x_1^2x_3)^2x_2-x_1^3,
x_1x_3^2-x_2^3,
x_0'x_2^2)\\
&=& (x_2^{-2}x_1^2x_3^2-x_1x_2,
x_2^{-3}x_1^4x_3^2-x_1^3,
x_1x_3^2-x_2^3,
x_0')\\
&=& (x_1x_3^2-x_2^3,
x_1x_3^2-x_2^3,
x_1x_3^2-x_2^3,
x_0')\\
&=& (x_1x_3^2-x_2^3,x_0')\end{eqnarray*}$$
giving us a much nicer set of generators to work with. From this we see
$$R' \cong (K[x_1,x_2])_{S_1S_2}[x_3]/(x_1x_3^2-x_2^3)\cong \left(K[x_1,x_2,x_3]/(x_1x_3^2-x_2^3)\right)_{S_1S_2}$$
which is an integral domain iff $K[x_1,x_2,x_3]/(x_1x_3^2-x_2^3)$ is, so we want to show that $x_1x_3^2-x_2^3$ is prime in $K[x_1,x_2,x_3]$. Since $K[x_1,x_2,x_3]$ is a UFD, we need only show that $x_1x_3^2-x_2^3$ is irreducible. There are algorithms for this in general, but we can simplify this case considerably by observing that $x_1x_3^2-x_2^3$ is homogeneous. If it is not irreducible, then it can be factored into two homogeneous polynomials of degree $1$ and $2$. This can be easily checked, as if
$$\begin{eqnarray*}x_1x_3^2-x_2^3&=&(a_1x_1+a_2x_2+a_3x_3)(b_1x_1^2+b_2x_1x_2+b_3x_1x_3+b_4x_2^2+b_5x_2x_3+b_6x_3^2)\\
&=&(a_1b_6+a_3b_3)x_1x_3^2+a_2b_4x^3+a_2b_3x_1x_2x_3+a_2b_6x_2x_3^2+\text{ other junk}\end{eqnarray*}$$
which implies $a_2b_3=a_2b_6=0$ and $a_1b_6+a_3b_3\neq 0,a_2b_4\neq 0$, which is easily seen to be a contradiction. Hence $R'$ is an integral domain, so $I$ is prime.
Note on Radical Ideals: Proving an ideal radical is, in my experience, generally more difficult. I try to figure out the primary decomposition of the ideal, and then show that each primary ideal in the decomposition is in fact prime by something like the above method. However, if it is true that inverting a non-zero-divisor does not change whether or not a ring has only nilpotent zero-divisors, then maybe the method above can be modified to prove that ideals are radical.
Lemma: If $x\in R$ is not a zero-divisor and we let $S = \{1,x,x^2,\ldots\}$, then $R$ has only nilpotent zero-divisors iff $R_S$ does.
Proof: Suppose $R$ has some non-nilpotent $y,z$ such that $yz=0$. Since $x$ is not a zero-divisor, the homomorphism $R\to R_S$ is injective so $\frac{y^n}{1},\frac{z^n}{1}\neq 0$ for any $n$ so $\frac{y}{1},\frac{z}{1}$ are not nilpotent, yet $\frac{y}{1}\frac{z}{1}=\frac{0}{1}$ so $R_S$ has only non-nilpotent zero-divisors. Suppose now that $R_S$ has some non-nilpotent $\frac{y}{1},\frac{z}{1}$ such that $\frac{y}{1}\frac{z}{1}=\frac{0}{1}$. Clearly $y,z$ must be non-nilpotent, but $yz=0$ so $R$ has non-nilpotent zero-divisors.
Computations for $(I : (x_1))$: I forgot to mention that the Groebner basis must be computed with respect to a lexicographical ordering which places $t$ first. According to Sage, the following is a Groebner basis for $tI+(1-t)(x_1)$:
$$\{x_1^4 - x_0^2x_1x_2, tx_0^2x_2 - x_1^3, tx_0x_2^2 - x_1^2x_3, x_0x_1x_2^2 -x_1^3x_3, tx_2^3 - x_1x_3^2,$$
$$ x_1x_2^3 - x_1^2x_3^2, x_1^2x_2 - x_0x_1x_3,
tx_0x_3 - x_1x_2, tx_1 - x_1\}$$
which after eliminating the terms involving $t$ and dividing by $x_1$ gives us
$$(I : (x_1))=(x_1^3 - x_0^2x_2, x_0x_2^2 -x_1^2x_3, x_2^3 - x_1x_3^2, x_1x_2 - x_0x_3)$$
which is equal to $I$ because
$$\begin{eqnarray*} x_1^3 - x_0^2x_2 &=& -(x_0^2x_2-x_1^3)\\
x_0x_2^2-x_1^2x_3 &=& x_0x_2^2-x_1^2x_3\\
x_2^3 - x_1x_3^2 &=& -(x_1x_3^2-x_2^3)\\
x_1x_2 - x_0x_3 &=& x_1x_2 - x_0x_3
\end{eqnarray*}$$
hence $x_1$ is not a zero-divisor in $R$.
Computations for $(I : (x_2))$: Sage gives the following Groebner basis for $tI+(1-t)(x_2)$:
$$\{tx_1^3 - x_0^2x_2, x_1^3x_2 - x_0^2x_2^2, x_0x_2^3 - x_1^2x_2x_3, x_2^4 -x_1x_2x_3^2, tx_1^2x_3 - x_0x_2^2,$$
$$ tx_1x_3^2 - x_2^3, x_1x_2^2 - x_0x_2x_3,tx_0x_3 - x_1x_2, tx_2 - x_2\}$$
which after eliminating the terms involving $t$ and dividing by $x_2$ gives us
$$(I : (x_2))=(x_1^3 - x_0^2x_2, x_0x_2^2 - x_1^2x_3, x_2^3 -x_1x_3^2, x_1x_2 - x_0x_3)$$
which is identical to $(I : (x_1))=I$ hence $x_2$ is not a zero-divisor in $R$, which means it is not a zero-divisor in $R_{S_1}$.
This is actually fairly simple. Let $I$ be a radical ideal in $k[x_1, \dots, x_n]$, $k$ algebraically closed, then it corresponds to a unique algebraic set $V \subseteq \mathbb{A}^n(k)$, $V=Z(I)$.
So then $V = V_1 \cup \dots \cup V_k$, where $V_k$ are algebraic varieties corresponding uniquely to prime ideals $I_1, \dots, I_k$. Thus $Z(I) = Z(I_1) \cup \dots \cup Z(I_k)$. Now one can show that when taking the union of algebraic sets corresponding to arbitrary ideals, that $$Z(I_1)\cup \dots \cup Z(I_k) = Z(I_1 \cdots I_k) $$ where the product is the product of ideals. Also one has in general that $Z(J)= Z(\sqrt{J})$, so that $$Z(I_1)\cup \dots \cup Z(I_k) = Z(I_1 \cdots I_k)=Z\left(\sqrt{I_1 \cdots I_k}\right) $$ Now in general one only has that $I_1\cdots I_k \subseteq I_1 \cap \dots \cap I_k$, but taking the radical of the product allows us to conclude (see Lemma 1.7 here and then use induction): $$Z\left(\sqrt{I_1 \cdots I_k}\right)= Z\left( \sqrt{I_1 \cap \dots \cap I_k} \right) = Z\left( \sqrt{I_1} \cap \dots \cap \sqrt{I_k} \right). $$ However, since the $V_i$ are algebraic varieties, all of the $I_i$ were prime, thus radical, i.e. for all $i$ we have $\sqrt{I_i}=I_i$, therefore we have actually shown that: $$Z(I_1) \cup \dots \cup Z(I_k) =Z\left(\sqrt{I_1 \cdots I_k}\right) = Z(I_1 \cap \dots \cap I_k) $$ Now since $\sqrt{I_1 \cdots I_k} = I_1 \cap \dots \cap I_k$, the former obviously being radical since taking the radical of an ideal is idempotent operation, the latter must also be a radical ideal, and therefore our one-to-one correspondence between zero sets and radical ideals (from the Strong Nullstellensatz) applies, in particular we have that $$Z(I)=Z(I_1)\cup \dots \cup Z(I_k)=Z(I_1 \cap \dots \cap I_k) \iff I = I_1 \cap \dots \cap I_k.$$ Thus starting with an arbitrary radical ideal in $k[x_1,\dots,x_n]$, we have returned a (in fact unique) intersection of prime ideals equaling it, thereby proving this special case of the Lasker-Noether theorem. $\square$
Proof of the lemma (in case the link goes dead in the future)
The fact that $\sqrt{I \cdot J} \subseteq \sqrt{I \cap J}$ follows from the fact that $I \cdot J \subseteq I \cap J$.
If $ a \in \sqrt{I \cap J}$ then $a^n \in I \cap J$ for some $n\in \mathbb{N}$, thus both $a^n \in I$ and $a^n \in J$ are true, and thus $a \in \sqrt{I} \cap \sqrt{J}$. Thus $\sqrt{I \cap J} \subseteq \sqrt{I} \cap \sqrt{J}$.
Now if $a \in \sqrt{I} \cap \sqrt{J}$, then $a^n \in I$, $a^m \in J$, for some $m,n \in \mathbb{N}$. Therefore $a^{m+n} \in I \cdot J$, so $\sqrt{I} \cap \sqrt{J} \subseteq \sqrt{I \cdot J}$. $\square$
Best Answer
We have $$J=(x_2-x_1^2, \dots, x_n-x_1x_{n-1},x_1x_n-x_2x_{n-1},\dots, x_{n-2}x_n-x_{n-1}^2).$$
Moreover, $\Bbb C[x_1,\ldots, x_n]/J\simeq\mathbb C[x_1]$, so $J$ is a prime ideal. (Notice that $x_2$ is replaced by $x_1^2$, $x_3$ by $x_1^3$, and so on. Thus all the $2\times 2$ minors of the matrix $M$ with the first column removed become zero.)