Untitled

global R
paris = [48.8566, 2.3522]; % Paris
coordinates = [
    40.7128, -74.0060;  % New York
    -33.8688, 151.2093; % Sydney
    51.5074, -0.1278;   % London
    -23.5505, -46.6333; % San Paolo
    35.6762, 139.6503   % Tokyo
];

R = 6371; % Raggio della Terra in km

% Funzione per calcolare la distanza con la formula dell'emisenoverso


% Calcolo delle distanze da Parigi
num_cities = size(coordinates, 1);
distances = zeros(num_cities, 1);

for i = 1:num_cities
    distances(i) = calculate_distance(paris, coordinates(i, :));
end

% Identificazione della città più distante
[distance_from_paris, furthest_location_from_paris] = max(distances)

function distance = calculate_distance(coord1, coord2)
global R
    lat1 = deg2rad(coord1(1));
    lon1 = deg2rad(coord1(2));
    lat2 = deg2rad(coord2(1));
    lon2 = deg2rad(coord2(2));

    delta_lat = lat2 - lat1;
    delta_lon = lon2 - lon1;

    a = sin(delta_lat / 2)^2 + cos(lat1) * cos(lat2) * sin(delta_lon / 2)^2;
    c = 2 * atan2(sqrt(a), sqrt(1 - a));

    distance = R * c; % Distanza in km
end