sim_t_tworegime function

Recursive simulation (root-to-tip) of two-regime models