f <- function(x,h,n) { return (x^n/(h^n + x^n))}
g <- function(x,h,n) { return (1-exp(-log(2)*(x/h)^n))}

curve(f(x,0.25,1),from=0,to=1,ylim=c(0,1),lwd=2,col="red",ylab="f(x)")
curve(g(x,0.25,1),lwd=2,col="blue",add=T)

curve(1-f(x,0.25,1),from=0,to=1,ylim=c(0,1),lwd=2,col="red",ylab="f(x)")
curve(1-g(x,0.25,1),lwd=2,col="blue",add=T)

curve(f(x,0.25,2),from=0,to=1,ylim=c(0,1),lwd=2,col="red",ylab="g(x)")
curve(g(x,0.25,2),lwd=2,col="blue",add=T)

pdf("funca.pdf",width=size,height=size)
par(mar=c(0.1,0.1,0.1,0.1),xaxt="n",yaxt="n",ann=F,xaxs="i",yaxs="i")
curve(f(x,0.25,1),from=0,to=1,ylim=c(0,1),lwd=2,col="red",ylab="f(x)")
curve(g(x,0.25,1),lwd=2,col="blue",add=T)
dev.off()

pdf("funcb.pdf",width=size,height=size)
par(mar=c(0.1,0.1,0.1,0.1),xaxt="n",yaxt="n",ann=F,xaxs="i",yaxs="i")
curve(f(x,0.25,2),from=0,to=1,ylim=c(0,1),lwd=2,col="red",ylab="g(x)")
curve(g(x,0.25,2),lwd=2,col="blue",add=T)
dev.off()