# 
# Proie prédateur delay

f(x)=x*(1-x/K)
x'=r*f(x)-m*x*y
y'=-s*y+Z*exp(-s*tau)*m*delay(x,tau)*delay(y,tau)

aux ytau=delay(y,tau)

x(0)=0.1
y(0)=0.1
init x=0.1,y=0.1

par tau=0,k=1,m=1,r=1,Z=0.6,s=0.02

# set maxdelay
@ delay=180
# done
d

# vvvvvvvv ORIGINAL vvvvvvvv
# Proie prédateur delay

f(x)=x*(1-x/K)
x'=r*f(x)-m*x*y
y'=-s*y+Z*exp(-s*tau)*m*delay(x,tau)*delay(y,tau)

x(0)=0.1
y(0)=0.1
init x=0.1,y=0.1

par tau=0,k=1,m=1,r=1,Z=0.6,s=0.02

# set maxdelay
@ delay=180
# done
d