#

#
#       parameter values defining model
#
par alpha=0.5,gamma=1,delta=0.5,epsilon=0.01
#
#       definice funkcí
#
x'=x-epsilon*x^2-x*y/(1+alpha*x)
y'=-gamma*y-delta*y^2+x*y/(1+alpha*x)

#     
#
x(0)=0
y(0)=0
#
#
# Numerické parametry:
#
#
@ total=30,xlo=0,xhi=101,ylo=0,yhi=30,dt=0.01,maxstor=200000,xp=x,yp=y
#
done