function [time,y,nt]=euler(time,y,nt,dt);
global c;
 
y(nt+1)=y(nt)+dt*(y(nt)-y(nt)^2-c*y(nt)^3);
 
nt=nt+1;
