function(x, y, m)
{
	n <- length(x)
	b0 <- c(1:m)
	b1 <- c(1:m)
	ser <- c(1:m)
	for(i in 1:m) {
		a <- sample(c(1:n), n, 
			replace = T)
		sx <- x[a,  ]
		sy <- y[a,  ]
		b1[i] <- (cor(sx, sy) * std(
			sy))/std(sx)
		b0[i] <- mean(sy) - b1[i] * 
			mean(sx)
		ser[i] <- std(sy) * sqrt(1 - 
			cor(sx, sy)^2)
	}
	res <- cbind(b0, b1, ser)
	res
}
