pexact <- function(x,mu0) {
    n <- length(x); s0 <- sd(x)
    t0 <- (mean(x)-mu0)/(s0/sqrt(n))
    2*(1-pt(abs(t0),n-1))
}