function [time,y,nt]=euler(time,y,nt,dt);
global c ad bd delay;
 
if nt-delay > 0
  y(nt+1)=y(nt)+dt*(-y(nt)-y(nt)^2-c*y(nt)^3+ad*y(nt-delay)-bd*y(nt-delay)^2);
else
  y(nt+1)=y(nt)+dt*(-y(nt)-y(nt)^2-c*y(nt)^3+ad*y(nt)-bd*y(nt)^2);
%  y(nt+1)=y(nt)+dt*(-y(nt)-y(nt)^2-c*y(nt)^3);
end
nt=nt+1;

