make_quad_form_matrix function

Get Variance Estimator's Quadratic Form Matrix