pp.tripleden <-
function (x, y, z, lwd=3, colors=c("gold", "steelblue", "black"), ...) 
{
	# placed in the public domain 2012 by Burns Statistics

	dx <- density(x)
	dy <- density(y)
	dz <- density(z)
	plot(dx, xlim=range(dx$x, dy$x, dz$x), ylim=range(dx$y, dy$y, dz$y),
		col=colors[1], yaxt="n", lwd=lwd, ...)
	
	lines(dy, col=colors[2], lwd=lwd)
	lines(dz, col=colors[3], lwd=lwd)
}