p <- c(d=0.1,h=0.5,m=5,delta=1)

# Make Figures
size<-5#inch

pdf("deathA.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")
plot(1,1,type="n",xlim=c(0,1.2),ylim=c(0,1.2),xlab="N",ylab="Hill")
with(as.list(p),curve(d+delta*N^m/(h^m+N^m),from=0,to=1.2,xname="N",ylab="Hill",col="red",lwd=2,add=TRUE))
lines(c(0,1.2),c(0.6,0.6))
dev.off()

pdf("deathB.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")
plot(1,1,type="n",xlim=c(0,1.2),ylim=c(0,8),xlab="N",ylab="Hill")
with(as.list(p),curve(d*(1+(N/h)^m),from=0,to=1.2,xname="N",ylab="Power",col="red",lwd=2,add=TRUE))
lines(c(0,1.2),c(2,2))
dev.off()