powerCurve <- function(mu,sigma,n,alpha0) {
  c <- qnorm(1-(alpha0/2),0,sigma/sqrt(n))
  x <- seq(mu-2,mu+2,0.01)
  power <- (1-pnorm(mu+c,x,sigma/sqrt(n))) + (pnorm(mu-c,x,sigma/sqrt(n)))
  temp <- paste("N =",n," sigma =",sigma," alpha0 =",alpha0);
  plot(x,power,type="l",ylab="power",xlab=temp)
}

par(mfrow=c(4,1))
#powerCurve(5,0.1,10,0.05)
#powerCurve(5,1,10,0.05)
#powerCurve(5,5,10,0.05)
#powerCurve(5,10,10,0.05)

powerCurve(5,1,10,0.001)
powerCurve(5,1,10,0.01)
powerCurve(5,1,10,0.05)
powerCurve(5,1,10,0.5)

#powerCurve(5,1,10,0.05)
#powerCurve(5,1,30,0.05)
#powerCurve(5,1,100,0.05)
#powerCurve(5,1,1000,0.05)