M31 = M1(1:10001, 2);
M37 = M7(1:10001, 2);