Is there a way to find the symbolic matrix associated to an expended quadratic form?
Suppose I have a complex quadratic form
Q(px,py)=p2yy2α2−2pxpyxyαβ+p2xx2β2−p2yy2α−p2xx2β
Of course all parameters has been declared as variables, that is
var('p_x p_y x y \alpha \beta')
(I use unicode α and β)
Is there a way to find the matrix A such that
Q(p)=p⊤Ap ?