import numpy as np

def procrustes(X, Y):
  X = X - X.mean(axis=0)
  Y = Y - Y.mean(axis=0)
  (U, S, Vt) = np.linalg.svd(np.dot(Y.T,X))
  R = np.dot(U, Vt)
  Y = np.array(np.matrix(Y) * np.matrix(R))
  return X, Y