Srange<-0:2000*100 Rhat<-alpha*Srange/(1+Srange/K) plot(Srange,Rhat,type='l',xlab="S ('000 t)",ylab="R (millions)",ylim=c(0,10000)) sr2<-srfun(selF2,sa,M,wa,pa) sr3<-srfun(selF3,sa,M,wa,pa) sr.crash<-srfun(Fcrash,sa,M,wa,pa) lines(Srange,Srange/sr.crash) text(150000,8500,"Beverton-Holt curve") text(150000,8000,expression(R=alpha*S/(1+S/K))) text(15000,6500,"Slope for\n Fcrash is") text(15000,6000,expression(alpha)) S0<-10000 R0<-alpha*S0/(1+S0/K) S1<-R0*sr2 R1<-alpha*S1/(1+S1/K) S2<-R1*sr2 R2<-alpha*S2/(1+S2/K) S3<-R2*sr2 R3<-alpha*S3/(1+S3/K) S4<-R3*sr2 R4<-alpha*S4/(1+S4/K) lines(c(S0,S0,S1,S1,S2,S2,S3,S3,S4),c(0,R0,R0,R1,R1,R2,R2,R3,R3)) lines(Srange,Srange/sr2)