Iterative vector methods for computing geodetic latitude and height from rectangular coordinates

Two iterative vector methods for computing geodetic coordinates (phi, h) from rectangular coordinates (x, y, z) are presented. The methods are conceptually simple, work without modification at any latitude and are easy to program. Geodetic latitude and height can be calculated to acceptable precision in one iteration over the height range from -10(6) to +10(9) m.