code: def fit(self, order=3):
          
          x = self.volumes**(-2./3.)
          self.eos_params = np.polyfit(x, self.energies, order)
          self._set_params()
docs: Overriden since this eos works with volume**(2/3) instead of volume.