- sage: x = sum( i*J.gens()[i] for i in range(6) )
- sage: basis = tuple( x^k for k in range(J.rank()) )
- sage: K = FiniteDimensionalEuclideanJordanSubalgebra(J,basis)
- sage: [ K(x^k) for k in range(J.rank()) ]
- [f0, f1, f2]
+ sage: X = matrix(QQ, [ [0,0,1],
+ ....: [0,1,0],
+ ....: [1,0,0] ])
+ sage: x = J(X)
+ sage: basis = ( x, x^2 ) # x^2 is the identity matrix
+ sage: K = FiniteDimensionalEuclideanJordanSubalgebra(J, basis)
+ sage: K(J.one())
+ f1
+ sage: K(J.one() + x)
+ f0 + f1