docs: def fit_quadratic(X, y):
      """
      Uses OLS with Polynomial order 2.
      """
      model = make_pipeline(
          PolynomialFeatures(2), linear_model.LinearRegression()
      )
      model.fit(X, y)
      return model