sage: B = matrix([[0,1,2],[-1,0,3],[-2,-3,0]])
sage: q = ZZ['q'].gen()
sage: R.<x,y,z> = algebras.qCommutingPolynomials(q, B)
sage: R.formal_series_ring()
Lazy completion of q-commuting polynomial ring in x, y, z over Univariate Polynomial Ring in q over Integer Ring with matrix:
[ 0 1 2]
[-1 0 3]
[-2 -3 0]
One could maybe make good use of
?