function(x, y)
{
	n <- length(x)
	x <- cbind(jay(n, 1), x)
	s <- t(x) %*% x
	b <- inv(s) %*% t(x) %*% y
	yhat <- x %*% b
	vyx <- sum((y - yhat)^2)/(n - 2)
	ser <- sqrt(vyx)
	res <- list(b, ser)
	names(res) <- c("LSE", "SE")
	res
}
