n'=n*exp(r*(1-n))*exp(-a*p)
p'=n*(1-exp(-a*p))
@ meth=discrete
@ total=400
init n=.2,p=.3
par a=0.1,r=1.5
done