function(x, mu, m)
{
	n <- length(x)
	avg <- mean(x)
	sig <- std(x)
	tn <- ((avg - mu) * sqrt(n))/sig
	tsx <- c(1:m)
	for(i in 1:m) {
		xs <- sample(x, n, replace
			 = T)
		avgs <- mean(xs)
		sigs <- std(xs)
		tsx[i] <- ((avgs - avg) * 
			sqrt(n))/sigs
	}
	prop <- 0 * jay(m, 1)
	for(i in 1:m) {
		if(abs(tsx[i]) >= abs(tn))
			prop[i] <- 1
	}
	pvalue <- sum(prop)/m
	res <- list(avg, sig, tn, pvalue)
	names(res) <- c("S.Mean", "S.STD", 
		"t-value", "P-Value")
	res
}
