clear all
close all

delta = 0.1;
dt    = 0.01;

x(1)=0.1;
y(1)=0;

i=2;

for t=dt:dt:15
    x(i) = x(i-1) + dt * y(i-1);
    y(i) = y(i-1) - dt * ( x(i-1) - x(i-1)*x(i-1)*x(i-1) - delta * y(i-1) );
    i=i+1;
end

t=0:dt:15;
plot (x,y,'g') % Grafico del diagrama de fases
xlabel('x')
ylabel('y')
hold off

figure
plot(t,x,t,y) % Grafico de la evolucion temporal


% %%% Grafico de las nuclinas %%%
% vv=0:.01:1;
% u1=A./(A+vv.*vv);
% u2=(B*vv-F)./(vv.*vv);
% 
% hold on
% plot(vv,u1,'r',vv,u2,'b')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%