xbar<-c(95.52,164.38,55.69,93.39,17.98,31.13)

  s<-matrix(c(3266.46,  1343.97, 731.54,  1175.50, 162.68, 238.37,  1343.97,  721.91, 324.25, 537.35, 80.17, 117.73,  731.54,  324.25,  179.28, 281.17, 39.15, 56.8,  1175.50,  537.35,  281.17,  474.98,  63.73, 94.85,  162.68,  80.17,  39.15,  63.73,  9.95, 13.88, 238.37,  117.73,  56.80,  94.85,  13.88,  21.26),6,6,byrow=T)

# part(a)

n<-61

p<-6

c<-sqrt(p*(n-1)/n/(n-p)*qf(0.95,p,n-p))

for (i in 1:6) {

  lower<-xbar-c*sqrt(s[i,i])

  upper<-xbar+c*sqrt(s[i,i])

  cat("(",lower,",",upper,")","\n")

}