Symbolic matrix multiplication to polynomial

Collect the array variables immediately

julia> x,y = collect.((x,y))
(Num[x[1], x[2], x[3]], Num[y[1], y[2], y[3]])

julia> y'A*x
(14y[1] - 3y[2] - 2y[3])*x[3] + (8y[3] + y[1] - 5y[2])*x[1] + (2y[1] + 2y[2] + 19y[3])*x[2]
1 Like