normal.lib
LIB "normal.lib"; ring R = 0, (x,y), dp; ideal I1 = x2 - y3; isNormal(I1); ==> 0 ideal I2 = x - y3; isNormal(I2); ==> 1