Traté de hacer una repetición de la secuencia de este modelo de la siguiente manera.
En primer lugar, podemos hacer las siguientes definiciones:
$N_{d,m}^{(k)}$ es el número de hombres la CONDENACIÓN de los transportistas en el $k$-ª generación.
$N_{d,f}^{(k)}$ es el número de mujeres DOOM portadores en el $k$-ª generación.
$N_{n,m}^{(k)}$ es el número de hombres DOOM no portadores de la $k$-ª generación.
$N_{n,f}^{(k)}$ es el número de mujeres que DOOM no portadores de la $k$-ª generación.
$T^{(k)}$ es el número de parejas en el $k$-ª generación.
$T_{d,d}^{(k)}$ es el número de parejas cada uno consta de 2 DOOM portadores en el $k$-ª generación.
$T_{d,n}^{(k)}$ es el número de parejas cada uno consta de 1 macho DOOM portador y 1 hembra DOOM no es portador de la $k$-ª generación.
$T_{n,d}^{(k)}$ es el número de parejas cada uno consta de 1 macho no-DOOM portador y 1 hembra DOOM transportista en el $k$-ª generación.
$T_{n,n}^{(k)}$ es el número de parejas cada uno consta de 2 DOOM no portadores de la $k$-ª generación.
Entonces tenemos:
$T^{(k)}=min((N_{d,m}^{(k)}+N_{n,m}^{(k)}) , (N_{d,f}^{(k)}+N_{n,f}^{(k)}))$
$T_{d,d}^{(k)}=T^{(k)}\frac{N_{d,m}^{(k)}}{N_{d,m}^{(k)}+N_{n,m}^{(k)}}\frac{N_{d,f}^{(k)}}{N_{d,f}^{(k)}+N_{n,f}^{(k)}}$
$T_{d,n}^{(k)}=T^{(k)}\frac{N_{d,m}^{(k)}}{N_{d,m}^{(k)}+N_{n,m}^{(k)}}\frac{N_{n,f}^{(k)}}{N_{d,f}^{(k)}+N_{n,f}^{(k)}}$
$T_{n,d}^{(k)}=T^{(k)}\frac{N_{n,m}^{(k)}}{N_{d,m}^{(k)}+N_{n,m}^{(k)}}\frac{N_{d,f}^{(k)}}{N_{d,f}^{(k)}+N_{n,f}^{(k)}}$
$T_{n,n}^{(k)}=T^{(k)}\frac{N_{n,m}^{(k)}}{N_{d,m}^{(k)}+N_{n,m}^{(k)}}\frac{N_{n,f}^{(k)}}{N_{d,f}^{(k)}+N_{n,f}^{(k)}}$
Como se supone en el modelo, una pareja con al menos 1 DOOM transportista tiene la probabilidad de $P_m$ a dar a luz al niño.
En el $k$-ésima generación, el número de parejas con al menos 1 DOOM transportista es $T_{d,d}^{(k)}+T_{d,n}^{(k)}+T_{n,d}^{(k)}$.
Por lo tanto, tenemos:
$N_{d,m}^{(k+1)}=P_m (T_{d,d}^{(k)}+T_{d,n}^{(k)}+T_{n,d}^{(k)}) N$
$N_{d,f}^{(k+1)}=(1-P_m) (T_{d,d}^{(k)}+T_{d,n}^{(k)}+T_{n,d}^{(k)}) N$
$N_{n,m}^{(k+1)}=N_{n,f}^{(k+1)}=0.5 T_{n,n}^{(k)} N$
A continuación, podemos obtener:
$\frac{(N_{d,m}^{(k+1)}+N_{d,f}^{(k+1)})}{(N_{n,m}^{(k+1)}+N_{n,f}^{(k+1)})}\frac{(N_{n,m}^{(k)}+N_{n,f}^{(k)})}{(N_{d,m}^{(k)}+N_{d,f}^{(k)})}=\frac{N_{d,m}^{(k)}N_{d,f}^{(k)}+N_{d,m}^{(k)}N_{n,f}^{(k)}+N_{n,m}^{(k)}N_{d,f}^{(k)}}{N_{n,m}^{(k)}N_{n,f}^{(k)}}\frac{(N_{n,m}^{(k)}+N_{n,f}^{(k)})}{(N_{d,m}^{(k)}+N_{d,f}^{(k)})}>1$
lo que significa que la proporción de la PERDICIÓN de los transportistas siempre va en aumento, lo que indica que el destino de los portadores tomar la totalidad de la población con el tiempo. Así obtenemos
$\displaystyle \lim_{k \rightarrow\infty}N_{n,m}^{(k)}=0$ $\displaystyle \lim_{k \rightarrow\infty}N_{n,f}^{(k)}=0$
Además, podemos tener:
$\frac{T^{(k+1)}}{T^{(k)}}=\frac{min((N_{d,m}^{(k+1)}+N_{n,m}^{(k+1)}) , (N_{d,f}^{(k+1)}+N_{n,f}^{(k+1)}))}{T^{(k)}}= \frac{min((N_{d,m}^{(k+1)}+N_{n,m}^{(k+1)}) , (N_{d,f}^{(k+1)}+N_{n,f}^{(k+1)}))}{T^{(k)}}=\frac{min((\frac{P_{m}T^{(k)}(N_{d,m}^{(k)}N_{d,f}^{(k)}+N_{d,m}^{(k)}N_{n,f}^{(k)}+N_{n,m}^{(k)}N_{d,f}^{(k)})N}{(N_{d,m}^{(k)}+N_{n,m}^{(k)})(N_{d,f}^{(k)}+N_{n,f}^{(k)})}+N_{n,m}^{(k+1)}) , (\frac{(1-P_{m})T^{(k)}(N_{d,m}^{(k)}N_{d,f}^{(k)}+N_{d,m}^{(k)}N_{n,f}^{(k)}+N_{n,m}^{(k)}N_{d,f}^{(k)})N}{(N_{d,m}^{(k)}+N_{n,m}^{(k)})(N_{d,f}^{(k)}+N_{n,f}^{(k)})}+N_{n,f}^{(k+1)}))}{T^{(k)}}=min((\frac{P_{m}(N_{d,m}^{(k)}N_{d,f}^{(k)}+N_{d,m}^{(k)}N_{n,f}^{(k)}+N_{n,m}^{(k)}N_{d,f}^{(k)})N}{(N_{d,m}^{(k)}+N_{n,m}^{(k)})(N_{d,f}^{(k)}+N_{n,f}^{(k)})}+N_{n,m}^{(k+1)}) , (\frac{(1-P_{m})(N_{d,m}^{(k)}N_{d,f}^{(k)}+N_{d,m}^{(k)}N_{n,f}^{(k)}+N_{n,m}^{(k)}N_{d,f}^{(k)})N}{(N_{d,m}^{(k)}+N_{n,m}^{(k)})(N_{d,f}^{(k)}+N_{n,f}^{(k)})}+N_{n,f}^{(k+1)}))$
Como $k \rightarrow\infty$, $N_{n,m}^{(k)}$ y $N_{n,f}^{(k)}$ $\rightarrow 0$, tenemos
$\displaystyle \lim_{k \rightarrow \infty}\frac{T^{(k+1)}}{T^{(k)}} =\lim_{k \rightarrow \infty} min(\frac{P_{m}(N_{d,m}^{(k)}N_{d,f}^{(k)}+N_{d,m}^{(k)}N_{n,f}^{(k)}+N_{n,m}^{(k)}N_{d,f}^{(k)})N}{(N_{d,m}^{(k)}+N_{n,m}^{(k)})(N_{d,f}^{(k)}+N_{n,f}^{(k)})},\frac{(1-P_{m})(N_{d,m}^{(k)}N_{d,f}^{(k)}+N_{d,m}^{(k)}N_{n,f}^{(k)}+N_{n,m}^{(k)}N_{d,f}^{(k)})N}{(N_{d,m}^{(k)}+N_{n,m}^{(k)})(N_{d,f}^{(k)}+N_{n,f}^{(k)})}) = min(P_{m}N,(1-P_{m})N)$
Si N<2, $\frac{T^{(k+1)}}{T^{(k)}}<1$, lo que significa que el número de parejas que se disminuirá de manera consistente en el largo plazo, como resultado, la población se extinguieron sin importar el valor de $P_m$.
Cuando N>=2,
Si $(\frac{1}{N}-P_m)(1-\frac{1}{N}-P_m)<0$, $\displaystyle \lim_{k \rightarrow \infty}\frac{T^{(k+1)}}{T^{(k)}}>1$, y la población va a aumentar.
Si $(\frac{1}{N}-P_m)(1-\frac{1}{N}-P_m)>0$, $\displaystyle \lim_{k \rightarrow \infty}\frac{T^{(k+1)}}{T^{(k)}}<1$, y la población se extinguió.
Si $(\frac{1}{N}-P_m)(1-\frac{1}{N}-P_m)=0$, $\displaystyle \lim_{k \rightarrow \infty}\frac{T^{(k+1)}}{T^{(k)}}=1$, y la población alcanzará un tamaño estable con el tiempo.
Las condiciones iniciales ($N_{d,m}^{(1)}$, $N_{d,f}^{(1)}$, $N_{n,m}^{(1)}$, $N_{n,f}^{(1)}$) no va a afectar a las conclusiones sobre si y sólo si estos valores son positivos.