python stats models - quadratischer Term in der Regression

Ich habe die folgende lineare Regression:

import statsmodels.formula.api as sm

model = sm.ols(formula = 'a ~ b + c', data = data).fit()

Ich möchte in diesem Modell einen quadratischen Term für b hinzufügen.

Gibt es eine einfache Möglichkeit, dies mit statsmodels.ols zu tun? Gibt es ein besseres Paket, das ich verwenden sollte, um dies zu erreichen?