n'=r*n*(1-n-a*p)
p'=p*(d+c*n)
par r=1.5,a=.1,d=.8,c=.8
init n=.6,p=.2
@ meth=discrete,total=300
done