# plot histogram of D
t<-read.table("Dcoeffs.dat")
D<-t[[1]]*1.e5

postscript(file="Dhist.eps",horizontal=FALSE,width=5,height=5)
hist(D,xlim=c(0,5),xaxp=c(0,5,5),
     main="",xlab=expression('D (10'^-5*'cm'^2*'/s)'))