function resulting_vector = N_age_plus(popul,ageUP,rho)

last = numel(popul);
ages_plus = ageUP+1:last;
people = popul.*exp(-rho*(0:last-1));
resulting_vector = sum(people(ages_plus))/sum(people);