png("samplingdist.png")
hist(spropDist,breaks=c(0,.24,seq(.33,.75,.06),1),xaxt="n",xlab="",ylab="",main="Sampling Distribution")
axis(1,at=c(0,.24,seq(.33,.75,.06),1),las=2)
curve(dnorm(x,.54,sqrt(.54*.46/80)),col="red",add=T)
dev.off()

show(icosahedron(color='green', opacity=.5, mesh=3), spin=1)