Abstract: A new algorithm for computing geodetic coordinates from Cartesian coordinateswais is derived by solving a modification of Levin's quartic equation with Newton's iteration from two different initial values.Experiments with two given initial values show that the new algorithm is sufficiently precise,fast,and stable.