function y=esercizio(s,m,j,k);
A=leslie(s,m);
B=A/sum(m);
n=size(s);
x=ones(n+1,1);
y=zeros(k+1,1);
y(1)=x(j);
for i=1:k
   x=B*x;
   y(i+1)=x(j);
end
