numerDecom.lib
LIB "numerDecom.lib"; ring r=(complex,16,I),(x,y,z),dp; poly f1=(x2+y2+z2-6)*(x-y)*(x-1); poly f2=(x2+y2+z2-6)*(x-z)*(y-2); poly f3=(x2+y2+z2-6)*(x-y)*(x-z)*(z-3); ideal J=f1,f2,f3; list w=list(0.5372775295412116,-0.7105339291010922,-2.2817700129167831); def D=JuReTopDim(J,w,2,2); setring D; t;