4.1 La siguiente rutina tiene asignados los valores de la siguientes forma: tiempo (t=20), capacidad de carga (k=2), población inicial (n=5) y tasa de crecimiento (r=7).
printf("Logistico oscilatorio\n");
t=20;
k=2;
n=5;
r=7;
x=zeros(64,1);
tp=(0:t-1);
Nt=zeros(t,1);
for j=1:t
Nt(j)=n;
Nt;
tablero=zeros(64,1);
for i=1:n
x=randi(64);
tablero(x)=tablero(x)+1;
endfor
s=tablero;
rr=zeros(64,1);
for i=1:64
if (s(i)<k)
rr(i)=s(i)*r;
end
end
rr;
r1=sum(rr);
so=zeros(64,1);
for i=1:64
if (s(i)==k)
so(i)=s(i)*1;
end
end
so;
so1=sum(so);
pasan=r1+so1;
n=pasan;
endfor
plot(tp,Nt)
xlabel ("t");
ylabel ("Nt");
title ("Crecimiento logistico Nt vs pt");
t_vs_Nt=[tp',Nt]4.2 En la siguiente rutina el usuario debe asignar los valores de: tiempo (t), capacidad de carga (k), población inicial (No) y tasa de crecimiento (r). Se recomienda usar tasas de crecimiento de 7 o mayores, y una capacidad de carga igual a 2.
printf("Logistico oscilatorio\n");
t=input('Ingrese el numero de repeticiones, t: ');
k=input('Ingrese la capacidad de carga, k: ');
n=input('Ingrese la poblacion inicial, No: ');
r=input('Ingrese la tas de crecimiento, r: ');
x=zeros(64,1);
tp=(0:t-1);
Nt=zeros(t,1);
for j=1:t
Nt(j)=n;
Nt;
tablero=zeros(64,1);
for i=1:n
x=randi(64);
tablero(x)=tablero(x)+1;
endfor
s=tablero;
rr=zeros(64,1);
for i=1:64
if (s(i)<k)
rr(i)=s(i)*r;
end
end
rr;
r1=sum(rr);
so=zeros(64,1);
for i=1:64
if (s(i)==k)
so(i)=s(i)*1;
end
end
so;
so1=sum(so);
pasan=r1+so1;
n=pasan;
endfor
plot(tp,Nt)
xlabel ("t");
ylabel ("Nt");
title ("Crecimiento logistico Nt vs pt");
t_vs_Nt=[tp',Nt]