code: def expected(self, observable, beta=1e5):
          
          return expected_value(observable,
                           self.eig_energies,
                           self.eig_states,
                           beta)
docs: Wrapper to the expected_value function to fix the eigenbasis