clear all
close all

    A=0.02; %
    B=0.07; % Parametros
    F=0; %
    ci=[1 .1]; % Condiciones iniciales


alfa = 1;

dt = 1;

x(1)=ci(1);
y(1)=ci(2);

i=2;


for t=dt:dt:1000
    x(i) = x(i-1) + dt * ( -x(i-1)*y(i-1)*y(i-1) + A*(1-x(i-1)) );
    y(i) = y(i-1) - dt * (  x(i-1)*y(i-1)*y(i-1) - B*y(i-1) + F );
    i=i+1;
end

t=0:dt:1000;
%%% 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')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
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