powered by
calculates the p matrix
get_p_matrix(branch_length, eig = phangorn::edQt(), rate = 1)
p matrix
branch length
eigen object
rate