showRadiation <- function(x,lambda,...) {
mx <- max(x)
mn <- min(x)
    result <- histogram(~x,type="density", xlim=c((mn-2.5),(mx+2.5)),
	xlab='Count Rate',
        panel=function(x,y,...){
	      panel.histogram(x,alpha=0.8,breaks=seq((mn-2.5),(mx+2.5),by=1,...))
	      panel.points((mn-2):(mx+2),dpois((mn-2):(mx+2),lambda),pch=19,alpha=0.8)
	}
    )
  return(result)
}

showRadiation0 <- function(x,lambda,...) {
mx <- max(x)
mn <- 0
    result <- histogram(~x,type="density", xlim=c((mn-2.5),(mx+2.5)),
	xlab='Count Rate',
        panel=function(x,y,...){
	      panel.histogram(x,alpha=0.8,breaks=seq((mn-2.5),(mx+2.5),by=1,...))
	      panel.points((mn-2):(mx+2),dpois((mn-2):(mx+2),lambda),pch=19,alpha=0.8)
	}
    )
  return(result)
}