library(tigerstats)

(m<-matrix(c(1,2,3,3), ncol=2, byrow=T))
layout(mat=m)
qnormGC(0.95,region="below",m=40,s=10, graph=T)
qnormGC(0.95,region="below",m=60,s=20, graph=T)
qnormGC(0.95,region="below",m=0,s=1, graph=T)

qnormGC(0.975,region="below",m=0,s=1, graph=T)
qnormGC(0.995,region="below",m=0,s=1, graph=T)

set.seed(1234)
n<-100000
z<-rnorm(n,0,1)
(mean(z))
(sd(z))
par(mfrow=c(2,1))
hist(z, freq=F,xlim=c(-3,3),breaks=1000)
par(new=T)
plot(density(z), axes=F, main="", xlim=c(-4,4), lwd=2, col="blue")
qnormGC(0.975,region="below",m=0,s=1, graph=T)