solve(crossprod(X), crossprod(X, Y))