You want to compute the critical values of a
polynomial mappingfrom C^3-> C.
The primary decomposition is not what you really need to compute.
You have just to eliminate the variables x,y,z from the ideal.
Code:
> ring r=0,(x,y,z,t),dp;
> poly f = 3*x^2-3*x*y^2+y^3-3*y+3*z^2;
> ideal i = jacob(f),f-t;
> eliminate(i,xyz);
_[1]=64t3+176t2-24t-665
This polynomial is irreducible over Q.
There are three isolated singularities each with Milnor number 1
Code:
> ring rdp =0,(x,y,z),dp;
> poly f = 3*x^2-3*x*y^2+y^3-3*y+3*z^2;
> vdim(std(jacob(f)));
3
The coordinates (x,y,z) expressed in terms by t
can be read off from the primdec result.
Code:
> setring r;
> primdecSY(i);
[1]:
[1]:
_[1]=64t3+176t2-24t-665
_[2]=z
_[3]=575y-16t2+248t+80
_[4]=1150x-144t2-68t+145