Fw = epsilon*(V - gam*W);
Fv = -V.*(V-1).*(V-alf) - W ;
