make_quad_form_matrix function

Represent a variance estimator as a quadratic form