第一矩阵
M1=structure(list(X1 = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L),
X0 = c(0.5, 0.5, 0, 0, 0, 0.5, 0.5, 1, -0.5, 1, -0.5, 0,
1, 0, -4, 0, 0.5, 1, 0, 0.5, 0, -2, 0.5, 0, 0), X0.1 = c(1.2,
0, 0, 0.5, 0, 0.5, 0, 0, 0, 0, 0,