function [w lambda]=ll(s,m,j,k)
%esercizio del foglio template1.pdf
%costruisce la matrice di Lotka e Leslie e ci applica il metodo delle potenze
%ritorna un vettore con le componenti j-esime dei vettori ottenuti

n=length(s);
A=zeros(n+1,n+1);

for j=1:n+1
    A(1,j)=m(j);
end
for i=2:n+1
    A(i,i-1)=s(i-1);
end

gamma=sum(m);
B=1/gamma*A;

e=ones(n+1,1);
v=e/norm(e,inf);

w=zeros(k,1);

for h=1:k
   v=B*v;
   v=v/norm(v,inf);
   w(h)=v(j);
end

u=B*v;
lambda=u(j)/v(j);

%eig(B)
