added function to simulate heterochronous data under the coalescent