clear all
close all
clc

disp('Ecuación: ')
 
disp('dp/dt = p*(m-(n*p))')
 
disp('___________________________________________________________________')
m=input('Intrinsic Growth Rate, m(1): ');
p=input('Population Size, p(10): ');
n=input('Density-dependent Crowding Effect (Intraspecific Competition), n(0.005): ');
t=input('Time, t(10): ');

f=@(t,p) p*(m-(n*p)); 
 
tspan=[0 t];
p0=p;
[t,p]=ode45(f,tspan,p0);

hold on
plot(t,p,'sr','MarkerEdgeColor','k',...
             'MarkerFaceColor','r')%,'LineWidth',2)
grid on
xlabel('t');
ylabel('p');
title('Crecimiento Logístico')
%__________________________________________________________________________