```
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 with degree 2
poly = PolynomialFeatures(degree=2)
# 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)
```