Estoy usando la librería sympy <code>import sympy as sp sp.init_printing() x1, x2, x3 = sp.symbols(&quot;x1 x2 x3&quot;) y1, y2, y3 = sp....

Estoy usando la librería sympy

import sympy as sp
sp.init_printing()
x1, x2, x3 = sp.symbols("x1 x2 x3")
y1, y2, y3 = sp.symbols("y1 y2 y3")
z1, z2, z3 = sp.symbols("z1 z2 z3")
G1, G2, G3 = sp.symbols("G1 G2 G3")
A = sp.Matrix([[x1, y1, z1], [x2, y2, z2], [x3, y3, z3]])
A_T = A.T
T_Rot = sp.Matrix([[G1, 0, 0], [0, G2, 0], [0, 0, G3]])
T = A_T * T_Rot * A
psi = sp.symbols("psi")
T1 = T.subs(G1, sp.cos(psi))
T1 = T1.subs(G2, -(1/2)*sp.cos(psi)-(sp.sqrt(3)/2)*sp.sin(psi))
T1 = T1.subs(G3, -(1/2)*sp.cos(psi)+(sp.sqrt(3)/2)*sp.sin(psi))

Cuando intento agrupar los terminos de cada elemento de la matriz por cos(phi) y sin(phi)

for i in range(9):
  T1[i] = sp.collect(T1[i],[sp.cos(psi), sp.sin(psi)])
T1

No obtengo ningun cambio, siguen agrupados por las variables x1, x2 y x3

Me gustaria me pudieron colaborar en darme alguna idea para poderlos agrupar por cos(phi) y sin(phi)

Gracias