:fn(cm, t) =sum over m((-i/h)*exp(-i(Em -En)t/h))Cm)