I have the following question...
Consider the ideal $J:= (x^2y-x^2y^2,\ x^2z-z^2yx,\ x^2+xz) \subset \mathbb{Q}[x,y,z]$
Is $x \in J?$ Is $x \in \sqrt{J} $?
I know finding if $x$ is in the radical of $J$ I need to compute the Gröbner basis and work it out from there... But I don't know the method to use to show $x \in J $
Groebner basis for $J$ in lexicographic ordering $x \prec y \prec z$ is $$ J := \left(x z^2, x y z(y-1), x(x+z) \right) $$ By applying the multivariate division algorithm relative to the same ordering you find $$ x = 0 \cdot p_1 + 0 \cdot p_2 + 0 \cdot p_3 + x $$ therefore $x\not\in J$, but $$ x^3 = 1 \cdot p_1 + 0 \cdot p_2 + (x-z) \cdot p_3 + 0 $$ hence $x \in \operatorname{Rad}(J)$.
Here are relevant Mathematica commands: