F=[x^2+y^2+z^2-1, x*y*z-1,x^3-y*z];
V=[x,y,z];
G=gr(F,V,0);
print(zero_dim(G,V,0))$
N=length(dp_mbase(map(dp_ptod,G,V)))$
print("Number of solution is ",0);
print(N);