% population model

% J_n -> J_{n+1} proportion of juveniles who 
% neither die nor reach age 15 during the (n+1)st year after 1990
M(1,1) = (1-0.0008)*14/15

% J_n -> A_{n+1} proportion of juveniles who 
% survive and reach age 15 during the (n+1)st year after 1990
M(2,1) = (1-0.0008)/15

% A_n -> J_{n+1} proportionate birth rate of adults
M(1,2) = 0.0172

% A_n -> A_{n+1} proportion of adults who do not die during the year
M(2,2) = 1-0.0136

init=[10.92;46.49];

x=init;
for k=1:20
    x=M*x;
end

M^40*init

[vecs, vals]=eigs(M)

AtoJratio=vecs(2,1)/vecs(1,1)

