function [geoid] = wgs84();

% WGS84 Datum
semimajor_axis = 6378137;        % equatorial radius
flattening     = 1/298.25792;
eccentricity   = sqrt(2*flattening - flattening^2);
geoid          = [semimajor_axis eccentricity];

