R code to calculate maximum likelihood branch lengths for a star-like tree with a gamma rate distribution.