Допишите код из прошлой задачи result = np.linalg.inv(X.T @ X) так, чтобы в переменной result был результат операции (X^T * X)^(−1) * X^T. Выведите форму массива, который получился.

# необходимые импорты
import numpy as np
# считывание npy-файла
X = np.load('matrix_X_train_preprocessed.npy')
# матричное умножение на транспонированную матрицу
result = np.linalg.inv(X.T @ X) @ X.T
# вывод формы результата
print(result.shape)