set.seed(100)
x <- rnorm(1000, 0, 1)
y <- x + rnorm(1000)+1

Linear_Regression<-lm(y~x)

plot(x,y,col="cornflowerblue")
abline(coef(Linear_Regression),lwd=3,col="orange")

Linear_Regression