import numpy as np
from sklearn.preprocessing import PolynomialFeatures
X = np.array([2, 3])
poly = PolynomialFeatures(3)
Y = poly.fit_transform(X)
print(Y)
print(poly.powers_)
sklearn: how to get coefficients of polynomial features