Transformation of geocentric to geodetic coordinates without approximations

An exact and relatively simple analytical transform of the rectangular coordinates to the geodetic coordinates is presented. It does not involve any approximation and the accuracy of practical calculations depends exclusively on the round-off errors. The algorithm is based on one solution to the quartic equation in tg(45°-ψ/2), where ψ is the parametric (or eccentric) latitude.