reconstruct_past_diversification function

Reconstruct past diversification dynamics from a diversity time series.