@@ -317,7 +317,7 @@ def lqe(A, G, C, QN, RN, NN=None):
     A, G, C = np.array(A, ndmin=2), np.array(G, ndmin=2), np.array(C, ndmin=2)
     QN, RN = np.array(QN, ndmin=2), np.array(RN, ndmin=2)
     P, E, LT = care(A.T, C.T, np.dot(np.dot(G, QN), G.T), RN)
-    return _ssmatrix(LT.T), _ssmatrix(P), _ssmatrix(E)
+    return _ssmatrix(LT.T), _ssmatrix(P), E
 
 
 # Contributed by Roberto Bucher <roberto.bucher@supsi.ch>
 
@@ -302,7 +302,7 @@ def test_LQR_3args(self):
     def check_LQE(self, L, P, poles, G, QN, RN):
         P_expected = np.array(np.sqrt(G*QN*G * RN))
         L_expected = P_expected / RN
-        poles_expected = np.array([-L_expected], ndmin=2)
+        poles_expected = np.array([-L_expected])
         np.testing.assert_array_almost_equal(P, P_expected)
         np.testing.assert_array_almost_equal(L, L_expected)
         np.testing.assert_array_almost_equal(poles, poles_expected)