simUltra function

Simulate ultrametric birth and death branching trees