# EE model
S'=r*S*(1-(S/K))-((lambda*E)/(1+E))*S-((alpha1*S*Y)/(beta+S+(m*E)))
E'=((lambda*E)/(1+E))*S-((alpha2*E*Y)/(beta+S+(m*E)))-mi1*E
Y'=((theta1*S+theta2*E)/(beta+S+(m*E)))*Y-mi2*Y
par r=1,K=500,lambda=1,alpha1=1,alpha2=1,beta=50,m=.75,theta1=.75,theta2=.5,mi1=.3,mi2=.1
init S=10,E=10,Y=10
done