n'=n*exp(r*(1-n/k)-a*p)
p'=n*(1-exp(-a*p))
par k=.5,r=.5,a=1
init n=.5,p=0
@ meth=disc,total=100
done
# r=1.6,k=5,n=.3766,p=.198  -- period 7 point