distance <- function(p1,p2){
    xdiff = p1["x"] - p2["x"]
    ydiff = p1["y"] - p2["y"]
    diff = sqrt(xdiff^2 + ydiff^2)
    as.numeric(diff)
}