diff --git a/adapt/feature_based/_tca.py b/adapt/feature_based/_tca.py
index 45847ec..73a72b5 100644
--- a/adapt/feature_based/_tca.py
+++ b/adapt/feature_based/_tca.py
@@ -123,7 +123,7 @@ def fit_transform(self, Xs, Xt, **kwargs):
         H = np.eye(n+m) - 1/(n+m) * np.ones((n+m, n+m))
         
         # Compute solution
-        a = np.eye(n+m) + self.mu * K.dot(L.dot(K))
+        a = K.dot(L.dot(K)) + self.mu * np.eye(n+m)
         b = K.dot(H.dot(K))
         sol = linalg.lstsq(a, b)[0]
         
@@ -166,4 +166,4 @@ def transform(self, X, domain="tgt"):
 
         K = np.concatenate((Kss, Kst), axis=1)
         
-        return K.dot(self.vectors_)
\ No newline at end of file
+        return K.dot(self.vectors_)