Write an \(\mathrm{R}\) function to compute the energy test statistic for a goodness-of-fit test of the simple hypothesis \(H_{0}: X \sim \operatorname{Beta}(\alpha, \beta)\), where \(\alpha, \beta\) are specified. How would you modify the function for the composite hypothesis?