Create a matrix that represents the following transformation
T_center ((x, y, z)ᵀ) = (x - Cₓz, y - Cᵧz, z),
where Cₓ, Cᵧ are the coordinates of the center found in the previous part. i. e. Translation in the y-direction by one unit. You can test that you have the correct matrix by applying the transformation (multiplying the transformation matrix) to each column in the matrix of homogeneous coordinates and then plotting the new points using the plot points() function from above. Use a for-loop to apply this transformation.