diff --git a/p2/stats/stats.go b/p2/stats/stats.go
index 7813547..f53a283 100644
--- a/p2/stats/stats.go
+++ b/p2/stats/stats.go
@@ -123,6 +123,7 @@ func ImpulseFunction(ruu []float64, ruy []float64) ([]float64, error) {
 
 	// g.Mul(invRuuMatrix, ruyVect)
 	g.Mul(ruyVect, invRuuMatrix)
+	// g.Product(ruyVect, invRuuMatrix)
 
 	// calculate h - the impulse function.
 	h, err := estimateImpulseFunc(g, lruu)