from sklearn.preprocessing import PolynomialFeatures
import numpy as np

# Create a sample dataset with two features
X = np.array([[1, 2], [3, 4], [5, 6]])

# Initialize the PolynomialFeatures object
poly = PolynomialFeatures()

# Transform the original features to polynomial features
X_poly = poly.fit_transform(X)

# Print the original features and their interactions
print('Original Features:\n', X)
print('Polynomial Features:\n', X_poly)