Creating new function in R equivalent to the lm() for finding coefficients b0 and b1

Question

I have a task to estimate the model to obtain b0 and b1 on the generated data (yt, zt and wt) using my created function an