X-Git-Url: http://gitweb.michael.orlitzky.com/?a=blobdiff_plain;f=mjo%2Feja%2Feuclidean_jordan_algebra.py;h=1426d5e16be4b6acc7c68a4494e2c6f1c4d61819;hb=4156fccce1265f500fe432b0f5567e43fbbc23d6;hp=3a98324cbb7c9a7f18a86677a53785a81988d28b;hpb=2edbb6042b34b4c366a52fa9be0d03bada22969d;p=sage.d.git

diff --git a/mjo/eja/euclidean_jordan_algebra.py b/mjo/eja/euclidean_jordan_algebra.py
index 3a98324..1426d5e 100644
--- a/mjo/eja/euclidean_jordan_algebra.py
+++ b/mjo/eja/euclidean_jordan_algebra.py
@@ -124,6 +124,18 @@ class FiniteDimensionalEuclideanJordanAlgebra(FiniteDimensionalAlgebra):
                 sage: (x*x)*(x*x*x) == x^5
                 True
 
+            We also know that powers operator-commute (Koecher, Chapter
+            III, Corollary 1)::
+
+                sage: set_random_seed()
+                sage: x = random_eja().random_element()
+                sage: m = ZZ.random_element(0,10)
+                sage: n = ZZ.random_element(0,10)
+                sage: Lxm = (x^m).matrix()
+                sage: Lxn = (x^n).matrix()
+                sage: Lxm*Lxn == Lxn*Lxm
+                True
+
             """
             A = self.parent()
             if n == 0: