Procedure from library sing.lib (see sing_lib).
sing.lib
slocus(i); i ideal
ideal of singular locus of i
Example:
LIB "sing.lib"; ring r = 0,(u,v,w,x,y,z),dp; ideal i = wx,wy,wz,vx,vy,vz,ux,uy,uz,y3-x2;; slocus(i); → _[1]=x → _[2]=w → _[3]=v → _[4]=u → _[5]=y2