function logistic(a,b,C0,n)
C=zeros(1,n+1);
C(1)=C0;
if C(1)>a/b
   disp(['C0 no puede ser mayor de:  ' num2str(a/b)])
   return
end
for k=1:n
   C(k+1)=max((a-b*C(k))*C(k),0);
end
figure(1)
clf
plot([0:n],C)
hold on
plot([0:n],C,'x')
title(['Modelo de Crecimiento Logistico: a=' num2str(a) ',   b=' num2str(b)])
xlabel('Generaciones');ylabel('Poblacion')