rn_cs_gen function

Decompose the $\mathbf{G}$-matrix from a character-state model