A<-matrix(c(1842,25440,3167,8522,359803,114617,4023,85978,189892), nrow=3, ncol=3, byrow=T)
A
d<-matrix(c(7838,167051,424756), nrow=3)
d
X<-matrix(c(38287,649993,704649), nrow=3)
X

IC<-matrix(data=NA, nrow=3, ncol=3, byrow=T)
for(i in 1:3) {
  for(j in 1:3) {
    IC[i,j]<-(A[i,j]/X[j])
  }
}
IC

I3<-matrix(c(1,0,0,0,1,0,0,0,1), nrow=3)
I3

TM<-I3-IC
TM

PIC<-solve(TM)
PIC

dd<-1.5*d
Xd<-PIC%*%dd
Xd

de<-dd-d
de

ie<-Xd-X-de
ie

(rp<-rowSums(PIC))
(cp<-colSums(PIC))
(m<-mean(rp))

(senstivity<-rp/m)
(impact<-cp/m)