> // LOCAL DIMENSION
> ring R = 0, (x,y,z), ls;
> poly f = z3-2zx2+y2;
> ideal I1 = f;
> ideal I2 = x,z;
> ideal I = I1*I2;

> // find the local dimension of V(I) at (0,1,0)
> ideal J = I;
> J = subst(J, y, y+1); //move (0,1,0) to the origin
> dim(std(J));
1
> std(J);
_[1]=z
_[2]=x


> // find the local dimension at the origin
> dim(std(I));
2

> // you can determine Tjurina number either using "sing.lib"
> // or manually 
> ideal J = jacob(f);
// ** redefining J **
> J;
J[1]=-4xz
J[2]=2y
J[3]=3z2-2x2
> J = J,f;
> kbase(std(J));
_[1]=z
_[2]=x2
_[3]=x
_[4]=1
> // Tjurina number of f at the origin is 4.