Add integrate method for hierarchy initial cond
Add integrate method for hierarchy initial cond

Might be useful to be able to supply the hierarchy initial
conditions as a density matrix as well.

CQuIC/pysme