  An analytically derived 'integral operator' approach is introduced to
estimate the expectation value of a quantum operator for an evolving state
weighted with an exponential function. This allows to compute quantities useful
in Nuclear Magnetic Resonance or Muon Spectroscopy experiment simulations with
greater speed and precision than the standard numerical integration approach.
