normal.lib
LIB "normal.lib"; ring r = 0,(x,y,z),dp; poly f = z5+y4+x3+xyz; ideal i = jacob(f),f; i=std(i); locAtZero(i); ==> 1 i= std(i*ideal(x-1,y,z)); locAtZero(i); ==> 0