def plot(self):
        """Plot fitted EOS curve."""
        import matplotlib.pyplot as plt

        ep = self.get_parameters()
        vols = self._volumes
        volume_points = np.linspace(min(vols), max(vols), 201)
        fig, ax = plt.subplots()
        ax.plot(volume_points, self._eos(volume_points, *ep), "r-")
        ax.plot(vols, self._energies, "bo", markersize=4)
        return plt