clear all
clc

disp('Ecuación Logística: ')
 
disp('dN/dt = r*N*((K-N)/K)')
 
disp('___________________________________________________________________')

disp('                                                                   ')

r=input('Tasa de Crecimiento Poblacional per capita, r (1): ');
N0=input('Población Inicial, N0 (100): ');
K=input('Máximo Valor de N (Capacidad de Carga del Sistema), K (10000): ');
t=input('Tiempo, t (30): ');

disp('                                                                  ')

% Código Básico de Colores en Matlab 
disp('Red [1 0 0],      Green [0 1 0],  Blue [0 0 1],   Yellow [1 1 0]')
disp('Magenta [1 0 1],  Cyan [0 1 1],   Black [0 0 0],  White [1 1 1]')

disp('                                                                  ')

C1=input('1er. Valor, Elija 0 o 1: ');
C2=input('2do. Valor, Elija 0 o 1: ');
C3=input('3er. Valor, Elija 0 o 1: ');

f=@(t,N) r*N*((K-N)/K); 

tspan=[0 t];
[t,N]=ode45(f,tspan,N0);

% N

% Al cambiar los valores de R, G, B en cada nuevo estudio es posible
% identificar las nuevas curvas en el gráfico.
R=C1;
G=C2;
B=C3;
 
hold on
plot(t,N,'o','MarkerFaceColor',[R,G,B],...
             'MarkerEdgeColor','k',...
             'MarkerSize',3)
grid on
xlabel('t');
ylabel('N');
title('Crecimiento Logístico')
%__________________________________________________________________________