implement the normal equation method for linear regression in python

# Sub-task 1) # Obtain model parameters using normal equation Get matrix X, adding a column of 1s to X, and vector 1 X = np. ones-Like (data, train) X:, 1:31-data_train, 0:21 y = data-train:,-1 # START TODO # Hint: theta = (XT X)^{-1} theta_method1- X^T y print('theta_methodl:', theta_method1)