#####   R-code for making Figure 2.9

	dat  <- read.table("soil.dat",header=F)
	x    <- dat$V2

	postscript("fig29.ps",width=5,height=5,horizontal=F)

	hist(x, breaks=c(0,0.25,0.5,0.75,1,3.75), freq=F, right=F,
	     xlab=expression(SiO[2]),mgp=c(2,1,0),cex=0.8,main="")

	fit <- density(x)

	lines(fit,lty=1,cex=0.8)

	graphics.off()