seqtrate function

Compute transition rates between states