编辑代码

Location<-array(1:22,dim=c(11,2))
Location
Length<-array(dim=c(11,11))

for (i in 1:11){
  for( j in 1:11){
    
     if(i >= j)
        next

Lng1<-Location[i,1]
Lat1<-Location[i,2]
Lng2<-Location[j,1]
Lat2<-Location[j,2]

Length[i,j]<-2*asin(sqrt(sin((Lat1-Lat2)/2)*sin((Lat1-Lat2)/2)+cos(Lat1)*cos(Lat2)*sin((Lng1-Lng2)/2)*sin((Lng1-Lng2)/2)))*6378.137

}
}
Length