code: def poly(self, return_coeffs=False):
          
          p = self.bpoints()
          coeffs = (p[0] - 2*p[1] + p[2], 2*(p[1] - p[0]), p[0])
          if return_coeffs:
              return coeffs
          else:
              return np.poly1d(coeffs)
docs: returns the quadratic as a Polynomial object.