function dydt = fitzhugh(t,y,e,g,b,a,d,c,f)
dydt = zeros(2,1); 
dydt(1) = y(1) - (y(1)^3)/3 - y(2);
dydt(2) = y(1) + B - G*y(2);