from sage.rings.rational_field import QQ
from sage.rings.polynomial.polynomial_ring_constructor import PolynomialRing
Z = PolynomialRing(QQ, "z")
z = Z.gen()
coefficients = [1, 14, 87, 148, 87, 14, 1]
sum(QQ(c) * z**i for i, c in enumerate(coefficients)) # h^*_{B_4}(z)The $h^*$-polynomials are transcribed from OEIS A259473 [5].