Hi,
Say I want to compute the centralizer of a matrix.
sage: P.=PolynomialRing(QQ)
sage: A=matrix(P,2,2,[2,1,1,1])
sage: B=matrix(P,2,2,[a,b,c,d])
sage: C=A*B-B*A
sage: C
[-b + c -a + b + d]
[ a - c - d b - c]
sage: var('a b c d')
(a, b, c, d)
sage: solve([-b + c==0,-a + b + d==0,a - c
Hello,
I am trying to solve a system of matrix equations of the form:
transpose(B)*(Mi)*B==Ni
where i is in range(p) for some fixed p > 9, Mi and Ni have integer
entries, and all matrices are of size p by p. I want to solve for B
(it need not have integer entries).
I have tried a brute force
