Simulates node depths in the birth-death model, conditionned by the number of tips and the age of the root (using the expression in "The shape and probability of reconstructed phylogenies", by Amaury Lambert and Tanja Stadler (Theoretical Population Biology, 2013))
yule_lengths(N, b, d, tmax)
Tip number
Birth rate
Death rate
Maximum age of the root
Lambert, Amaury, and Tanja Stadler. Birth-death models and coalescent point processes: The shape and probability of reconstructed phylogenies. Theoretical population biology 90 (2013): 113-128.