function bif_hnnmap
bif_hnn(@weifen4,'bif_odef_sign.mat');

function bif_hnn(odef,filename)
op=odeset('RelTol',1e-6,'AbsTol',1e-6);
ope=odeset('RelTol',1e-6,'AbsTol',1e-6,'Events',@events);
x0=[0 0 1];
tic;
p=(0.46:0.01:0.54);
xp=[];
for i=1:length(p);
a=p(i);
[t,x]=ode45(odef,[0,50],x0,[],a);
[t,x,te,xe,ie]=ode45(odef,[0,150],x(end,:),ope,a);
xp=[xp;xe(ie==1,[1,3]) ones(sum(ie==1),1)*p(i)];
save(filename,'xp');
disp(xp);
end
plot(xp(:,end),xp(:,1),'.','markersize',2);

%------------------------------
function [value,isterminal,direction]=events(t,x,a,w)
value=[x(2)];
isterminal=[0];
direction=[1];
%------------------------------

function dy=weifen4(t,x,a)

if x(1)< -a
    K=-x(1)-1;
elseif x(1)<a&&x(1)>-a 
    K=(1-a)*x(1)/a;
else
    K=-x(1)+1;
end

dy=[x(2);...
   x(3);...
   -0.6*x(3)-x(2)+K;];