Use zsolve to compute Hilbert basis of a cone