rm(list=ls())

pdf("cdf1.pdf")
y = seq(-5,5,length=250)
x = c(seq(-5,0,length=100),rep(0,50),seq(0,5,length=100))
F = pnorm(x)
plot(y,F,type="l",xlab="",ylab="",col="purple",lwd=4,cex.axis=1.3)
abline(h=.5,col="red",lwd=7,lty=3) 
abline(v=-1.024,col="red",lwd=5)
# lines(y,F,lwd=4,col="purple")
# abline(h=.5,col="red",lwd=7,lty=3) 
title(ylab=title(ylab=expression(F(x)),line=2.4,cex=1.8,cex.lab=1.7)) 
title(xlab=expression(F(1/2)^{-1}),line=3.5,cex.lab=1.7)
axis(side=2,at=.5,"0.5",line=-.8,las=2,col.axis="red",cex.axis=1.4,tck=-.01)
axis(side=1,at=-1.02,"-1.02",line=-.7,las=1,col.axis="red",tck=-.01,cex.axis=1.4)
dev.off()

pdf("cdf2.pdf")
y = seq(-5,5,length=200)
x1 = seq(-5,-.15,length=100)
x2 = seq(.15,5,length=100)
plot(y[1:100],pnorm(x1),ylim=c(0,1),type="l",xlim=c(-5,5),
               lwd=4,xlab="",ylab="",col="purple",cex.axis=1.3,xaxt="n")
lines(y[101:200],pnorm(x2),lwd=4,col="purple")
abline(h=.5,col="red",lwd=3)
abline(v=0,col="red",lwd=7,lty=3)   
points(0,.56,cex=1.5,pch=19,col="purple")
points(0,.44,pch=19,cex=1.5,col="purple")  
axis(side=2,at=.5,"0.5",line=-.8,las=2,col.axis="red",cex.axis=1.3,tck=-.01)
where = c(-4,-2,2,4)
axis(side=1,at=where,line=0,col.axis="black",cex.axis=1.4,tck=-.01)
axis(side=1,at=0,"0",line=0,col.axis="red",cex.axis=1.4,tck=-.01)
title(ylab=expression(F(x)),line=1.9,cex.lab=1.5)
title(xlab=expression(F(1/2)^{-1}),line=3.5,cex.lab=1.5)

dev.off()