%Example on two stages and two ages %ns=number of stages ns=4; m1 = 0.5; %%0.7; fertility age 1 s1 = 0.55;% survival age 1 to age 2 m2 = 1.1; %%0.75*1.02; fertility age 2 x1=0; y1=0; p1=0; F1 = zeros(ns,ns); D1 = F1; m11 = m1; %fertility for age 1 size small m12 = m1; %fertility for age 1 size large M1 = diag([m11 m12 0 0]); d11 = 1; D1(1,1) = d11; D1(2,1) = (1-d11); d12 = 0.8; D1(1,2) = d12; D1(2,2) = (1-d12); F1 = D1*M1; P1 = zeros(ns,ns); G1 = P1; %p1 = 0.5; %neutral growth, half to each %b1g = 1.5; %ranges from 2 (relative size stays forever) % to 0 (size ranking switches) p11 = 0.5; s11 = s1; %fert for age 1 size small s12 = s1; %fert for age 1 size large S1 = diag([s11 s12 0 0]); G1(3,1) = p11; G1(4,1) = (1-p11); G1(3,2) = (1-p11); G1(4,2) = p11; P1 = G1*S1; F2 = zeros(ns,ns); D2 = F2; %m2 = mult*0.75*1.02; %x2 = 0.25; %sel on fert d21 = 0.5; %no heritability m21 = m2; %fert for age 1 size small m22= m2; %fert for age 1 size large M2 = diag([0 0 m21 m22]); D2(1,3) = d21; D2(2,3) = (1-d21); D2(1,4) = (1-d21); D2(2,4) = d21; F2 = D2*M2; outvals = smallmat(ns, F1, F2, P1);%this is (2*ns+8) long save small_neutral %%outvals = [uu' uu2' vv' R0 Tc r0 a1 a1+a2 a1+a2+a3 xi1 xi phi]