%TMMC 2021 - Sex Ratios

clear all
close all
clc


% 0 -> tipo 0
% 1 -> tipo 1

p = 0.5 % Propoción inicial de los tipos.
t = 100 % Generaciones 
n = 100 % Tamaño inicial de la población.
N = 2 % Cantidad de hijos por hembra.
s1 = 0.6 % Proporción de machos del tipo 1
s0 = 0.45 % Proporción de machos del tipo 2 
s = p*s1 + (1-p)*s0 % Propoción de machos en la población 

C1 = zeros(1,t); % Vectores con la cantidad de cada tipo a través de las generaciones.
C0 = zeros(1,t);

cant1 = p*n % Cálculo de la cantidad de cada tipo
cant0 = (1 - p)*n

C1(1) = cant1;
C0(1) = cant0;

for i = 1:t
  
  Wss1 = (N^2)*(1 - s1 + s1*((1-s)/s)) % Fitness tipo 1.
  Wss0 = (N^2)*(1 - s0 + s0*((1-(s))/(s))) % Fitness tipo 2. 
  
  p = p*(Wss1)/(p*(Wss1) + (1-p)*(Wss0))  % Cálculo de la proporción de los tipos
  s = p*s1 + (1-p)*s0 % Cálculo de la proporción de machos en la población
  cant1 = p*n % Cálculo de las nuevas cantidades
  cant0 = (1-p)*n
  
  
  C1(i+1) = cant1
  C0(i+1) = cant0
  
endfor 

plot(1:t+1,C1,'b',C0,'r') 

cant1/n
cant0/n
s