# smalltumor_immunemodel.ode

dx/dt=r*x*(1-x) -dx*x - g*x*w*exp(-c*x)
dw/dt=s + a*x*w*exp(-c*x) - dw*w

p s=0.05,a=1,dw=0.2
p r=1.0,dx=0.05,g=1,c=2.0

init x=1,w=1

@ meth=qualrk,total=1000

@ maxstor=100000
@ xp=x,yp=w,xlo=0,ylo=0,xhi=2,yhi=2

# AUTO bifurcation settings
@ autoxmin=0,autoxmax=0.25,autoymin=-0.1,autoymax=2
@ dsmax=0.1,dsmin=.001,parmin=0,parmax=0.25

done


d