set.seed(1)
n <- 25
z <- cbind(runif(n,min=1,max=10),runif(n,min=1,max=10))
d <- compute_distance_matrix(z)
compute_gain_transp(sample(n),c(4,23),d) # -6.661
compute_gain_transp(sample(n),c(17,3),d) # 4.698
Run the code above in your browser using DataLab