7. Presa-depredador

printf("Presa-depredador\n")
n=20;
m=10;
k=10;
t=50;
tp=(0:t-1);
Nt=zeros(t,1);
Nt2=zeros(t,1);
rd=3;
rp=2;
for j=1:t
  Nt(j)=n;
  Nt2(j)=m;
  Nt;
tablero=zeros(64,1);
  for i=1:n
    x=randi(64);
    tablero(x)=tablero(x)+1;
endfor
tablero2=zeros(64,1);
for i=1:m
    x1=randi(64);
    tablero2(x1)=tablero2(x1)+1;
  endfor
s=tablero;
s2=tablero2;
r=zeros(64,1);
r1=zeros(64,1);
r2=zeros(64,1);
r3=zeros(64,1);
so=zeros(64,1);
for i=1:64
  if (s2(i)==0)
    choice=1;
  endif
  if (s2(i)>0)
    choice=2;
  endif
switch choice
  case 1
    if (s(i)<k)
      r3(i)=s(i)*rp;
    endif
    if (s(i)==k)
      so(i)=s(i);
    endif
  case 2
    if (s(i)>s2(i))
      r(i)=s2(i)*rd;
    endif
    if (s(i)==s2(i))
      r1(i)=s2(i);
    endif
    if (s(i)<s2(i))
      r2(i)=s2(i)*0;
    endif
endswitch
end
rs=sum(r);
rs1=sum(r1);
rs2=sum(r2);
rs3=sum(r3);
so1=sum(so);
pasanP=rs3+so1;
pasanD=rs+rs1+rs2;
if (pasanD==0)
  pasanD=m;
endif
if (pasanP==0)
  pasanP=n;
endif
n=pasanP;
m=pasanD;
endfor
plot(tp,Nt,tp,Nt2)
xlabel ("t");
ylabel ("Nt");
title ("Interaccion presa-depredador");
legend ("Presa","Depredador");
t_vs_Nt_Nt2=[tp',Nt,Nt2]