Q<-matrix(c(0,1/3,0,0,1/3,0,0,0, 1/2,0,1/2,0,0,0,0,0, 0,1/3,0,1/3,0,0,0,0, 0,0,1/2,0,0,0,0,1/2, 1/2,0,0,0,0,1/2,0,0, 0,0,0,0,1/3,0,1/3,0, 0,0,0,0,0,1/2,0,1/2, 0,0,0,1/3,0,0,1/3,0), nrow=8,ncol=8)

I<-matrix(c(1,0,0,0,0,0,0,0, 0,1,0,0,0,0,0,0, 0,0,1,0,0,0,0,0, 0,0,0,1,0,0,0,0, 0,0,0,0,1,0,0,0, 0,0,0,0,0,1,0,0, 0,0,0,0,0,0,1,0, 0,0,0,0,0,0,0,1), nrow=8,ncol=8)