function out1 = threevctrfld(x,param)
%
% 

a=param(1);
eps=param(2);
gamma=param(3);
c=param(4);
scl=1;

v=x(1);u=x(2);w=x(3);

f=v.*(v-a).*(1-v);
out1 = [u; 
        c*u/scl-f+w;
        scl*eps*(v-gamma*w)/c];
end
