quadratic_form function

Computes quadratic form x' A x