kk=ZZ
R=kk[x0,x1,x2,M]
I=ideal(M*x0*x1*x2-(x0^2+x1^2+x2^2))
J=saturate(I,x0*x1*x2)
S=R/J
d=x0*x1*x2
D=ideal(x0*x1*x2,(x1*x2*M-x0)*x1*x2,x0*(x0*x2*M-x1)*x2,x0*x1*(x0*x1*M-x2))
n=1
D'=ideal(d^n)
M=saturate(D',D)
D'==M