Solutions to the ellipsoidal Clairaut constant and the inverse geodetic problem by numerical integration