现在有一段程序如下：
r<-6371.004
d=0
for(i in 1:2){
x2=b[i,2]
y2=b[i,3]
{for(j in 1:10){
x1=a[j,6]
y1=a[j,5]
z=x1-x2
d[i+j-1]=r*acos(cos(pi*y1/180)*cos(pi*y2/180)*cos(pi*z/180)+sin(pi*y1/180)*sin(pi*y2/180))
}
}
}
d
我想实现的是i循环两次，j循环10次，然后一共得到20个d的数据，但是d[i+j-1]最后出来的是11个，现知道此处有错误，请问怎样修改？此外，还想实现一个功能：想用if语句判断当d值小于某一数时（例如25），输出所有满足的b中第一列对应的数，这一步应该怎样实现，谢谢！