How to construct a polynomial for numerical data?

[Alalalalaki]

For example I have some variables "x1, x2, x3, ...", and want all polynomials with a simple choice on the degree. The only solution I know so far is by using preprocessing.PolynomialFeatures in sklearn. But I can't see why it cannot be done simply in patsy.

Thanks.