fit_musse function

Fit a discrete-state-dependent diversification model via maximum-likelihood.