For a given linear regression model, generate a sample of size n= 100 with parameter values ? = 5, ? = 1,? = 3, ?^(2) = 0.2, and constants xi = i/100, z i= (100?i)/100, i= 1,...,100. Write a function in R that finds the maximum likelihood estimates of the model parameters ?, ?, ??R1 and ?^(2) ? (0,10).

Yi = at (Bxcityzi) + BEi, ci ~N(0, 02), i = 1, ..., n