Agree. Given the power ordering, `numpy.polynomial.polynomial.polyvander` is a better fit but will need a transpose (as does np.vander).