Add polygrid2d() for consistency

Add a wrapper for `polygrid2d()` in `numpy.polynomial.polynomial`.
Improve doc