code: def density(self):
          
          h = self.histogram.astype(np.float)
          bindifs = np.array(np.diff(self.bin_edges), float)
          return h / (bindifs * self.n)
docs: Gives emprical PDF, like np.histogram(...., density=True)