Plotting mechanism for polynomial zonotopes

[mforets]

In principle the mechanism would be a replica of overapproximate(X::TaylorModelReachSet, Zonotope; kwargs...), but without the time domain, ie. the relevant keywords are dom (to choose a subdomain of [-1, 1]^n, and nsdiv, partition to partition such domain. The result is either a zonotope or a union of zonotopes. We can have a similar method with Hyperrectangle as output.

cc: @lucaferranti

[mforets]

This will be basically implemented in #2972; it will only remain to add a plot recipe here, which first overapproximates the PZ with a union of Z and then calls the appropriate plot recipe.