为什么diag（a[b[,2],b[,3]]） <-  b[,1] 运算出来的结果并不是将b[,1]赋值给a[b[,2],b[,3]]的对角线上呢？
具体情况见下图