Hi Jose, thank you again for your response!

I will try to explain myself better. Having a list of cities and their populations, I have gotten a matrix of distances between them:

cities = Entity["City", #] & /@ {{"Istanbul", "Istanbul", 
     "Turkey"}, {"Moscow", "Moscow", "Russia"}, {"London", 
     "GreaterLondon", "UnitedKingdom"}, {"Berlin", "Berlin", 
     "Germany"}, {"Madrid", "Madrid", "Spain"}};

populations = QuantityMagnitude@EntityValue[cities, "Population"];

R = QuantityMagnitude@Table[GeoDistance[cities[[i]], cities[[j]]], {i, 1, 
     Length[cities]}, {j, 1, Length[cities]}];

What I want is to obtain a second matrix, where every element is given by the population contained within a circle, defined by center in city i, and radius R[[i,j]], excluding the populations of cities i and j. This is my own attempt, based on previous replies:

S = Table[Table[Boole[0 < R[[i, k]] < R[[i, j]]], {k, 1, 5}] . 
populations, {i, 1, 5}, {j, 1, 5}]

It may not be the most efficient solution, but it seems to work, at least for small cases. Any additional ideas would be much appreciated!

Thanks.

Thank you both @Jose M. Martin-Garcia and @Henrik Schachner, your replies have been very helpful!

I am trying now to construct a matrix S, where the element ij contains the total population inside the previous circle, excluding the ones of cities i and j.

I have found that the populations can be obtained by using:

population = CityData[#, "Population"] & /@ cities

But I am not sure how to construct that matrix. Any idea would be much appreciated.

Hello David,

I'd suggest to use entities for the cities:

cities = {
   Entity["City", {"Istanbul", "Istanbul", "Turkey"}], 
   Entity["City", {"Moscow", "Moscow", "Russia"}], 
   Entity["City", {"London", "GreaterLondon", "UnitedKingdom"}], 
   Entity["City", {"Berlin", "Berlin", "Germany"}], 
   Entity["City", {"Madrid", "Madrid", "Spain"}]
}

Define also the two cities of reference:

madrid = Entity["City", {"Madrid", "Madrid", "Spain"}]
berlin = Entity["City", {"Berlin", "Berlin", "Germany"}]

Compute the distance between them (units depend on whether you are, so in the US we get miles):

In[]:= dist = GeoDistance[madrid, berlin]
Out[]= Quantity[1162.07, "Miles"]

Now you can query:

In[]:= Select[cities, GeoDistance[madrid, #] < dist &]
Out[]= {
   Entity["City", {"London", "GreaterLondon", "UnitedKingdom"}], 
   Entity["City", {"Madrid", "Madrid", "Spain"}]
}

This performs an independent call to WA per city, which will be slow for long lists. We can do better in those cases by fetching all positions at once in advance.