density matrix has output now.  looks like it has negative elements on the diagonal...
kind of working now?
added output for populations, can see they are broken
still broken
found error in f
grr sign error
found sign error, now it seems to work for the general (k-b) two-level system