function out1=chemovctfld(X,parameters)
%
%
%

a1=parameters(1);
a2=parameters(2);
a3=parameters(3);
m1=parameters(4);
m2=parameters(5);
m3=parameters(6);
eps=parameters(7);
del=parameters(8);

x=X(1).*(X(1)>0);
y=X(2).*(X(2)>0);
z=X(3).*(X(3)>0);
f1=m1*(1-x-y-z)./(a1+(1-x-y-z));
f2=m2./(a2+x);
f3=m3./(a3+y);

  out1=[(x+1e-6*(x<=0)).*(f1-1-y.*f2)/eps;
    (y+1e-6*(y<=0)).*(x.*f2-1-z.*f3);
    del*(z+1e-6*(z<=0)).*(y.*f3-1)];