%LamdaS Approximation for Ardisia escallonioides (eg., Tulja et al., 2003) in Markovian environments
%load demographic projection matrices and transition prob. matrices
%load COSMBASE
%load hurset


K = 7;  %7 habitat states

B = [m1,m2,m3,m4,m5,m6,m7]; %m1, m2, .., m7 are the projection matrices
S = 8; %number of life-stages

for hfreq = 1:length(phurset)
        Bdiag=blkdiag(m1,m2,m3,m4,m5,m6,m7);
        P = squeeze(bigcmat(hfreq,:,:));
        mu(hfreq)=max(abs(eig(Bdiag*kron(P,eye(S)))));
        [lamS1(hfreq),lam0(hfreq),W1(hfreq),W2(hfreq),b(:,:,hfreq),CV(:,:,hfreq),w0(:,hfreq),u0(:,hfreq),v0(:,hfreq)]=StocGrow_Markov1(K,S,P,B);
    end