Computes the parameter lambda_epsilon for n in 1:N
lambda_N(epsilon, beta, N, upper = 10000)Minimum size of unsampled splits (see appendix 1)
The beta parameter quantifying the balance of the tree
Tip number
Upper bound from the integral (see lambda_n expression, default to 10000)
\(\lambda_\epsilon = 2 \int_(\epsilon,\infty) exp(-(\beta+n+1)x)(1-exp(-x))^\beta dx\)
Maliet O., Gascuel F., Lambert A. (2018) Ranked tree shapes, non-random extinctions and the loss of phylogenetic diversity, bioRxiv 224295, doi: https://doi.org/10.1101/224295