Can someone help me with R programming for this problem?

Consider the following data: X 3.04 4.55 7.47 6.33 5.00 5.20 9.66 11.11 8.34 7.10 y 10.48 1.76 -19.38 50 -.63 0.08 -16.88 -27.66 -11.95 -9.05 Z 21.34 15.14 -.10 19.26 15.60 16.14 7.31 .51 9.81 10.06 Consider the data in the previous problem. Use the matrix form of the least-squares coefficients to build a linear model for z using x and y as predictor variables. To do this, let the data matrix X have three columns: the first is a column of ten Is, the second has all the x-values, and the third column has all the y- values. Carry out the process using matrix multiplication in R: B = (XTX) -1XTY, then verify your coefficients are correct with Im(z ~ x + y)