Given the hamiltonian of an open quantum system modeled by the Feynman-Vernon model or Leggett-Caldeira model, and a spectral density, temperature and initial system density matrix, the program calculates the system density matrix at future points in time.