x <- c(0.49487,0.52802,0.28667,0.62058,0.14704,0.18519,
       0.17889,0.14554,0.29480,0.46317)
xi <- max(x)
gamma <- length(x)
p <- 0.9
xtop <- xi/(1 - p)^(1/gamma)
cat("HDR is",c(xi,xtop),"\n")