
 eps3 = eps1+ mu1*W./(V+mu2); 
 
 Fw = eps3.*(-W - gs*(V-eps2).*(V-alf-1)); 
 Fv =Istim - ga*((1-Isc).*(V-eps2).*(V-alf).*(V-1) +Isc.*(V-1)) - V.*W ;
