function(filename)
{
	cat("Input alpha and beta parameters for the beta density, on separate lines, followed by return key :\n"
		)
	parms <- scan(what = numeric(), sep = "\n", strip.white = T)
	X11()
	cat("Input number of successes and number of failures, on separate lines, followed by return key :\n"
		)
	results <- scan(what = numeric(), sep = "\n", strip.white = T)
	#	postscript(filename, horizontal = F)
	x <- seq(0.005, 0.995, by = 0.01)
	toplim1 <- max(dbeta(x, parms[1] + results[1], parms[2] + results[
		2]))
	toplim2 <- max(dbeta(x, parms[1], parms[2]))
	toplim3 <- max(dbeta(x, results[1] + 1, results[2] + 1))
	toplim <- max(c(toplim1, toplim2, toplim3))
	print(c(toplim1, toplim2, toplim3, toplim))
	plot(x, dbeta(x, parms[1] + results[1], parms[2] + results[2]), type = 
		"l", lty = 3, ylim = c(0, 1.1 * toplim), xlab = "Value", ylab
		 = "Density", main = paste("Beta density,\nalpha = ", parms[
		1], ", beta = ", parms[2], " \n  Successes", results[1], 
		"Failures ", results[2]))
	lines(x, dbeta(x, results[1] + 1, results[2] + 1), type = "l", lty = 2)
	lines(x, dbeta(x, parms[1], parms[2]), type = "l", lty = 1)
	legend(0.6, 6, c("Prior", "Likelihood", "Posterior"), lty = c(1, 2,
		3))
	cat("To remove the graphics window, enter \n dev.off() \n at the prompt. \n"
		)
}