qmatrix function

Transition intensity matrix