bd_ML function

Maximization of the loglikelihood under the diversity-independent, possibly time-dependent diversification model