convert x,y,z to latitude and longitude in matlab

main.m
% define constants for the WGS 84 ellipsoid
a = 6378137.0;
f = 1/298.257223563;

% define the x, y, z coordinates 
x = 1234567.89;
y = 2345678.90;
z = 3456789.01;

% calculate the longitude and latitude from the x, y, z coordinates
r = sqrt(x^2 + y^2 + z^2);
lon = atan2(y, x);
lat = atan2(z, sqrt(x^2 + y^2));

% convert the longitude and latitude from radians to degrees
lon = rad2deg(lon);
lat = rad2deg(lat);

% apply the WGS 84 ellipsoid correction to the latitude
e2 = f*(2-f);
N = a./sqrt(1-e2*sind(lat).^2);
alt = r./cosd(lat) - N;
lat = atan2d(z, sqrt(x^2 + y^2 - alt.^2));

% output the resulting longitude and latitude 
fprintf('Longitude: %f\nLatitude: %f\n', lon, lat);
683 chars
27 lines