pdf("cdfexample.pdf")

n = 20
Y = rnorm(n)


par(mfrow=c(1,1))

grid = seq(-3,3,length=100)
plot(grid,pnorm(grid),type="l",lwd=3,xlab="",ylab="F(y)")

Y = sort(Y)
Y = c(-3,Y,3)
F = c(0,(1:n)/n,1)
lines(Y,F,type="s",lwd=3,col="red")
dev.off()


cat("Median = ",median(Y),"\n")

a = round(n*(.5 - 2*sqrt(1/(4*n))))
b = round(n*(.5 + 2*sqrt(1/(4*n))))
print(a)
print(b)
print(c(Y[a],Y[b]))

dev.off()