counting goes through years)

for(k in 0:(limsup-1)){prod\_2\_init = prod\_2\_init\*((sumalphacount-1)-k)^(-1)}

count\_Km <- matrix(0, 2^K, 2) # matrix counting the transitions from m
