function xp = fhn()

epsilon = 0.08; % 0.08;
b = 0.7;
c = 0.8;
I = 0.32;

xp = [-1/3, 0, 1-1/c, I-b/c];