source("hdr.txt")
x=c(9,14,15,15,16,18,18,19,19,20,21,22,22,24,24,26,27,29,30,32)
n <- length(x)
S <- sum((x-mean(x))^2)
cat("Values of chi^2 corresponding to HDR for log chi^2",
    hchisq(0.95,n),"\n")
cat("Interval for phi                                  ",
    sort(S/hchisq(0.95,n)),"\n")
cat("HDR for inverse chi^2                             ",
    hchisq(0.95,n,log=FALSE,inverse=TRUE),"\n")
cat("Interval for phi                                  ",
    S*hchisq(0.95,n,log=FALSE,inverse=TRUE),"\n")