function [lat,lon,h,slat,slon,sh,covl]=xyz2latlon(x,y,z,sx,sy,sz,sxy,syz,sxz,fast) in matlab

main.m
function [lat,lon,h,slat,slon,sh,covl]=xyz2latlon(x,y,z,sx,sy,sz,sxy,syz,sxz,fast)
% Convert Cartesian coordinates (x,y,z) to geodetic coordinates (latitude, longitude, and height)
% and calculates a 3x3 covariance matrix for the geodetic coordinates.
%
% Inputs:
%   x, y, z: Cartesian coordinates (meters)
%   sx, sy, sz: standard deviations of x, y, z (meters)
%   sxy, syz, sxz: covariances between x,y; y,z; x,z (meters^2)
%   fast (optional): if true, use a faster but less accurate algorithm
%
% Outputs:
%   lat, lon, h: geodetic coordinates (degrees, degrees, meters)
%   slat, slon, sh: standard deviations of lat, lon, and h (meters)
%   covl: 3x3 covariance matrix for lat, lon, h
693 chars
15 lines

main.m
% Define Cartesian coordinates and covariance
x = 100;
y = 200;
z = 300;
sx = 1;
sy = 2;
sz = 3;
sxy = 0.1;
syz = 0.2;
sxz = 0.3;

% Convert Cartesian coordinates to geodetic coordinates and covariance
[lat,lon,h,slat,slon,sh,covl] = xyz2latlon(x,y,z,sx,sy,sz,sxy,syz,sxz);

% Print results
disp(['Latitude: ' num2str(lat) ' +/- ' num2str(slat) ' deg']);
disp(['Longitude: ' num2str(lon) ' +/- ' num2str(slon) ' deg']);
disp(['Height: ' num2str(h) ' +/- ' num2str(sh) ' m']);
disp('Covariance matrix:');
disp(covl);
516 chars
21 lines