max(X1, X2, Y) :- X1 >= X2,
    Y=X1.
max(X1, X2, Y) :- X2 >= X1,
    Y=X2.


findall(Populacja, miasto(_, Populacja, _), Populacje),
foldl(max, Populacje, 0, maxPop),
miasto(Miasto, maxPop, _), Miasto.