cities = [
  -262 1240
  -1219 1375
  -623 1080

  1002 1624

  64 2309
  -552 1967
  757 4588

  970 1000
  1666 742

  -1539 -873
  -4139 -5062
  -4751 -7025

  1650 750
  1965 2257

  -2505 334
  -2873 375
  -3484 419

  -3223 1004
  -6661 725
  -7005 1092
  -10191 726
]

portals = [
  -313 949
  -600 1546
  880 1395

  -1875 450
  -1763 -865
  -6808 -6683

  -1819 521
  -6962 970
  -6447 618
]

scatter(cities(:,1),cities(:,2), [], 'b', 's', 'filled')
hold on
scatter(0, 0, [], 'g', 's', 'filled')
scatter(portals(:,1),portals(:,2), [], 'r', 's', 'filled')



axis('equal', 'ij')
grid on