Обратная матрица вероятностей с функцией решения

Question

Я хочу взять обратную матрицу i.p с solve функция. это i.p матрица является результатом вычитания единичной матрицы p, Вот код, чтобы получить i.p:

p<-rbind(c(0,0,0,0,0,0,3), c(0,0,7,9,0,0,0), c(3,6,0,75,0,0,0), c(0,0,20,0,1,35,34), c(0,0,1,0,0,0,0), c(0,0,10,1,0,0,27), c(0,3,31,2,0,0,0))
dev<-c(sum(0,0,0,0,0,0,3), sum(0,0,7,9,0,0,0), sum(3,6,0,75,0,0,0), sum(0,0,20,0,1,35,34), sum(0,0,1,0,0,0,0), sum(0,0,10,1,0,0,27), sum(0,3,31,2,0,0,0))
p<-p/dev
p[7,]<-c(0,3,31,2,0,0,0)/sum(0,3,31,2,0,0,0, 30)
i<-diag(1, 7, 7)    #make a 7X7 identity matrix
i.p<-i-p         #identity matrix substract probability matrix  I-P
i.p.s<-solve(i.p, i)   #inverse of I-P

Вот p, i.p, а также i.p.s матрицы:

Это обратное i.p.s умножает матрицу i.p должна быть единичная матрица.

i.p %*% i.p.s

Однако этот множитель, рассчитанный с помощью приведенного выше кода, явно не был идентичной матрицей. Это почему?