# predator prey with allee effect for prey
x'=x*(a1+a2*x-a3*x^2-b*y)
y'=y*(-d+g*x)
par d=1,g=1,a1=.1,a2=1,a3=1,b=1
@ xp=x,yp=y,xlo=-.05,xhi=1.2,ylo=-.05,yhi=1.2
@ nmesh=100
@ total=200,nout=10
done