Tuning the Polynomial Feature for Logistic Regression in Python

Question

If I want to incorporate degree two polynomials into my logistic model (which has two predictor variables), such as I have tried below:

df_poly = df[[\'Y\',\'