-e1 = computeError(data, p_nz_mf, H)
-e2 = computeError(data, p_nz_mfl, H)
-e3 = computeError(data, p_a, H)
-e4 = computeError(data, p_p, H)
+e1 = computeError(data, p1, H)
+e2 = computeError(data, p2, H)
+e2 = computeError(data, p3, H)
+e3 = computeError(data, p4, H)
+e4 = computeError(data, p5, H)