How to construct a polynomial for numerical data?

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.