6. Competencia interespecífica

printf("Competencia interespecífica \n")
# Especie A
n=20;
# Especie B
m=15;
# capacidad de carga especie A
k1=2;
# capacidad de carga especie B
k2=3;
# Tiempo de repeticiones
t=200;
tp=(0:t-1);
NtA=zeros(t,1);
NtB=zeros(t,1);
# Fecundidad de especie A
fd=3;
# Fecundidad de especie B
fp=4; 
for j = 1:t
    NtA(j) = n;
    NtB(j) = m;
    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

    sA = tablero;
    sB = tablero2;
    rA1 = zeros(64,1);
    rA2 = zeros(64,1);
    rA3 = zeros(64,1);
    rB1 = zeros(64,1);
    rB2 = zeros(64,1);
    rB3 = zeros(64,1);
 for i = 1:64
        if (sB(i) == 0)
            choice = 1;
        endif
        if (sB(i) > 0)
            choice = 2;
        endif
 switch choice
            case 1
                if (sA(i) < k1)
                    rA1(i) = sA(i) * fp;
                endif
                if (sA(i) == k1)
                    rA2(i) = sA(i);
                endif

            case 2
                if (sA(i) > sB(i))
                    rA3(i) = sA(i);
                elseif (sA(i) == sB(i))
                    rA4(i) = 0;
                elseif (sA(i) < sB(i))
                    if (sA(i) == 0)
                        if (sB(i) == k2)
                            rB1(i) = sB(i);
                        elseif (sB(i) < k2)
                            rB2(i) = sB(i) * fp;
                        endif
                    elseif (sA(i) > 0)
                        rB3(i) = sB(i);
                    endif
                endif
        endswitch
    endfor
rA1p = sum(rA1);
    rA2p = sum(rA2);
    rA3p = sum(rA3);
    rB1p = sum(rB1);
    rB2p = sum(rB2);
    rB3p = sum(rB3);

    pasanA = rA1p + rA2p + rA3p;
    pasanB = rB1p + rB2p + rB3p;

    if (pasanA == 0)
        pasanA = m;
    endif
    if (pasanB == 0)
   pasanB = n;
    endif
    n = pasanA;
    m = pasanB;
endfor
plot(tp,NtA,tp,NtB)
xlabel("t");
ylabel("Nt");
title("Competencia interespecífica");
t_vs_Nt_Nt2=[tp',NtA,NtB];