刚刚开始学习R 不久，在对大矩阵A分块运行的时候遇到这么一个问题，代码运行老是报错，不知道错误出在哪里,以下是程序：
i<-1 #第一行处理
m<-matrix(0,nrow=30,ncol=30)
I<-diag(30)
for (j in 1:31) {
  if (j==i) Int1_<-matrix(0,nrow=30,ncol=1)
  else {
    Aij<-A[(1+30*(i-1)):(30*i),(1+30*(j-1)):(30*j)]
    Ajjn<-solve(I-A[(1+30*(j-1)):(30*j),(1+30*(j-1)):(30*j)])
    Int1_<-Aij%*%Ajjn%*%Y_pro[(1+30*(j-1)):(30*j)]
  }
  Int1_<-as.matrix(Int1_)
  m[,j]<-Int1_
}


运行后报错：
Error in `[<-`(`*tmp*`, , j, value = Int1_) : subscript out of bounds

运行之后有结果，同时也报错了，不知道啥原因，
请问哪位前辈知道哪里出错了嘛？非常感谢