Procedure from library sing.lib (see sing_lib).
sing.lib
tjurina(id); id=ideal or poly
id=ICIS (isolated complete intersection singularity)
int = Tjurina number of id
Tjurina number = -1 implies that id is not an ICIS
Example:
LIB "sing.lib"; ring r=32003,(x,y,z),(c,ds); ideal j=x2+y2+z2,x2+2y2+3z2; tjurina(j); → 5