# hostpar.ode
# host-parasitoid model
hp=f*h/(1+a*p/k)^k
h'=hp
p'=h-hp/f
par f=1,a=1,k=2
@ meth=discrete,total=100
done