Question: Primary Decomposition over C

How can i specify the ground field C= complex numbers and calculate a primary decomposition for an ideal defined over C? For example

with(PolynomialIdeals):
J:= <x^2 +1 >;
PrimaryDecomposition(J);

should return x+I, x-I, but instead returns x^2+1

Please Wait...