f(x,y,z) = (y*z,x*z,x*y)
var('a1 a2 a3 a4 a5 a6 a7 a8 a9') A = Matrix([[a1,a2,a3],[a4,a5,a6],[a7,a8,a9]]); show(A)
(a1, a2, a3, a4, a5, a6, a7, a8, a9)
a1a4a7a2a5a8a3a6a9
var('b1 b2 b3 b4 b5 b6 b7 b8 b9') B = Matrix([[b1,b2,b3],[b4,b5,b6],[b7,b8,b9]]); show(B)
(b1, b2, b3, b4, b5, b6, b7, b8, b9)
b1b4b7b2b5b8b3b6b9
var('x11 x12 y11 y12') X1 = Matrix([[x11],[x12], [1]]) Y1 = Matrix([[y11],[y12], [1]])
(x11, x12, y11, y12)
(A*X1)[2]
(a7*x11 + a8*x12 + a9)
E = B*f((A*X1)[0][0],(A*X1)[1][0],(A*X1)[2][0]) show(E[0].simplify_full()) show(E[1].simplify_full()) show(E[2].simplify_full())
a6a9b1+a3a9b2+a3a6b3+(a4a7b1+a1a7b2+a1a4b3)x112+(a5a8b1+a2a8b2+a2a5b3)x122+((a6a7+a4a9)b1+(a3a7+a1a9)b2+(a3a4+a1a6)b3)x11+((a6a8+a5a9)b1+(a3a8+a2a9)b2+(a3a5+a2a6)b3+((a5a7+a4a8)b1+(a2a7+a1a8)b2+(a2a4+a1a5)b3)x11)x12
a6a9b4+a3a9b5+a3a6b6+(a4a7b4+a1a7b5+a1a4b6)x112+(a5a8b4+a2a8b5+a2a5b6)x122+((a6a7+a4a9)b4+(a3a7+a1a9)b5+(a3a4+a1a6)b6)x11+((a6a8+a5a9)b4+(a3a8+a2a9)b5+(a3a5+a2a6)b6+((a5a7+a4a8)b4+(a2a7+a1a8)b5+(a2a4+a1a5)b6)x11)x12
a6a9b7+a3a9b8+a3a6b9+(a4a7b7+a1a7b8+a1a4b9)x112+(a5a8b7+a2a8b8+a2a5b9)x122+((a6a7+a4a9)b7+(a3a7+a1a9)b8+(a3a4+a1a6)b9)x11+((a6a8+a5a9)b7+(a3a8+a2a9)b8+(a3a5+a2a6)b9+((a5a7+a4a8)b7+(a2a7+a1a8)b8+(a2a4+a1a5)b9)x11)x12
var('k') k*Y1
[k*y11]
[k*y12]
[ k]