A more streamlined approach than nested for loops to fill an array in R

Question

I\'m trying to write a function in R that models population sizes (of different age brackets) across a network of connected populations. As a simplified example, imagine I want