larstiq

try:
    import numpy
except ImportError:
    import numpypy as numpy



data = [['a', 0]]
data = numpy.array(data)

larstiq

gbTrace = 3
load "realroots.m2"
W = QQ[a, b, c, d, z]

I = ideal(2*a^3+a^2*b+a*b^2+2*b^3+6*a^2*c+2*a*b*c+b^2*c+6*a*c^2+b*c^2+2*c^3+a^2*d+2*a*b*d+6*b^2*d+2*a*
      c*d+2*b*c*d+c^2*d+a*d^2+6*b*d^2+c*d^2+2*d^3+(2/3)*a^2+(5/3)*a*b+b^2+(4/3)*a*c+(5/3)*b*c+(2/3)*c^2
      +(5/3)*a*d+2*b*d+(5/3)*c*d+d^2+a+(2/5)*b+c+(2/5)*d+3/5,2*a^3+a^2*b+a*b^2+2*b^3-6*a^2*c-2*a*b*c-b^2
      *c+6*a*c^2+b*c^2-2*c^3-a^2*d-2*a*b*d-6*b^2*d+2*a*c*d+2*b*c*d-c^2*d+a*d^2+6*b*d^2-c*d^2-2*d^3+(2/3
      )*a^2+(5/3)*a*b+b^2-(4/3)*a*c-(5/3)*b*c+(2/3)*c^2-(5/3)*a*d-2*b*d+(5/3)*c*d+d^2+a+(2/5)*b-c-(2/5)*
      d+3/5,2*a^3+a^2*b+a*b^2+2*b^3-a^2*c-2*a*b*c-6*b^2*c+a*c^2+6*b*c^2-2*c^3+6*a^2*d+2*a*b*d+b^2*d-2*a*