SpDist1 <- function(data) {

  dist <- spDists(meuse)
  n <- length(data)
  for (i in 1:(n - 1)) {
    for (j in (i+1):n) {
      d <- dist[i,j]
      dd <- (c(i,j,d))
    }

}
  return(dd)
}