



clear t
clear y
tstep=.01
t=0:tstep:20;
y(:,1)=[-1;1];p=[.2,.2,3];
for lp=2:length(t)
    dxdt= fhnfunode(t(lp),y(:,lp-1),p);
    y(:,lp)=dxdt*tstep+y(:,lp-1);
end
hold on
plot(t,y,'k')

%%

[t,yode45]=ode45(@fhnfunode,t,y(:,1),[],p);

hold on
plot(t,yode45,'b')