d1 <- read.table("distance116.dat")
d2 <- read.table("distance129.dat")


pdf("dist_hist.pdf", w=10, h=5)
par(mfrow=c(1,2))
hist(d1[,2])
hist(d2[,2])
dev.off()