R: Create a function equivalent to lm()

Question

I need to create a function equivalent to lm() in R.

Model is: Yt = B0 + B1*Zt + Ut

My data is:

set.seed(122)
n <- 100
v <- rnorm(n ,mean=